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General Relativity receives quantum corrections relevant at macroscopic distance scales 
and near event horizons. These arise from the conformal scalar degrees of freedom in the ex- 
tended effective field theory of gravity generated by the trace anomaly of massless quantum 
fields in curved space. The origin of these conformal scalar degrees of freedom as massless 
poles in two-particle intermediate states of anomalous amplitudes in flat space is exposed. 
These are non-local quantum pair correlated states, not present in the classical theory. At 
event horizons the conformal anomaly scalar degrees of freedom can have macroscopically 
large effects on the geometry, potentially removing the classical event horizon of black hole 
and cosmological spacetimes, replacing them with a quantum boundary layer where the effec- 
tive value of the gravitational vacuum energy density can change. In the effective theory, the 
cosmological term becomes a dynamical condensate, whose value depends upon boundary 
conditions near the horizon. In the conformal phase where the anomaly induced fluctuations 
dominate, and the condensate dissolves, the effective cosmological "constant" is a running 
coupling which has an infrared stable fixed point at zero. By taking a positive value in 
the interior of a fully collapsed star, the effective cosmological term removes any singular- 
ity, replacing it with a smooth dark energy interior. The resulting gravitational condensate 
star configuration resolves all black hole paradoxes, and provides a testable alternative to 
black holes as the final state of complete gravitational collapse. The observed dark energy 
of our universe likewise may be a macroscopic finite size effect whose value depends not on 
microphysics but on the cosmological horizon scale. The physical arguments and detailed 
calculations involving the trace anomaly effective action, auxiliary scalar fields and stress 
tensor in various situations and backgrounds supporting this hypothesis are reviewed. Origi- 
nally delivered as a series of lectures at the Krakow School, the paper is pedagogical in style, 
and wide ranging in scope, collecting and presenting a broad spectrum of results on black 
holes, the trace anomaly, and quantum effects in cosmology. 



I. INTRODUCTION: GRAVITATION AND QUANTUM THEORY 



Although it has been clear for nearly a century that quantum principles govern the microscopic 
domain of atomic, nuclear and particle physics, and certainly the Standard Model of Strong & Elec- 
troweak interactions is a fully quantum theory of matter, gravitational phenomena are still treated 
as completely classical in Einstein's General Relativity. Perhaps just as significant as this formal gap 
between gravitation and quantum principles, most of our intuitions about gravity remain essentially 
classical, particularly in the macroscopic domain. Although it is generally agreed that at the funda- 
mental microscopic Planck scale, a theory of gravitational interactions must come to terms with the 
quantum aspects of matter and even of spacetime itself, it is usually assumed that quantum effects 
are negligible on the scale of macroscopic phenomena, at astrophysical or cosmological distance scales, 
where classical General Relativity (GR) is presumed to hold full sway. 

In non-gravitational physics quantum effects are present on a wide range of scales in a variety of 
ways, some of them striking, others more subtle and less immediately appreciated. Semi-conductors, 
superfluids, superconductors and atomic Bose-Einstein condensates are unmistakable macroscopic mani- 
festations of an underlying quantum world. On astrophysical scales, the degeneracy pressure of fermions, 
which at first seemed an esoteric feature of quantum statistics is now fully accepted as the basis for the 
stability of such macroscopic objects as white dwarves and neutron stars, both of which are ubiquitous 
throughout the universe. As a posteriori consequences of quantum statistics one may note the periodic 
table, the foundations of chemistry itself and hence of biological processes, which being familiar in ordi- 
nary experience seem far less exotic than neutron stars or superfluidity. However chemical bonding, the 
structure and function of hemoglobin and DNA in the human body, and the overall stability of matter 
itself at ordinary temperatures and densities are every bit a consequence of quantum principles as a 
sample of superfluid 4 He climbing up the walls of its dewar. 

It was not only in the microscopic world of the atom but experiments on macroscopic matter and the 
puzzles they generated for classical mechanics, such as the ultraviolet problem of blackbody radiation 
and the specific heat of solids, that led to the development of quantum mechanics [HE]- Since quantum 
effects play a role in the properties of bulk matter and macroscopic phenomena in most every other 
area of physics, there is no reason why gravity, which couples to the energy content of quantum matter 
at all scales, should be immune from quantum effects on macroscopic scales. 

If the effects and predictions of a quantum theory of gravity can be tested only at Planck lengths 
or energies, the quest for such a theory would be mostly academic, an exercise better left for future 
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generations possessing more complete and accurate information about the ultra microscopic world. Thus 
the important questions at the outset are: Are quantum effects anywhere relevant or distinguishable in 
the macroscopic domain of gravitational phenomena? And if so, can one say anything reliable about 
macroscopic quantum effects in gravity, without necessarily possessing a complete, fundamental and 
well tested theory valid down to the microscopic Planck scale? 

There are two problems at the forefront of current research where there are indications that quantum 
effects may play a decisive role in gravitational physics at macroscopic distance scales. The first of these 
is that of ultimate gravitational collapse, presumed in classical GR to lead a singularity of spacetime 
called a black hole, which generates a number of theoretical paradoxes and challenges for quantum 
theory. The second problem of great current interest is the apparent existence of cosmological dark 
energy, which is causing the expansion of the universe to accelerate, and which has the same equation 
of state p = — p as that of the quantum vacuum itself. 

Both of these problems are at the intersection of gravitation and quantum theory at macroscopic 
length scales. In both cases an approximate or effective theory of quantum effects in gravity far from 
the Planck scale should be the appropriate and only framework necessary. In these notes, such an 
effective field theory (EFT) framework is developed according to the same general principles now widely 
recognized in other areas of physics. The essential technical tools involve the use of relativistic quantum 
field theory in curved spacetime, and the key ingredient of our analysis is the conformal or trace 
anomaly of the stress-energy- momentum tensor T a b of massless fields, a quantum effect with low energy 
implications. 

Within this essentially semi-classical EFT framework, the principal qualitative result is that Ein- 
stein's General Relativity can and does receive quantum corrections from the effects of the trace anomaly 
which are significant and in certain circumstances may even dominate at macroscopic distance scales, 
much larger than the Planck scale. No assumptions about the extreme short distance degrees of free- 
dom or the precise nature of fundamental interactions at that scale will be used or needed. Instead our 
analysis will rely only upon the assumption that the Principle of Equivalence in the form of coordinate 
invariance of the effective action of a metric theory of gravity under smooth coordinate transformations 
applies in the quantum theory at sufficiently low energy scales far below the Planck scale. With this 
moderate theoretical input, and without invoking unknown or esoteric physics beyond the Standard 
Model, we shall investigate the macroscopic effects of the quantum conformal anomaly on gravitational 
systems at the astrophysical scale of the event horizon of the collapse of massive stars, and on the very 
largest Hubble scale of the visible universe itself. 
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II. THE CHALLENGE OF BLACK HOLES 



The first problem in which classical General Relativity is challenged by quantum theory is in the 
physics of black holes. Before discussing quantum effects, let us review the standard classical theory of 
the gravitational collapse of massive stars. 

A normal main sequence star sustains itself by nuclear fusion of hydrogen into helium. The star 
itself formed when enough predominantly hydrogen gas collapsed to a high enough density for nuclear 
fusion reactions to occur. The energy generated by the fusion of H into He nuclei generates heat and 
pressure, which supports the star against further gravitational collapse. A main sequence star remains 
in this stable steady state, producing radiant energy for typically several billions of years, depending 
upon its mass. Eventually, the hydrogen is exhausted, and the star goes through a sequence of less 
exothermic nuclear reactions, fusing nuclei of heavier and heavier elements to extract energy from the 
difference in rest masses. Since iron is the most stable nucleus, this process eventually exhausts all the 
available sources of nuclear fusion energy. At that point the star can no longer sustain itself against the 
force of gravity, and its matter must resume gravitational collapse upon itself. 

Classically, nothing can halt this collapse. However, quantum matter obeys quantum statistics. 
Because of Fermi-Dirac statistics, if the mass of the star is not too great, and it has cooled sufficiently, 
a new stable configuration, a white dwarf star held up by its quantum degeneracy pressure can be 
formed. In other cases, the collapse of the outer envelope onto the Fe core produces a violent explosion, 
a stellar nova or supernova, in which prodigious amounts of mass and energy are ejected. This leaves 
behind a even more compact object in which the electrons and protons are forced under high pressure 
to become neutrons. A neutron star, sustained against further collapse by the quantum degeneracy 
pressure of neutron matter, rotating very rapidly at nuclear densities and beaming out radiation guided 
by its strong magnetic fields may be observed by astronomers as a pulsar. 

If the mass of the stellar remnant core exceeds a certain value, called the Tolman-Oppenheimer- 
Volkoff (TOV) limit of 1.5M to 3.2M (depending upon the eq. of state of dense nuclear matter, 
which is not very accurately known), not even the neutron degeneracy pressure is enough to prevent 
final and inexorable collapse due to gravity [3]. Since we have no direct observations of these final 
stages of complete gravitational collapse, it is here that the reliance upon Einstein's theory of General 
Relativity becomes critical, and the discussion takes on a decidely more mathematical flavor. 
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A. Black Holes in Classical General Relativity 



Just a year after the publication of the field equations of General Relativity (GR), K. Schwarzschild 
found a simple, static, spherically symmetric solution of those equations, with the line element [3J, 

dr 2 

ds 2 = -f(r) dr 2 + — - + r 2 (dO 2 + sin 2 9 # 2 ) , (2.1) 
h[r) ' 

where in this case the two functions of r are equal: 

f{r) = h (r) = 1 - ^ . (2.2) 
c z r 

This Schwarzschild solution to the vacuum Einstein's equations, with vanishing Ricci tensor R a b = 
and stress tensor T a b = for all r > describes an isolated, non-rotating object of total mass M. In 
that sense it is the gravitational equivalent of the Coulomb solution, 

4> = - (2.3) 

r 

for the electrostatic potential of an isolated, static charge e in Maxwell's theory of electromagnetism. 
Just as in the Coulomb case, the Schwarzschild solution has a singularity at the origin of the spherical 
coordinates at r = 0, where the gauge invariant field strengths (measured in gravity by the Riemann 
curvature tensor and its contractions) diverge, and there is a delta function source. 

In classical electromagnetism at the finite scale of the classical electron radius, r c = e 2 /mc 2 , where 
the electrostatic self-energy becomes comparable to the rest mass energy, some deviation from the 
simple picture of a structureless point particle is to be expected. In quantum electromagnetism (QED) 



the classical linear divergence of (2.3) is softened somewhat into a logarithmic ultraviolet divergence 
of the self-energy of a charged Dirac particle. This logarithmically divergent self-energy is absorbed 
into a renormalization of its total observable mass. However already at the larger scale of the electron 
Compton wavelength h/mc > r c , the single electron description has to be replaced by the many body 
description of a quantum field theory with vacuum polarization effects. Hence the pointlike singularity 



and linear divergence of the classical Coulomb potential (2.3) is not present in the more accurate many 
body quantum theory. 



Apart from the singularity at r = 0, analogous to (2.3) the Schwarzschld line element also possesses 
another kind of mathematical singularity at the finite Schwarzschild radius, 



2GM „ / M 

^0 



c2 _ 2.953 — J km, (2.4) 
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where the function f(r) = h(r) vanishes. This macroscopic radius is the location of the Schwarzschild 
event horizon, the locus of points which defines the sphere dividing the exterior region from the interior 
region. It is the analog of the classical radius r c where the Newtonian gravitational self-energy GM 2 jr 
becomes comparable to the total rest mass energy Mc 2 . Thus r s is the length scale at which some 
substructure should be expected. 

Classical General Relativity does not give much hint of this substructure. Instead, the change of sign 



of the functions f(r), h(r) for r < r s in (2.2) indicates that for the interior region t becomes a spacelike 
variable, while r becomes timelike. Hence any radiation, even a light ray emanating from a point in the 
interior cannot propagate outward and is drawn inexorably toward the singularity at r = 0, giving rise 
to the popular name black hole. 



If the Schwarzschild solution (2.1 )-(2.2) for r < r s is taken seriously, the singularity at r = is present 
in Einstein's theory for any mass M > 0, including the certainly macroscopic mass of a collapsed star 
with the mass of the sun, Mq ~ 2 x 10 33 gm. or even that of supermassive objects with masses 10 6 
to 1O 9 M0. The collapse of such enormous quantities of matter with vastly more degrees of freedom 
than that of a single electron to a single mathematical point at r = certainly presents a challenge to 
the imagination, and one that it seems Einstein himself sought arguments to avoid [5]. The situation 
is scarcely more acceptable if the singularity is removed only by the intervention of quantum effects at 
the extremely tiny Planck length (G7i/c 3 )2 ~ 1.6 x 10~ 33 cm. 

A light wave emitted from any r > r s with local frequency cvi oc outward towards infinity is redshifted 
according to the redshift relation, 



1 T 

Woo = C^oc/ 5 = wzoci/l , (2.5) 

showing that a light wave emitted at the horizon becomes redshifted to zero frequency and cannot 
propagate outward at all. Conversely and equivalently, a light wave with the finite frequency oj = ojoo 
far from the black hole is blueshifted to an infinite local frequency at the horizon. This gravitational 
redshift /blueshift is purely a kinematic consequence of the classical time dilation effect of a gravitational 
field, which has been tested in a number of experiments [6-8j. The event horizon is therefore a kind of 
critical surface for the propagation of light rays, and hence all other matter interactions. 

Unlike the central singularity at r = 0, the scalar invariant quantities that can be constructed from 
the contractions of the Riemann curvature tensor remain finite as r — > r s . Thus the fully contracted 
quadratic Riemann invariant 

12r 

r 



nabcd n S fo R\ 



which diverges at the origin remains finite at r 



Moreover, although the time for an infalling 



particle to reach the horizon is infinite for any observer remaining fixed outside the horizon, the proper 
time measured by the particle itself during its infall remains finite as r — > r s [3 [7] . Thus despite the 
singularity of the Schwarzschild coordinates at r = r s , physics must continue onto smaller values of r in 



the interior region. Since the line element (2.1) is again non-singular for < r < r s , and in the absence 



of clear evidence to the contrary, the most straightforward possibility would seem to be to assume that 
this non-singular vacuum interior (up to the origin or at least some extreme microscopic scale much less 
than r s ) can be matched smoothly to the non-singular exterior Schwarzschild solution. 

This matching was achieved by the coordinate transformations and analytic continuation of the 
Schwarzschild solution found by Kruskal and Szekeres [9j. The Kruskal maximal analytic extension of 
the Schwarzschild geometry is pictured in the Carter-Penrose conformal diagram, Fig. [TJ 

r = 




r = Q 

FIG. 1: The Carter-Penrose conformal diagram of the maximal Kruskal analytic extension of the Schwarzschild 
geometry. Radial light rays are represented in this diagram as 45° lines. The angular coordinates 0,<j) are 
suppressed. 



The analytic extension of the Schwarzschild geometry relies on finding a judicious change of the 



time and radial coordinates (t,r) of (2.1) to new ones (T, X), which are regular on the horizon, and 
therefore can be used to describe the local geometry there without singularities. Explicitly, for r > r s 
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this coordinate transformation to Kruskal-Szekeres coordinates (T, X) is given by 



The inverse transformation is 



i 

-^-ijV^sinhQ^) (2Ja) 
X=(f -l)%^co S h(A) . (2.7b) 



1 ) e r / r s = X 2 — T 2 . (2.8b) 



In the new (T,X) coordinates, the Schwarzschild line element (2.2) becomes 



4r 3 , 

ds 2 = — ^ e r/2r s (-dT 2 + dX 2 ) + rW , (2.9) 
r 



where r is to be regarded as the function of (T, X) given implicitly by d2.8bh, and dVL 2 = dO 2 + sin dm 
is the usual spherical line element on § 2 . 

In these Kruskal-Szekeres coordinates, the event horizon at r = r s is seen actually to be comprised 
of two distinct null surfaces, the future and past horizons for T = ±X, respectively. The entire exterior 
region r > r s is mapped into the region X > |T| > of the (X,T) plane. The X = +T boundary of 
this region gives the coordinates of a particle trajectory infalling into the black hole, a depicted in Fig. 
[TJ whereas the X = —T boundary corresponds to the time reversed coordinates of a particle trajectory 
outgoing from a white hole. The time reversed case must be present mathematically because of the 



second order nature of Einstein's eqs., and the static solution (2.1) which admits the time reversal 
symmetry t — > — t and T — > —T. Whether or not this time reversed white hole case corresponds to any 
real macroscopic body is of course another question. 



Since the Schwarzschild line element in Kruskal-Szekeres coordinates (2.9) is completely regular 
at T = ±X, one can equally well consider the extension of the coordinates to the interior regions, 
T > X > and T < —X < 0, into the black hole and white hole interior regions at the top and bottom 
of Fig. [TJ at least as far and until the true curvature singularity at r = is reached. Inspection of (2.8) 
shows that this implies an analytic continuation of the original Schwarzschild (t, r) coordinates around 
a logarithmic branch cut to complex values. The relation of the (T, X) coordinates to the original (t, r) 
coordinates is singular at r = r s , although the (T, X) coordinates themselves are regular at T = ±X. 
Similarly, by a further complex analytic continuation, one can continue to the parallel exterior region 
on the left of Fig. [TJ with X < —\T\ < 0. Thus, by the relatively simple but singular change of 



coordinates (2.7)-(2.8), we seem to have reached the conclusion that the simplest static spherically 
symmetric Schwarzschild solution to the vacuum Einstein's equations predicts the existence not only of 
a true singularity at r = but also of an entirely separate and macroscopically large, asymptotically 
flat region in addition to the original one. 

Of course, one is free to assert that the white hole interior region and the parallel asymptotically 
flat universe do not exist, and excise them in a gravitational collapse from realistic initial conditions, 
replacing the excised region of Fig. [I] with the non-singular interior of a collapsing matter distribution. 
However, this is a matter of initial conditions, and nothing in the equations of General Relativity 
themselves forces us to do this. Dirac expressed skepticism of the interior Schwarzschild solution on 
physical grounds [TO]. 

The apparently unphysical features of the Schwarzschild solution appear as soon as we admit complex 
analytic continuation of singular coordinate transformations. Based on the Principle of Equivalence 
between gravitational and inertial mass, Einstein's theory possesses general coordinate invariance under 
all regular and real transformations of coordinates. It is the appending to classical General Relativity of 
the much stronger mathematical hypothesis of complex analytic continuation through singular coordinate 
transformations that leads to the global aspects of the Schwarzschild solution which may be unrealized 
in Nature. 

The point which is often left unstated is that the mathematical procedure of analytic continuation 
through the null hypersurface of an event horizon actually involves a physical assumption, namely 
that the stress-energy tensor T a b is vanishing there. Even in the purely classical theory of General 
Relativity, the hyperbolic character of Einstein's equations allows generically for stress-energy sources 
and hence metric discontinuities on the horizon which would violate this assumption. Additional physical 
information is necessary to determine what happens as the event horizon is approached, and the correct 
matching of interior to exterior geometry. What actually happens at the horizon is a matter of this 
correct physics, which may or may not be consistent with complex analytic continuation of coordinates 



(2_7)-(|2_8). 

The static Schwarzschild solution of an isolated uncharged mass was generalized to include electric 
charge by Reissner & Nordstr0m and more interestingly for astrophysically realistic collapsed 
stars, to include rotation and angular momentum by Kerr [12J. The complete analytic extensions of the 
Reissner-Nordstr0m and Kerr solutions were found as well |13j . The global properties of these analytic 
extensions are more complicated and arguably even more unphysical than in the Schwarzschild case. For 
slowly rotating black holes with angular momentum J < GM 2 /c, there are an infinite number of black 
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hole interior and asymptotically flat exterior regions, and closed timelike curves in the interior region(s), 
which violate causality on macroscopic distance scales [Hj. Again these apprently unphysical features 
appear in GR only if the mathematical hypothesis of complex analytic extension and continuation 
through real coordinate singularities are assumed. This analytic continuation is generally invalid if 
there are stress-tensor sources encountered at or before the breakdown of coordinates. 

The Schwarzschild white hole and the analytic extension through the horizons also raises questions 
about macroscopic time reversibility. Once classical particles fall through the future event horizon, 
there is no way to retrieve them without violating causality, and something irreversible would seem to 
have occurred. This is somewhat troubling from the point of view of thermodynamics, since if matter 
disappears completely from view when it falls into a black hole, it carries any entropy it has in its internal 
states with it, and the entropy of the visible universe would apparently have decreased, violating the 
second law of thermodynamics, 

A5>0, (2.10) 

which states that the total entropy of an isolated system must be a non-decreasing function of time in 
any spontaneous process. 

At the same time the infall of matter into a black hole certainly increases its total energy. In the 
general black hole solution characterized by mass M, angular momentum J, and electric charge Q, one 
can define a quantity called the irreducible mass Mi rr by the relation, 

/ Q2 \2 2t2 

M2 = Mirr + J^) +77^72", (2.H) 



or 

M *r = f" " W + h ~ ~ ° 2 1 ' (2 ' 12) 

and show that in any classical process the irreducible mass can never decrease [15] : 

AM? r >0. (2.13) 

Since one can also show that the irreducible mass is related to the geometrical area A of the event 
horizon of the general Reissner-Nordstr0m-Kerr-Newman black hole via 

A= l ^-Ml rl (2.14) 

this theorem is equivalent to the statement that the horizon area is a non- decreasing function of time 
in any classical process \15\ [T6] . 
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Simply by taking the differential form of (2.11) one obtains |17) 

„2 



dE = dMc 2 



IttG 



KdA + VtdJ + ®dQ 



(2.15) 



which is just the differential form of Smarr's formula for a Kerr-Newman rotating, electrically charged 
black hole, in which 



9. 
A 



1 

M 
AttJ 

Wa ' 
Q_ 

M 
4vrG 



c 

AG 



c 

2G 



Att 2 G 
IfX 2 



2nQ 2 
~A~J~ 



(Q 4 + 4c 2 J 2 ) 



2GM 2 -Q 2 + 2y / G' 2 M 4 - GM 2 Q 2 - c 2 J 2 



(2.16a) 
(2.16b) 
(2.16c) 
(2.16d) 



are the horizon surface gravity, angular velocity, electrostatic potential and area respectively \17\ I18j. 



All dimensionful constants have been retained to emphasize that (2.15)-(2.16) are formulae derived from 
classical GR in which no h whatsoever appears. Notice also that the coefficient of dA in (2.15 ), c 2 k/8itG 
has both the form and dimensions of a surface tension. 



The classical conservation of energy is expressed by the first law of black hole mechanics (2.15). 



The classical area theorem (2.13) naturally evokes a connection to entropy and the second law of 



thermodynamics (2.10). //horizon area (or more generally any monotonic function of it) could somehow 
be identified with entropy, and this entropy gain is greater than the entropy lost by matter or radiation 



falling into the hole, then the second law (2.10) would remain valid for the total or generalized entropy 
of matter plus black hole horizon area. The simplest possibility would seem to be if entropy is just 
proportional to area. 

Motivated by these considerations, and as suggested by a series of thought experiments, Bekenstein 



proposed that the area of the horizon (2.14) should be proportional to the entropy of a black hole 



[19] . Since A does not have the units of entropy, it is necessary to divide the area by another quantity 
with units of length squared before multiplying by Boltzmann's constant k g , to obtain an entropy. 
However, classical General Relativity (without a cosmological term) contains no such quantity, G/c 2 
being simply a conversion factor between mass and distance. Hence Bekenstein found it necessary 
for purely dimensional reasons to introduce Planck's constant h into the discussion. Then there is a 
standard unit of length, namely the Planck length, 



L 



IHG 



pi 



1.616 x 10" 33 cm. 



(2.17) 
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S BH =-yk B -^, (2.18) 
1^1 



Bekenstein proposed that the entropy of a black hole should be 

A 

IT 
J Pl 

with 7 a constant of order unity [19]. He showed that i/such an entropy were assigned to a black hole, 
so that it is added to the entropy of matter, Stot = S m + S BH , then this total generalized entropy would 
plausibly always increase. In fact, this is not difficult at all, and the generalized second law AStot > is 
usually satisfied by a very wide margin, simply because the Planck length is so tiny, and the macroscopic 
area of a black hole measured in Planck units is so enormous. Hence even the small increase of mass and 
area caused by dropping into the black hole a modest amount of matter and concomitant loss of matter 
entropy AS m < is easily overwhelmed by a great increase in S BH , AS BH S> |A5 m |, guaranteeing that 
the generalized total entropy increases: AStot > 0. 

Since h has entered the assignment of entropy to a black hole horizon, the discussion can no longer be 
continued in purely classical terms, and we must discuss quantum effects in black hole geometries next. It 
is worth remarking that whereas Planck's constant enters the thermodynamics of macroscopic quantum 
systems, such as the formulae for black body radiation, by normalizing the volume of the integral over 



phase space, no such interpretation is available for (2.18). Instead h has been used to form a new quantity 



Lpi with units of length, which we would ordinarily associate with the microscopic length scale at which 
strong quantum corrections to Einstein's theory should become important. Why such a microscopic, 
fundamentally quantum length scale should be needed to determine the bulk thermodynamic entropy 
of a macroscopic object on the scale of kilometers where classical GR applies, and no large quantum 
corrections to GR are expected near the horizon, is far from clear. 

B. Quantum Black Holes and Their Paradoxes 

Since classically all matter is irretrievably drawn into a black hole, the idea that black holes can 
instead radiate energy seems quite counterintuitive. More remarkable still is Hawking's argument that 
this radiation would necessarily be thermal radiation, with a temperature [20] 



= _hK_ j=q= m 

H 2nck B %irGk B M' 1 ' 

where the first equality is general, and the second equality applies only for a Schwarzschild black hole 
with J = Q = 0. With the temperature inversely proportional to its mass assigned to a black hole by 
this formula, if we assume that the first law of thermodynamics in the form 

dE = dMc 2 = T H dS BH + QdJ + <Z>dQ. (2.20) 
12 



applies to black holes, then the coefficient 7 in (2.18) is fixed to be 1/4. This formula is simple and 
appealing, and has been generally accepted since soon after Hawking's paper first appeared. How- 
ever, simultaneously and from the very beginning, a number of problems with this thermodynamic 
interpretation made their appearance as well. 



The first curious feature of (2.20) is that h cancels out between T and dS BH . Of course, this is a 



necessary consequence of the fact that (2.20) is identical to the classical Smarr formula (2.15) in which 
h does not appear at all. The identification of S BH with the entropy of a black hole is founded on the 



purely classical dynamics of Christodoulou's area law (2.13), in which quantum mechanics played no 
part whatsoever. Clearly, multiplying and dividing by h does not necessarily make a classical relation 



a valid one in the quantum theory. On the other hand, the Hawking temperature (2.19) is a quantum 
spontaneous emission effect, analogous to the Schwinger pair creation effect in a strong electric field 
|21j which vanishes in the classical limit h — >■ 0. Temperature, usually accepted as a classical concept 
has no apparent meaning for a black hole in the strictly classical limit, unless it is identically zero. 
As a consequence, if the identification of the classical area rescaled by k B L~j^ with entropy and the 



thermodynamic interpretation of (2.20) is to be generally valid in the quantum theory, then the classical 
limit h — > (with M fixed) which yields an arbitrarily low Hawking temperature, assigns to the black 
hole an arbitrarily large entropy, completely unlike the zero temperature limit of any other cold quantum 
system. 

Closely related to this paradoxical result is the fact, pointed out by Hawking himself [22], that a 
temperature inversely proportional to the M = E/c 2 implies that the heat capacity of a Schwarzschild 
black hole, 

dE 8vrG/c D M 2 Me 2 



dT H he T H 



< 



(2.21) 



is negative. In statistical mechanics the heat capacity of any system (at constant volume) is related to 
the energy fluctuations about its mean value {E) by 



1 



k D T 2 



(E-(E)Y)>0 



(2.22) 



If pressure or some other thermodynamic variable is held fixed there is an analogous formula. Hence on 
general grounds of quantum statistical mechanics, the heat capacity of any system in (stable) equilibrium 



must be positive. The positivity of the statistical average in (2.22) requires only the existence of a well 
defined stable ground state upon which the thermal equilibrium ensemble is defined, but is otherwise 
independent of the details of the system or its interactions. 
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In fact, it is easy to see that a black hole in thermal equilibrium with a heat bath of radiation at 
its own Hawking temperature T = T H cannot be stable |22j. For if by a small thermal fluctuation it 
should absorb slightly more radiation in a short time interval than it emits, its mass would increase, 



AM > 0, and hence from (2.19) its temperature would decrease, AT < 0, so that it would now be 
cooler than its surroundings and be favored to absorb more energy from the heat bath than it emits in 
the next time step, decreasing its temperature further and driving it further from equilibrium. In this 
way a runaway process of the black hole growing to absorb all of the surrounding radiation in the heat 
bath would ensue. Likewise, if the original fluctuation has AM < 0, the temperature of the black hole 
would increase, AT > 0, so that it would now be hotter than its surroundings and favored to emit more 
energy than it absorbs from the heat bath in the next time step, increasing its temperature further. 
Then a runaway process toward hotter and hotter evaporation of all its mass to its surroundings would 
take place. In either case, the initial equilibrium is clearly unstable, and hence cannot be a candidate 
for the quantum ground state for the system. This is the physical reason why the positivity property 



of (2.22) is violated by the Hawking temperature (2.19). The instability has also been found from the 



negative eigenvalue of the fluctuation spectrum of a black hole in a box of (large enough) finite volume 



The time scale for this unstable runaway process to grow exponentially is the time scale for fluctu- 
ations away from the mean value of the Hawking flux, not the much longer time scale associated with 
the lifetime of the hole under continuous emission of that flux. This time scale for thermal fluctuations 
is easily estimated. It is the typical time between emissions of a single quantum with typical energy (at 
infinity) of k B T, of a source whose energy emission per unit area per unit time is of order (k B T H ) 4 /fr?c 2 . 
Multiplying by the area of the hole A ~ (GM) 2 /c 4 , and dividing by the typical energy k B T, we find 
the average number of quanta emitted per unit time. The inverse of this, namely 

1 (he) 3 r„ tl ( M \ , 

At ~ - v ' q ~ ~ 1(T 5 — sec (2.23 

is the typical time interval (as measured by a distant observer) between successive emissions of individual 
Hawking quanta (again as observed far from the black hole). This time scale is quite short: 10 //sec 
for a solar mass black hole. Any tendency for the system to become unstable would be expected to 
show up on this short a time scale, governing the fluctuations in the mean flux, which is of order of 
the collapse time itself and before a steady state flux could even be established. With the existence of 
a stable equilibrium in doubt, one may well question whether macroscopic equilibrium thermodynamic 
concepts such as temperature or entropy are applicable to black holes at all. 
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Another rather peculiar feature of the formula (2.18) for the entropy is that it is non- extensive, 
growing not like the volume of the system but its area. It is non-extensive in a second important 



respect. The fixing of 7 = 1/4 by (2.19) and (2.20) is independent of the number or kind of particle 
species. Normally, we would expect the entropy of a system to grow linearly with the number of distinct 
particle species it contains. For example, if the number of light neutrino species in the universe were 
doubled, the entropy in the primordial plasma would be doubled as well, because the available states 
of one species are independent of and orthogonal to the states of a second distinct species, and must 



be counted separately. The formula (2.18) with only fundamental constants and the pure number 
coefficient 7 = 1/4 does not seem to allow for this, 't Hooft found a way to compensate the number 
of species factor in the horizon "atmosphere," but this configuration with matter or radiation densely 
concentrated near r = r = r s is no longer then a Schwarzschild black hole [24J . It also remains singular 
at the origin. 

Finally, let us state the obvious: a solution of any set of classical field equations is simply one 
particular configuration in the space of field configurations. As such, one would not usually associate 
any entropy with it. Matter sources to the field equations which have internal degrees of freedom may 
carry entropy of course, but a vacuum solution to the equations such as the Schwarzschild solution 



(2.1)-(2.2), with at most a singular point source at the origin would not ordinarily be expected to carry 
any entropy whatsoever. What "entropy" should one associate with the structureless classical Coulomb 
field ggfr? 

In this connection it may be worth pointing out that although the Coulomb field does not have an 
event horizon, there would still be an "information paradox" if we allowed charged matter to be attracted 



into (or emitted from) the Coulomb field singularity in (2.3) at r = and disappear from (or appear 
in) the visible universe. Such disappearance or appearance processes would clearly violate unitarity as 
well. In quantum theory we exclude such a possibility by restricting the Hilbert space of states to those 
with wave functions that are normalizable at the origin, so that no energy or momentum flux can either 
vanish or appear at the Coulomb singularity. Note that we impose this boundary condition at r = 
without any detailed knowledge of the extreme short distance structure of a charged particle in QED, 
confident that whatever it is, unitarity must be respected. In the case of black hole radiance by contrast, 
the positive Hawking energy flux at infinity must be balanced by a compensating flux down the hole 
and eventually into the singularity at r = 0, which makes clear why problems with unitarity and loss of 
information must result from such a flux. Boundary conditions on the horizon have consequences for the 
behavior of the fluctuations at both the singularity and at infinity. As even the example of scattering off 
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the Coulomb potential shows, these boundary conditions require physical input. Unphysical boundary 
conditions can easily lead to unphysical behavior. 



Returning to the entropy (2.18), it is instructive to evaluate S BH for typical astrophysical black holes. 
Taking again as our unit of mass the mass of the sun, Mq ~ 2 x 10 33 gm., we have 

2 



S BH ~ 1.050 x W 77 k E 



M 



(2.24) 



This is truly an enormous entropy. For comparison, we may estimate the entropy of the sun as it is, 
a hydrogen burning main sequence star, whose entropy is given to good accuracy by the entropy of a 
non-relativistic perfect fluid. This is of the order Nk B where N is the number of nucleons in the sun 
N ~ M & /m N ~ 10 57 , times a logarithmic function of the density and temperature profile which may 
be estimated to be of the order of 20 for the sun. Hence the entropy of the sun is roughly 



5 ~ 2 x 10 58 k B , 



(2.25) 



or nearly 19 orders of magnitude smaller than (2.24). 



A simple scaling argument that the entropy of any gravitationally bound object with the mass of 



the sun cannot be much more than (2.25) can be made as follows. The entropy of a relativistic gas 



at temperature T in equilibrium in a box of volume V is of order VT^ . The total energy is of order 
VT 4 . Eliminating T from these relations gives S ~ V^E*. For a relativistic bound system the energy 
E ~ .Mc 2 while the volume is of order r 3 ~ (GM/c 2 ) 3 . Hence S ~ (VGM)2. Keeping track of k B , K 



and c in this estimate gives 



(2.26) 



where M Pi = ^hc/G = 2.176 x 10~ 5 gm. If there are v species of relativistic particles in the object 
then this estimate should be multiplied by v±. This estimate applies to the entropy of relativistic 



radiation within the body, and is lower than (2.25) because the radiation pressure in the sun is small 



compared to the non-relativistic fluid pressure. However, the entropy from the relativistic radiation 



pressure (2.26) grows with the 3/2 power of the mass, whereas the non-relativistic fluid entropy (2.25) 
grows only linearly with M . For stars with masses greater than about 50M Q which are hot enough for 



their pressure to be dominated by the photons' T 4 radiation pressure, (2.26) indeed gives the correct 
order of magnitude estimate of such a star's entropy at a few times 10 59 fc s [25]. In order to have such 
an entropy, the temperature of the star must be of order E/S ~ 10 MeV or 10 11 °K, while the Hawking 
temperature (2.19) for the same 5OM black hole is a very, very cold 10~ 9 °K. It is difficult to see how 
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the entropy of the black hole could be a factor of 10 20 larger while its temperature is a factor of 10 20 
lower than a relativistic star of the same mass. 



The point is that even for the extreme relativistic fluid the entropy (2.26) for a gravitationally bound 



system in thermal equilibrium (in which entropy is always maximized) grows only like the 3/2 power of 
the mass, and hence will always be much less than (2.24), proportional to M 2 for very massive objects. 



Moreover the discrepant factor between (2^4) and (2.26) is of order (M/M P {)* ~ 1O 19 (M/M 



i 

) 2 , no 



matter what the non-black hole progenitor of the black hole is. Since the formula (2.24) makes no 



reference to how the black hole was formed, and a black hole may always be theoretically idealized as 



forming from an adiabatic collapse process, which keeps the entropy constant, (2.24) states that this 
entropy must suddenly jump by a factor of order 10 19 for a solar mass black hole at the instant the 
horizon forms. When Boltzmann's formula, 



S = k B InW(E) 



(2.27) 



is recalled, relating the entropy to the total number of microstates in the system W(E) at the fixed 



energy E, we see that the number of such microstates of a black hole satisfying (2.24) must jump by 
exp(10 19 ) at that instant at which the event horizon is reached, a truly staggering proposition. 

This tremendous mismatch between the number of microstates of a black hole inferred from S BH 
and that of any conceivable physical non-black hole progenitor is one form of the information paradox. 
Another form of the paradox is that since black holes are supposed to radiate thermally at temperature 
T H up until their very last stages, when their mass falls to a value of order Mpi, there would seem to 
be no way to recover all the information apparently lost in the black hole formation and evaporation 
process. This difficulty with S BH is so severe that it led Hawking to speculate that perhaps even the 
quantum mechanical unitary law of evolution of pure states into pure states would have to be violated 
by black hole physics [26J. Although this speculation has currently fallen into disfavor [27J, it is still far 
from clear what the missing microstates of an uncharged black hole are, and how exactly unitarity can 



be preserved in the Hawking evaporation process if (2.19) and (2.24) are correct. 



For all of these reasons the thermodynamic interpretation of (2.20) remains problematic in quantum 



theory. On the other hand, if the h is cancelled and one simply returns to the differential Smarr relation 



(2.15), derived from classical GR, these difficulties immediately vanish. One would only be left to 



explain the relationship of surface gravity k to surface tension. In Refs. [28, 29 1 and Sec. VII a possible 



resolution is proposed in which area is not entropy at all but indeed the area of a physical surface and 
the surface gravity can be related to the surface tension of this surface. 
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Although the thermodynamic interpretation of (2.20) and (2.24), leads to myriad difficulties, it has 
been essentially universally assumed, and great ingenuity has been devoted to postulating the new 
physics of some kind which would be needed to account for the vast number of microstates in the 



interior of a black hole required by (2.24). It is not possible to do justice to the various approaches in 
detail here. Curiously the multitude of approaches all seem to give the same answer, despite the fact 
that the states they are counting are very different [3Qj. Of course, any theory that reduces to Einstein's 
theory in an appropriate limit will conserve energy and obey the first law of thermodynamics. If the 
effective action of the theory reduces to the Einstein- Hilbert action then the logarithm of its generating 



functional Z should produce the entropy (2.24), suggesting that these formal results are actually a 



feature of the classical theory, embodied in the Smarr formula (2.15), and independent of the model 
used or microstates counted. 

The lines of research involving exotic internal constituents to obtain S BH are all the more remarkable 
when one recalls where we began, with the classical GR expectation that all one has to do is change 



coordinates as in (2.7) to see that "nothing happens" at the event horizon to a particle falling through. 
If the horizon is really just a harmless coordinate singularity-the very assumption underlying the argu- 



ments leading to the Hawking temperature, and hence the entropy (2.24), how can the semi-classical 
assumption of no energy or stresses, and analyticity and regularity at the geometry at the horizon with 
no substructure in the interior then lead to the diametrically opposite conclusion of exotic new physics, 
with exp(10 19 ) additional microstates, and perhaps the radical alteration of classical spacetime itself 
the instant the horizon is reached? 

Counting the microstates hiding in electrically charged black holes in string theory or other models 
also leave unanswered the question of how a presumably well-defined quantum theory with a stable 
ground state (which always has a positive heat capacity) could ever yield the negative heat capacity 



(2.21) of the original, uncharged Schwarzschild black hole in 3 + 1 dimensions. 



Another line of thought has attempted to identify the entropy (2.24) with quantum entanglement 



entropy |31H33j . This is the entropy that results when a quantum system is divided into two spatial 
partitions and one sums over the microstates of one of the partitions, forming a mixed state density 
matrix from a pure state wave functional even at zero temperature. This has the attractive feature that 
on general grounds it is proportional to the area of the surface at which the two partitions are in contact. 
It has the unattractive feature that the coefficient of the area law in quantum field theory is infinite. 
The reason for this ultraviolet divergence is the same as the reason for the area law itself, namely that 
the largest contribution to the entanglement entropy comes from the ultraviolet components of the wave 
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functional within a vanishingly thin layer near the surface. If this divergence is cut off at a length scale 



of the order of Lpi, a large but finite entropy of order of (2.18) is obtained |32j . 

Counting these states at very high frequencies as contributing to the entropy is sensible only if those 
states are occupied. In quantum statistical mechanics the unoccupied states at arbitrarily high energies, 
no matter how many of them, do not contribute to the thermodynamic entropy of the system. In the 
black hole case, the standard classical assumption of the horizon as a harmless coordinate singularity 
and the Hawking-Unruh state corresponding to this classical assumption (discussed in more detail in 
the next section) treats these high frequency states as unoccupied vacuum states with respect to a 



locally regular coordinate freely falling system at the horizon, such as (2.7). Hence it is difficult to see 
why they should be counted as contributing to the entropy. If on the contrary one treats these modes 



as occupied with respect to the singular static frame (2.1), then it is difficult to see why their mean 
energy-momentum or fluctuations in (T a b ) should be negligible near the horizon. Since the Hawking 
thermal flux originates as radiation closer and closer to r s with arbitrarily high frequencies at late times, 
if these states are occupied it is also difficult to see why any cutoff at the Planck scale or otherwise 
should be imposed to compute the entanglement entropy. 

The attempt to count microstates near the horizon to account for the black hole entropy associated 
with the Hawking effect brings us face to face with questions about the structure and meaning of the 
"vacuum" itself at trans- Planckian frequencies. One way or another some physical input is needed 
to determine the precise boundary conditions on the near horizon modes upon which the entire set 
of results and physical consequences for macroscopic black hole physics hinge. At the horizon the 
classical supposition that nothing happens at a coordinate singularity is in tension with the behavior 
and assumptions of quantum field theory (in a fixed background) at very high frequencies. This tension 



is the source of the paradoxes, since it is the classical supposition that leads to (2.18) and (2.19) which 
seem themselves to lead to the opposite conclusion that either quantum states at arbitrarily short 
distances near the horizon are playing an important physical role (unlike in flat space), or entirely 
new physics and degrees of freedom must be invoked to explain black hole entropy, undermining the 
semi-classical assumption of mild behavior at the horizon. Once arbitrarily high frequency modes near 
the horizon are admitted into the discussion, one should reconsider whether it is reasonable to treat 
gravity classically with the background geometry known and completely determined and ask whether 
the near-horizon behavior is indeed mild, or whether there might be large quantum backreaction effects 
on the local geometry of spacetime in its vicinity. Notice also that this trans-Planckian problem for 



ultrashort distance modes arises near a black hole horizon despite the fact that the horizon radius (2.4) 
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itself is quite macroscopic and very large compared to the Planck scale. 

The trans-Planckian problem and the divergence of the entanglement entropy as the black hole 
horizon is approached is also reminiscent of the ultraviolet catastrophe of the energy density of classical 



thermal radiation. The cancellation of h from dE in (2.20) is similar to the cancellation in the energy 
density of modes of the radiation field in thermodynamic equilibrium, ftw n BE {ui)uj 2 dijj — > k B Tuj 2 du 
in the Rayleigh-Jeans limit of very low frequencies, fkv <C k B T where the Bose-Einstein distribution 
n BE {uj) = [exp(hw/k B T) — — > k B T/Huj. If this low energy relation from classical Maxwell theory 
is improperly extended into the quantum high frequency regime, Hlo 3> k B T it leads to a divergent 
integral over oj and hence an infinite energy density of the radiation field at any finite temperature. 
This ultraviolet catastrophe and the low temperature thermodynamics of solids led Planck and Einstein 
to take the first steps towards a quantum theory of radiation and bulk matter [H [2]. The analogous 
high frequency divergence of the entanglement entropy near a black hole horizon suggests that it results 



from a similar improper extension and misinterpretation of the classical formula (2.15) extended into 



the high frequency, low temperature regime where quantum effects become important. 

C. Quantum Fields in Schwarzschild Spacetime 

The preceding discussion indicates that quantum effects, particularly at short distances need to be 
treated very carefully when black hole horizons are involved. Given the high stakes of the possibility of 
fundamental revision of the laws of physics and/or vast numbers of new degrees of freedom and the role 
of ultrahigh frequency trans-Planckian modes to account for the Hawking effect and black hole entropy, 
it would seem reasonable to return to first principles, and re-examine carefully the strictly classical view 
of the event horizon as a harmless kinematic singularity, when h ^ and the quantum fluctuations of 
matter are taken into account. 

Consider the basic set up of a quantum theory of a scalar field in fixed Schwarschild spacetime. 
Although the locally high energy matter self-interactions and gravitational self-interactions themselves 
are almost certainly important near the event horizon, we ignore them here in order to simplify the 
discussion. Then the generally covariant free Klein- Gordon eq., 

1 ~ / ' „ab< 



{-a+fi 2 )^ = --^=d a [V^gg ab d b ^) +^$ = 12.2s, 



for a scalar field of mass [i in the static, spherically symmetric Schwarzschild geometry (2.2) is separable, 
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with eigenfunctions of the form 



<Pwtm(t,r,i 



'2uj 



Ye m (9, 



(2.29) 



Here Yi m is a spherical harmonic, and the radial function f U £ satisfies the ordinary differential eq. 

d 2 



dr 



* 2 



fut = U fu£: 



(2.30) 



in terms of the Regge- Wheeler ( "tortoise" ) radial coordinate, 



r + r g In 



(2.31) 



with the potential, 



V e 



£{£+!) r 



+ 



(2.32) 



which may be viewed as an implicit function of r* through the relation (2.31). Note that as r ranges 



from r s to oo, r* ranges over the entire real line from — oo to +oo, and that the potential Vg vanishes at 



the lower limit, but is otherwise everywhere positive. As a corollary note from (2.32) also that at the 
horizon r 



si 



the scalar field mass fi drops out. Since we are interested in the near horizon behavior 



we may concentrate on the massless case and set [x = 0. 



Eq. (2.30) defines a standard one dimensional scattering problem, with two linearly independent 
scattering solutions that for fj, = have the asymptotic forms, e =tlwr * as r — > r s , r* — > — oo, and as 
r, r* —7- +oo. Accordingly, we may define the two fundamental linearly independent scattering solutions 



of (|2.30|) f w ^ by their asymptotic behaviors as 




V 

1 

r* 



— oo 
+oo 

— oo 
+oo 



(2.33a) 



(2.33b) 



Because of the constancy of the Wronskian associated with eq. (2.30), the reflection and transmission 



coefficients of (2.33) obey 



At*(w)Bi(u) 



\Af( U )\ 2 + \Bt( U )\^ 
\A?(u)\. 



(2.34a) 
(2.34b) 
(2.34c) 
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Because these two scattering solutions are linearly independent, independent creation and destruction 
operators must be introduced for them in the canonical quantization of the Heisenberg field operator, 



J ° Im I=L,R 



(2.35) 



The independent canonical commutation relations, 

= 2irh5(uj — aj')5u'5 mm '5 IJ 



i J] 

a ulmi a u'l'm' 



with each of I, J taking the values L, R enforce the canonical equal time commutation relation, 



(9$ 

<f>(t,r,6,cl>),—(t,r',e',<t>') 



ih 



5(r* 



r*')5(9 



r 2 sin( 



on the field, provided the normalization condition 

/oo 
dr*f% f^ = 2n5 IJ 5(u J -u J ') 
-oo 



(2.36) 



(2.37) 



(2.38) 



is satisfied. 



L R 

From (2.33) and (2.34) the non- vanishing of the reflection coefficient A f ' implies that outgoing 



spherical waves at the black hole horizon are a linear superposition of outgoing and ingoing spherical 
waves at infinity, and similarly for ingoing spherical waves. Notice that this differs from flat space in 
spherical coordinates, both in the presence of a scattering potential and the existence of two linearly 
independent regular solutions (2.33) of the radial wave equation (2.30). In flat space, r s = and 
r* = r so the corresponding wave equation has a singular point at the origin r = within the finite 



range of the radial variable. This forces one to accept only the solutions of (2.30) which are regular at 



the origin, namely rji(kr) (with k = \f uj 2 — /i 2 ) and exclude the irregular solutions whose derivatives 
diverge there. The mass (i does not drop out and there remains a gap between the positive and negative 
energy solutions in flat space. 



In contrast, in the Schwarzschild case the change of variables (2.30) and wave eq. (2.31) shows that 



the equation and both its solutions are regular at the horizon r* 



-oo. The origin r = is not 



present within the range — oo < r* < +oo at all. Hence no particular linear combination of (2.33) is 



preferred a priori, and we need to retain both solutions. The frequency integral in (2.35) also extends 



down to uj = 0. At u) = the differential eq. (2.30) admits solutions which behave linearly in r* , hence 



logarithmically as ln(r/r s — 1) near the horizon. Whereas modes behaving this way over an infinite 
domain are excluded by initial conditions with compact support, the horizon is a finite distance away 



from any point of fixed r > r s in the physical Schwarzschild line element (2.1), and hence these modes 
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are no longer excluded a priori. In several important respects, the radial wave eq. (2.30) and Hilbert 



space spanned by its solutions are discontinuously different in the flat and Schwarzschild cases. 



Since the static Schwarzschild geometry (2.2) approaches ordinary flat space as r — > oo one natural 
definition of the "vacuum" would seem to be the state annihilated by all of the a^ m , viz., 



B) 



R 



B) = 0. 



(2.39) 



This state and its Green's functions were studied in detail by Boulware [35] and is denoted here by \B). 
If one calculates the expectation value of the stress-energy tensor of the massless (conformally coupled) 
scalar field, 

1 



\ (v a $) ( v 6 $) -Jn(v$) 2 -^ v a v fe $ . 

3 o 3 



(2.40) 



in the Boulware state, one finds the usual (quartic) divergence of the vacuum energy, obtained also in 
flat space, which must be removed, and a finite remainder, which vanishes as r — > oo (as r^/r 6 ) just 
as one would expect for the Minkowski vacuum far from the black hole. However the renormalized 



expectation value of (2.40) has the property that [36 

1 R 



(B\T%\B) 



He 
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1 



diag(-3,l,l,l), 



(2.41) 



as r 



where it diverges. Thus in this Boulware state the apparent coordinate singularity of the 



Schwarzschild horizon is now the locus of arbitarily high energy densities. Clearly, if such a state were 
realized in practice, its stress-energy would act as a physical source for the semi-classical Einstein's 
equations, 

R 



R\- -5 a b + A8 a b = 8irG(T a b ) R . 



(2.42) 



(with A = here) and necessarily influence the background spacetime (2.1) which assumed T a b = 0. 



Such large stresses as present in (2.41) would cause the solution of (2.42) to deviate markedly from the 



classical Schwarzschild geometry (2.2) near the horizon, and require re-evaluation of the entire starting 



point of the discussion, and certainly the analytic continuation (2.7). 



The stress-energy in (2.41) is the negative of that of a scalar field in a thermal state at the local 



blueshifted Tolman-Hawking temperature |37j . 

T loc = T H (l-'f)" 1 ■ (2-43) 
Mathematically, the stress-energy is proportional to an integral over frequencies whose finite part is 



proportional to Tf , and the behavior (2.41) as r — > r s is obtained. This contribution to the frequency 
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integral is dominated by frequencies uj ~ c/r s , defined by (2.29) with respect to the time at infinity, 



which is the only fixed scale entering the scattering potential in (2.32) if [i = 0. (Since near the 



horizon the potential Vi vanishes and the mass fi drops out of the leading behavior of {T a b ) R as r 



SI 



all fields behave essentially as massless fields there in any case). The local frequency of these finite u 
modes becomes arbitrarily large, even exceeding the Planck scale on the horizon, which is what leads 



to the divergence in (2.41). 

It is clear that the trans-Planckian issue arises because of the infinite blue shift of frequencies at 
the event horizon, a necessary consequence of the gravitational redshift of waves followed backwards to 



their origin at the horizon, expressed in the relation (2.5). Classically, this infinite blueshift presents no 
particular problem, since the energy of classical waves can be made arbitrarily small, no matter how 
high their frequency, simply by making their amplitude small enough. As soon as h ^ (no matter 
how small), the situation is quite different, as the amplitude of quantized wave modes is bounded from 



below by the Heisenberg uncertainty relation, encoded in the commutation rules ( 2.36 )-( 2.37). The 
local energy of the wave mode with local frequency 0Ji OC ai is 

\ he 



Elocal = hlOlocal = ?Ud [l 



1 



(2.44) 



which also diverges on the horizon. Since energy-momentum couples universally to gravity, the very 
large local vacuum zero point energy can affect the geometry there. Let us emphasize that this large 



effect derives from the choice of state in (2.39), and cannot be removed by a coordinate transformation, 
once the state has been specified. In the Boulware state the finite vacuum polarization effects and their 
backreaction on the geometry near the horizon are very large in any coordinates. 

0, Ei oca i 



The relation (2.44) shows that the limits H — > and r 



r s do not commute. If h 



0. 



for all r > r s , and one might entertain the logical possibility of analytically continuing the exte- 
rior Schwarzschild geometry into the interior region, by extending the notion of general coordinate 
invariance for real, differentiable coordinate transformations, x^ — > x /fl (x) to complex meromorphic 
transformations, and get around coordinate singularities on the real axis. However, the behavior of the 
quantum vacuum zero point energy near the horizon depends on arbitrarily high local frequencies and 



is not smooth. In the Boulware state \B) it diverges, (2.41). Hence analytic continuation around the 
coordinate singularity there may not be physically justified in the quantum theory, and certainly not if 
the matter field is in this state. 

Hawking and Unruh argued for a different state in the gravitational collapse problem, different from 
the Boulware state and one which is no longer time symmetric |20|. 138] . In that state, denoted by \U), 
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only the L ingoing modes are taken to be in the vacuum in the first half of (2.39), but the R outgoing 
modes in the final state are thermally populated at infinity with T = T . The additional finite thermal 



flux has a stress-energy tensor that just cancels the diverging negative energy density (2.41) of the 
Boulware state near the future horizon, T = +X (in any proper set of regular coordinates there), and 
gives a positive flux of Hawking radiation to infinity. Indeed, the Unruh state | U) is constructed by the 
requirement that its "vacuum" modes are analytic in the Kruskal null coordinate U = T — X across 
the future horizon. The Hawking flux in the Unruh state may be thought of as bringing the quantum 
expectation value up to its vacuum value at the future horizon where the Unruh state is locally similar 
to the flat space vacuum. This adjustment necessarily produces a non- vacuum state of flux of real 
quanta at infinity. 

In this way, the Hawking-Unruh state maintains the regularity (and smallness) of the stress-energy 
tensor on the future horizon, consistent with the assumption of negligible local backreaction of the 
radiation on the spacetime geometry itself. To obtain this result, however, one must use the same set of 



modes (2.33) in either case, and follow them to arbitrarily large local frequencies near the horizon, with 
a specific boundary condition of analyticity in U = T— X there. This produces an outward Hawking flux 
by a compensating negative energy flux through the future horizon and into the coordinate singularity 
at r = 0. In other words, one must assume that ordinary quantum field theory and the wave equation 



(2.28) holds on a fixed classical background geometry with arbitrary accuracy all the way down distances 
of the order than and even arbitrarily smaller than the Planck scale Lpi, at which short distances one 
would normally expect the semi-classical theory to break down completely. 

Whereas for example in the Coulomb scattering problem one assumes regular wave functions, uni- 
tarity and no flux into or out of the Coulomb singularity (and a formally divergent stress-energy density 
there), in the Schwarzschild case the Hawking-Unruh state assumes a regular stress-tensor on the future 
horizon which necessarily requires a negative energy flux through it and down to the singularity, vio- 
lating unitarity. Conversely, the Boulware state has no energy flux through the horizon but a diverging 



stress-energy there (2.41). Finally the Hartle-Hawking-Israel state is a thermal one for both the L and 
R modes, and therefore has no flux into or out of the singularity, and a regular stress tensor on the 
horizon, but it is a state thermally populated with both ingoing and outgoing quanta even at infinity 
[39H41] . As we already have seen, such a state is thermodynamically unstable. Unlike in flat space, 
there is no choice of boundary conditions which satisfies all three criteria of i) regularity on the future 
horizon, ii) zero flux there (and hence zero flux into the future singularity at r = 0), and iii) vacuum-like 
at infinity. So an inescapable conclusion is that at least one of these three criteria must be abandoned, 



25 



but pure mathematics cannot tell us which. 

Despite the apparently thermal expectation values of the two states \B) and \U), each are pure 
states related to each other by a Bogoliubov transformation. In the case of \U) the apparently thermal 
emission is consistent with a pure state because precise quantum mechanical phase correlations are set 
up and maintained between the modes outgoing to infinity and those infalling into the future singularity. 
The pure state becomes a mixed thermal state if and only if one sums over the modes localized behind 
the horizon as unobservable [42-44]. Of course, such an averaging procedure entails giving up any hope 
of keeping track of the correlations that might exist between the radiated quanta at different times. 
It is also somewhat paradoxical that although information seems to be "lost" by the pair creation 
process in which one member of the pair falls into the black hole, the mass of the hole and hence the 



Bekenstein-Hawking entropy (2.24) is decreased by the thermal emission. This is very different from a 
normal thermal emission process from a star such as the sun for example. In thermal equilibrium the 
star's radiant energy is supplied by its nuclear interactions in its core, and simply passed outward at 
a steady rate. Neither the temperature profile nor the total entropy of the sun changes in this steady 
state process, and the change is its mass is completely negligible. 

Later authors have shown that the stress-tensor for a thermally populated state at an arbitrary 
temperature T ^ T H behaves like 

{T%) * ^hj^^-^K'-'f) diag ( - 3 ' h h 1} ' (2 - 45) 

near the horizon. This divergence and its disappearance if the temperature equals the Hawking temper- 



ature can be understood geometrically. From the Kruskal coordinate transformation (2.7), we observe 
that if the Schwarzschild static time coordinate t is continued to imaginary values, then the resulting 
Euclidean signature Riemannian manifold has a conical singularity at r = r s unless the Euclidean time 
variable is periodic with periodicity an integer multiple of 47rr s /c = (3 H . Likewise in order to avoid 
a singularity in the Green's functions and stress-energy tensor of a quantum field on this background, 
they too must have this Euclidean periodicity. This is nothing other than the Kubo-Martin-Schwinger 
(KMS) Euclidean periodicity condition |46j for the thermal state of a quantum field theory at temper- 
ature h/k B f5 H =T H . In fact this is one way in which the Hawking temperature was intuitively arrived 
at in ref. |39^ 142]. If T ^ T H , the Euclidean periodicity of the thermal state does not match 2Trh/k B T H , 



and the conical singularity at r = r s leads to the divergence in (2.45). 

Although the divergence in the renormalized expectation value is cancelled if the modes are populated 



with a thermal distribution at a temperature precisely equal to the Hawking temperature, c.f. (2.45), 
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our discussion of fluctuations in the previous section leads us to expect that even a slight variation of 
the temperature away from the mean value of (T a b ) will produce very large fluctuations in the energy 
density near the horizon. Fluctuations are intrinsic to both thermal and quantum theory, and require 
calculating {T a b {x)T c d (x')) , for a linear response treatment of their effect on the spacetime for their 
quantitative analysis [47j. One would expect that the natural time scale for the instability at the 



horizon to develop in such an analysis is the only dynamical timescale available, given by (2.23) and 
therefore very rapid. 

The minimal conclusion from these considerations is that the macroscopic quantum physics of black 
holes is quite delicately dependent upon what one assumes for the population of trans-Planckian fre- 
quency modes in the near vicinity of the black hole horizon. Depending upon how they are treated, by 
boundary conditions, either these ultrahigh frequencies are responsible for the thermal evaporation, or 
they cause the stress-tensor to diverge and can produce significant backreaction. The results depend 
radically on the choice of state, and the correct physics can be determined only if fluctuations about 
typical states are studied in a systematic way in the collapse problem. To date this investigation has 
not been carried out. In any case it is clear that one cannot maintain vacuum boundary conditions 
at both the horizon and large distances from the black hole. Thus the obstruction to a global vacuum 
in Schwarzschild spacetime has a topological character, related to the possible appearance of a conical 
singularity on the horizon and involving singular coordinate transformations there. 

Examples of singular coordinate transformations, conical- like singularities and global obstructions 
associated with genuine physical effects are known in other areas of quantum physics. In classical 
electromagnetism the gauge potential is unmeasurable locally and may be gauged away. Only the 
field strength tensor F^ u = d^A u — d u A^ and functions of it are gauge invariant and locally measurable. 
However, the circulation integral § c A^dx^ is also gauge invariant, and quantum mechanically therefore 
so is the phase factor 

W(C) = exp (jt j> A^dxA (2.46) 

around any closed loop C. It measures the magnetic flux threading the loop. If W(C) ^ 1 , then the 
gauge potential cannot be set equal to zero everywhere along C by a regular gauge transformation, 
even if the local electromagnetic field evaluated at all points along C vanishes. The singular gauge 
transformation which is necessary to set A^ = physically corresponds to the creation or destruction 
of a magnetic vortex in a superconductor, which would be pointlike and have infinite energy (were 
it not for its normal, non-superconducting core of finite radius). The requirement that the complex 
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valued electron pair density be single valued around any closed loop leads to flux quantization of 

magnetic flux in superconductors [48], with q = 2e for the Cooper pair [4"9l [50]. 



A second example of the physical relevance of the non-local phase factor (2.46) is the Aharonov 



Bohm effect [52] , which shows that the interference pattern of electron waves passing around a solenoidal 
magnetic field confined to a certain region of space is affected by the presence of the field in the interior 
enclosed by the interfering trajectories, even if the electrons' classical trajectories are confined to the 
region where the local field strength tensor F^ v vanishes identically. The non-local gauge invariant phase 



factor (2.46) has physical consequences for interference of electron waves that do not depend upon the 



strength of the local field along the classical electron trajectory. 

Both of these physical effects of gauge potentials which can be gauged away locally but not globally 
are expressed mathematically by the statement that QED is properly defined as a theory of a U(l) fiber 
bundle over spacetime. Depending upon the topology of the base manifold of spacetime (for example 
whether or not it is "punctured" by excising a region where the magnetic field of the Aharonov-Bohm 
solenoid is non-zero), the topology of this fiber bundle may be non-trivial and non-local gauge invariant 



quantities such (2.46) can carry information about physical processes. 



The global quantum effect of blueshift near a black hole horizon has a topological aspect which is 
similar. Although the contractions of the local Riemann curvature tensor remain finite in the classical 



Schwarzschild geometry (2.6), this static geometry has a timelike Killing field with components, 

K a = (1,0,0,0) 



in the static coordinates (t,r,0,4>) of (2.2). The norm of this Killing vector is 



(2.47) 



(2.48) 



exactly the gravitational redshift (blueshift) factor appearing in (2.5) or (2.44), and to the inverse 



fourth power in (2.41). The quantum state of the system, specified in Fourier space by (2.39) retains 



this global information about the infinite blueshift at the horizon relative to the standard of time in the 



asymptotically flat region, r — > oo, because of the existence of the global timelike Killing field (2.47). 



This generator of time translation symmetry has been used in defining the Boulware state (2.39) to 



distinguish positive from negative frequencies, and hence distinguish particle-waves from antiparticle- 



waves in the quantum theory. The norm (2.48) is a completely coordinate invariant (albeit non-local) 



scalar quantity, not directly related to the local curvature. Hence there is no reason of coordinate 
invariance that precludes it from having physical effects, and in particular, large physical effects at the 
horizon in a state like \B). 
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The horizon where the norm (2.48) vanishes has topological significance. On the Euclidean section 



t — > it with it — > it + f3 H periodically identified as suggested by Hawking and the Unruh state boundary 
conditions, the Euler characteristic, defined in terms of the Riemann dual tensor *R a bcd = ^abefR e c J^, 
is 



Xi 
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1 
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— r 2 dr 



(2.49) 



where we have used /3 H = h/k B T H = 4-7rr s (with c = 1), (2.6) and the vacuum Einstein eqs. R a b = R = 



for the Schwarzschild solution. The Euclidean Schwarzschild manifold with period j3 H and r > r s has 



the topology IR 2 x § 2 , unlike flat M 4 . This is reflected in the Euler characteristic (2.49), and the doubling 



of regular solutions to the radial eq. ( 2.30 ) . General theorems in differential geometry relate the number 



of fixed points of a Killing field where the norm (2.48 ) vanishes to the Euler characteristic of the manifold 
|51j . so that x E = 2 is associated with the vanishing of (2.48) at r = r s . A periodicity condition on 



the orbits of the Killing field (2.47), particularly in the complexified domain, in order to eliminate the 



conical singularity at r s when the period of it is different than f}„ is completely non-trivial at the 
quantum level and an a priori unwarranted assumption. It is analogous to assuming the triviality of the 



U(l) bundle and of the phase factor (2.46), which would miss genuine physical effects such as Abrikosov 
vortices and the quantization of circulation and magnetic flux in superfluids and superconductors |50| . 
and the Bohm-Aharanov effect in QED |52j. 

In the gravitational case the possible non-triviality of the GL{A) bundle of General Relativity is 



encoded in the the fact that the Euler density in (2.49) can be expressed as the total divergence of a 



frame dependent topological current [53], dual to an anti-symmetric 3- form gauge potential, 



* r> * ryabcd 
i^abcd -K 



V a Q a = V a (e abcd A bcd ) , 



(2.50) 



and the surface integral of this gauge potential over a closed bounding 3-surface, 



£ abcd A a bc dT,d 



is coordinate (gauge) invariant, under the gauge transformation, 



(2.51) 



A-abc A.abc + V[ a 6> bc ] 



(2.52) 
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Thus if £ is the S 2 ® M. tube at fixed r from ii to t2, the integral (2.51) measures the topological 



Gauss-Bonnet- Chern charge residing in the interior of the tube, much as the circulation integral in the 



exponent of (2.46) measures the magnetic flux from a superconducting vortex or an Aharonov-Bohn 



solenoid threading its interior. 



Concerning the Riemann tensor itself, we note that in the general static metric of (2.1), the tensor 
component 



^-S(S-?-£f 



(where primes denote differentiation with respect to r) becomes —f"/2 and hence remains finite when 
the two functions / = h are equal. However, Einstein's eqs. in the static geometry, 

1 ri 

-G\ = ——[r(l-h)} = -8irG T\ = 8ttG P , (2.54a) 
r z ar 

G\ = — e %- + \ (h - 1) = 8vrG T. = 8nGp , (2.54b) 
rj ar r A 

give 

i(j) = -^ G{p+p) T (2 - 55) 

Hence if the non- vacuum stress-energy has p+p > in the region where / or h vanishes, in general h ^ f 
and the cancellation of the singularities at h = f = 0, special for static vacuum solutions to Einstein's 
eqs., will not occur. Any perturbation of the vacuum Schwarzschild spacetime with p + p ^ in a 
static frame in the vicinity of the horizon has the potential to produce Riemann tensor perturbations, 
5R tr 

tr ~ ( r ^f) i which are genericaUy large at the horizon, where / — > 0, and thus will produce a large 
change in the local geometry there. Further, the eq. of stress-energy conservation in a static, spherically 
symmetric spacetime, 

V a n = ^ + + = 0, (2.56) 

ar 2} ar r 

(where p± = T e d = is the transverse pressure while p = T r T is the radial pressure) . This shows that 
any matter with the effective eq. of state p = p± = wp must have a stress-energy which behaves like 



/-(l+«)/2w j which diver 

ges on the horizon if w > 0. Note that this is consistent with p.4l| for w — 
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producing a conical singularity. 

From these various considerations we conclude that the cancellation of divergences on a black hole 
event horizon is an extremely delicate matter, and there is no reason to expect it to be generic in 
quantum theory. This does not violate the Principle of Equivalence, if that principle is understood to 
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be embodied in the requirement that physics should respect general coordinate invariance under all real, 
non-singular coordinate transformations. Singular coordinate transformations are another matter. 

Mathematically, gauging away coordinate singularities on the horizon amounts to an additional as- 
sumption about the triviality of the GL(A) frame bundle of General Relativity, which may not be war- 
ranted. Experimentally well-established applications of quantum theory already teach us by a number 
of examples in other fields that such improper "gauge" transformations generally contain new physical 
effects. Physically these effects are associated with the quantum wavelike nature of matter which can- 
not be idealized as arbitrarily small, isolated pointlike particles, particularly in gravity where such an 
extreme local limit must produce infinite energies. Instead, matter fields obey wave equations such as 



(2.28) which require boundary conditions for their full solution. In quantum theory the dependence of 



local physics on these boundary conditions, through specification of the quantum state of the system 



as in (2.39) does not violate the Principle of Equivalence. It is only our uncritical and unexamined 
classical notions of strict and absolute locality that are violated in known quantum phenomena such 
as entanglement and macroscopic coherence. Expectation values of the stress-energy tensor T a b can 



perfectly well depend on non-local invariants such as (2.48) and (2.51), analogous to (2.46) in QED. A 
new set of calculation tools is needed in order to determine these quantum effects in a systematic way, 
and bring gravity into accordance with general quantum principles on macroscopic scales. This is what 
we seek to provide in the succeeding sections. 

III. THE CHALLENGE OF COSMOLOGICAL DARK ENERGY 

The second challenge for macroscopic quantum effects in gravity arise on the cosmological scale of 
the Hubble expansion itself, and in particular upon the discovery of the acceleration of the expansion 
rate of the universe. 

A. The Cosmological Constant and Energy of the Vacuum 

In classical General Relativity, the requirement that the field equations involve no more than two 
derivatives of the metric tensor allows for the possible addition of a constant term, the cosmological 
term A, to Einstein's equations, 

R\-^s\ + A n = ^n- (3.i) 
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If transposed to the right side of this relation, the A term corresponds to a constant energy density p\ = 
c 4 A/8irG and isotropic pressure p\ = —c 4 A/8ttG permeating all of space uniformly, and independently 
of any localized matter sources. Hence, even if the matter T a b = 0, a cosmological term causes spacetime 
to become curved with a radius of curvature of order |A|~2. 

With A = there is no fixed length scale in the vacuum Einstein equations, G/c 4 being simply a 
conversion factor between the units of energy and those of length. Hence in purely classical physics 
there is no natural fundamental length scale to which A can be compared, and A may take on any value 
whatsoever with no difficulty (and with no explanation) in classical General Relativity. 

As soon as we allow h 7^ 0, there is a quantity with the dimensions of length that can be formed 



from h,G, and c, namely the Planck length (2.17). Hence in quantum theory the quantity 

, 9 hGA 

A = AL 2 pl = -3- (3.2) 

becomes a dimensionless pure number, whose value one might expect a theory of gravity incorporating 
quantum effects to address. Notice that like the effects we have been considering in black hole physics 
this quantity requires both h and G to be different than zero. 

Some eighty years ago W. Pauli was apparently the first to consider the question of the effects of 
quantum vacuum fluctuations on the the curvature of space [33]. Pauli recognized that the sum of zero 
point energies of the two transverse electromagnetic field modes in vacuo, 

»» (f 3 k huik 1 he 

contribute to the stress-energy tensor of Einstein's theory as would an effective cosmological term A > 0. 



"*= 2 / w?— = ^ijr-" A (3 ' 3) 



Since the integral (3.3) is quartically divergent, an ultraviolet cutoff £~* n of (3.3) at large k is needed. 
Taking this short distance cutoff L m i n to be of the order of the classical electron radius e 2 /mc 2 , Pauli 
concluded that if his estimate were correct, Einstein's theory with this large a A would lead to a universe 
so curved that its total size "could not even reach to the moon" |54| . If instead of the classical electron 
radius, the apparently natural but much shorter length scale of L m i n ~ L p i is used to cut off the 



frequency sum in (3.3), then the estimate for the cosmological term in Einstein's equations becomes 
vastly larger, and the entire universe would be limited in size to the microscopic scale of L p i itself, in 
even more striking disagreement with observation. 

Clearly Pauli's estimate of the contribution of short distance modes of the electromagnetic field to 



the curvature of space, by using (3.3) as a source for Einstein's eqs. (3.1) is wrong. The question is 
why. Here the Casimir effect may have something to teach us [55 . The vacuum zero point fluctuations 
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being considered in (3.3) are the same ones that contribute to the Casimir effect, but this estimate 
of the scale of vacuum zero point energy, quartically dependent on a short distance cutoff L m i n , is 
certainly not relevant for the effect observed in the laboratory. In calculations of the Casimir force 
between conductors, one subtracts the zero point energy of the electromagnetic field in an infinitely 
extended vacuum (with the conductors absent) from the modified zero point energies in the presence 
of the conductors. It is this subtracted zero point energy of the electromagnetic vacuum, depending 
upon the boundary conditions imposed by the conducting surfaces, which leads to experimentally well 
verified results for the force between the conductors. In this renormalization procedure the ultraviolet 
cutoff £~* n drops out, and the distance scale of quantum fluctuations that determine the magnitude of 
the Casimir effect is not the microscopic classical electron radius, as in Pauli's original estimate, nor 
much less the even more microscopic Planck length L p i, but rather the relatively macroscopic distance 
d between the conducting boundary surfaces. The resulting subtracted energy density of the vacuum 
between the conductors is [55] 

7T 2 he 



This energy density is of the opposite sign as (3.3), leading to an attractive force per unit area between 
the plates of 0.013 dyne/cm 2 (^to/gT) 4 , a value which is both independent of the ultraviolet cutoff £~* n , 
and the microscopic details of the atomic constituents of the conductors. This is a clear indication, 
confirmed by experiment, that the measurable effects associated with vacuum fluctuations are infrared 
phenomena, dependent upon macroscopic boundary conditions, which have little or nothing to do with 



the extreme ultraviolet modes in (3.3). 

Actually, the original Casimir calculation of the force between exactly parallel flat plates of infinite 
conductivity hides some important features of the general case. As soon as the conducting plates have 
any finite curvature, the local stress-energy tensor diverges on the boundary. A general classification of 
these divergences has been given in Ref. [56J. This shows that in the presence of curved boundaries, 
there is residual sensitivity to ultraviolet effects, much as in the divergence of the stress-energy tensor 



(2.41) in the Boulware state in the Schwarzschild geometry. By now the meaning of these divergences 
and the correct physical method of handling them have been understood |57j . The divergences are 
artifacts of the theoretical over-idealization of the conductors as perfect at arbitrarily high frequencies. 
Any real metal has a finite conductivity which leads to a finite skin depth and vanishing reflection 
coefficient at arbitrarily high frequencies. When the effects of finite conductivity of real metals such as 
gold used in the experiments are taken into account, all local stresses and energy densities are finite, and 
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the theoretical predictions in accord with experiment |57j . Thus although the Casimir effect itself and 



in its original form (3.4) is a macroscopic effect of vacuum zero point energy, the details do retain some 
sensitivity to the physics of the surface, which is short distance compared to the macroscopic separation 
d (although of course completely unrelated to an ultrashort distance cutoff on the scale of the electron 
radius or Lpi). 

By the Equivalence Principle, local ultrashort distance behavior in a mildly curved spacetime is 
essentially equivalent to that in flat spacetime. Hence on physical grounds we should not expect the 



extreme ultraviolet cutoff dependence of (3.3) to affect the universe in the large any more than it affects 
the force between metallic conductors in the laboratory, although any possible surface boundary effects 
will have to be treated carefully. 

In the case of the Casimir effect a constant zero point energy of the vacuum, no matter how large, 
does not affect the force between the plates. In the case of cosmology it is usually taken for granted 
that any effects of boundary conditions can be neglected. It is not obvious then what should play the 
role of the conducting plates in determining the magnitude of p v in the universe, and the magnitude of 
any effect of quantum zero point energy on the curvature of space has remained unclear from Pauli's 
original estimate down to the present. In recent years this has evolved from a question of fundamental 
importance in theoretical physics to a central one of observational cosmology as well. Observations 
of type la supernovae at moderately large redshifts (z ~ 0.5 to 1) have led to the conclusion that the 
Hubble expansion of the universe is accelerating [58J. According to Einstein's equations this acceleration 
is possible if and only if the energy density and pressure of the dominant component of the universe 
satisfy the inequality, 

p + 3p = p (l + 3u>) < 0. (3.5) 

A vacuum energy with p > and w = p v /p v = —1 leads to an accelerated expansion, a kind of 
"repulsive" gravity in which the relativistic effects of a negative pressure can overcome a positive energy 



density in (3.5 ). Taken at face value, the observations imply that some 74% of the energy in the universe 
is of this hitherto undetected w = —1 dark variety. This leads to a non-zero inferred cosmological term 
in Einstein's equations of 

3ff 2 r 3 
A meas ~ (0.74) -l~1.4x 10- 56 cm" 2 ~ 3.6 x lO -122 — . (3.6) 
c z nG 

Here Hq is the present value of the Hubble parameter, approximately 73 km/sec/Mpc ~ 2.4x 10 -18 sec -1 . 



The last number in (3.6) expresses the value of the cosmological dark energy inferred from the SN la 
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data in terms of Planck units, L p 2 = ^j, i.e. the dimensionless number in (3.2) has the value 



A ~ 3.6 x 1CT 122 . (3.7) 

Explaining the value of this smallest number in all of physics is the basic form of the cosmological 
constant problem. 

B. Classical de Sitter Spacetime 

Just as in our discussion of black hole physics in Sec. [IlJ we begin our discussion of cosmological 
dark energy with a review of the simplest classical spacetime that serves as the stage for discussion of 
quantum effects of a cosmological term. The maximally symmetric solution to the classical Einstein 



equations (3.1) with a positive cosmological term is the one found by W. de Sitter in 1917 [59J. In fact, 



de Sitter expressed his solution of the vacuum Einstein eqs. (3.1) with A > and T a b = in a static, 
spherically symmetric form analogous to that of the Schwarzschild solution (2.1). In those coordinates 
de Sitter's solution can be presented as 

.,2 



f( r ) = h{r) = 1 - H 2 r 2 = 1 - ^ , (3.8) 



instead of (2.2). The length scale r H is related to the cosmological term by 



In this spherically symmetric static form, the de Sitter metric has rotational 0(3) symmetry about the 
point r = and an apparent horizon at r = r H . 

Actually the de Sitter solution has a larger 0(4, 1) symmetry which can be made manifest in its 
global analytic extension, analagous to the T, X Kruskal-Szekeres coordinates in the Schwarzschild 
case. Consider a five dimensional Minkowski space with the standard flat metric, 

ds 2 = rj AB dX A dX B = -dT 2 + dW 2 + dX 2 + dY 2 + dZ 2 , (3.10) 

subject to the condition, 

r] AB X A X B = —T 2 + W 2 + X 2 + Y 2 + Z 2 = r 2 H = H~ 2 . (3.11) 

with the indices A,B = 0, 1,2, 3, 4 raised and lowered with the five dimensional Minkowski metric rj AB = 
diag (— 1,1,1,1,1). Here and henceforward we generally set the speed of light c = 1, except when needed 
for emphasis. 
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FIG. 2: The de Sitter manifold represented as a single sheeted hyperboloid of revolution about the T axis, in 
which the X 1 , X 2 coordinates are suppressed. The hypersurfaces at constant T are three-spheres, § 3 . The 
three-spheres at T = ±00 are denoted by I± . 

If T were a spacelike coordinate rather than timelike so that the four dimensional manifold defined 
by (3.10) and (3.11) had a Euclidean signature, we would clearly be discussing a four-sphere S with 



an 0(5) symmetry group. Because of the Lorentzian signature, these relations instead define a single 
sheeted hyperbolid of revolution, depicted in Fig. [2} with the non-compact symmetry group 0(4, 1), 



the maximal possible symmetry for any solution of the vacuum Einstein field equations (3.1) with a 
positive cosmological term and T a b = in four dimensions. 
The curvature tensor of de Sitter space satisfies 



R 



ab 



cd 

J?" 
K b 
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3H 2 5 a 
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(3.12a) 
(3.12b) 
(3.12c) 



The Lie algebra so(4, 1) is generated by the generators of the Lorentz group in the 4 + 1 dimensional flat 



embedding spacetime (3.10). In the coordinate (adjoint) representation the 10 anti-Hermitian generators 
of this symmetry are 

d „ d 



L AB -X A —-X B — 



-L r 



These 10 generators satisfy 



\Lab-, Lcd] — ~Vac L bd + rj BC L AD — rj BD L AC + T] AD L BC , 



(3.13) 



(3.14) 



the Lie algebra so(4, 1). De Sitter space has 10 Killing vectors corresponding to these 10 generators. 
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The hyperbolic coordinates of de Sitter space are defined by 



T = — sinhu, (3.15a) 
H 

1 

W = — cosh u cos x, (3.15b) 
H 

1 

X 1 = -— cosh u sin x fi l , i = l,2, 3. (3.15c) 
H 



where 



n = (sin 9 cos (ft , sin 6 sin </> , cos ( 



is the unit vector on S 2 , and cast the de Sitter line element in the form 

ds 2 = — 2 [-du 2 + cosh 2 u (dx 2 + sin 2 x d£l 2 )] . 
The quantity in round parentheses is 

dQl = [d(sin x n 1 )} 2 + [d(cos x)] 2 = d X 2 + sin 2 X dtt 2 



(3.16) 



(3.17) 



(3.18) 



the standard round metric on S . Hence in the geodesically complete coordinates of (3.15) the de 



Sitter line element (3.17) is an hyperboloid of revolution whose constant u sections are three spheres, 



represented in Fig. [2j which are invariant under the 0(4) subgroup of 0(4, 1). It is sometimes convenient 



to define the hyperbolic conformal time coordinate v by secu = coshn, so that (3.17) becomes 

~ 2 - H- 2 sec 2 v (-dv 2 + dn 2 3 ) , 



ds 2 



(3.19) 



conformal to the Einstein static cylinder with —tt/2 < v < ir/2. 

In cosmology it is more common to use instead the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
line element with flat M 3 spatial sections, viz. 



ds 2 



-dr 2 + a 2 (r) dx • dx = -dr 2 + a 2 (r) (dx 2 + dy 2 + dz 2 ) 
-dr 2 + a 2 (T){dg 2 + g 2 dft 2 ) . 



De Sitter space can be brought in the FLRW form by setting 

If 1\ Ha 
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(3.20) 



(3.21a) 

(3.21b) 
(3.21c) 
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with 



o(t) 
Q 



\/x 2 + 



y 2 + z A 



(3.22a) 
(3.22b) 



From (3.21) and (3.22), T + W > in these coordinates. Hence the flat FLRW coordinates cover only 



one half of the full de Sitter hyperboloid, with the hypersurfaces of constant RW time r slicing the 
hyperboloid in Fig. [2] at a 45° angle. 

The change of time variable to the conformal time coordinate, 



V 



-H~ l e~ H ' ! 



1 



1 



is also often used to express the de Sitter line element is the conformally flat form 

1 



ds 2 = a 2 {-dr, 2 + dx 2 ) 



H 2 V 2 



/From (3.15) and (3.21), 



cosh u sin \ = Hq a 
sinh u + cosh u cos x = a 



(-dr] 2 + dx 2 ). 



Q 
V 

1 

"Hrj 



which gives the direct relation between hyperbolic coordinates and flat FLRW coordinates. 
The de Sitter static coordinates (t, r, 6, <fi) are defined by 

T = \z\J\- H 2 r 2 smh(Ht) , 
H 



W 

X 1 



H 



r n 



(3.23) 



(3.24) 



(3.25a) 
(3.25b) 



(3.26a) 

(3.26b) 
(3.26c) 



They bring the line element (3.10) into the static, spherically symmetric form (2.1) with (3.8), i.e. 

dr 2 



-(1 - H 2 r 2 )dt 2 + 



(1 - H 2 r 2 ) 



+ r 2 dn 2 . 



(3.27) 



Just as in the Schwarzschild case (2.2) these static coordinates cover only part of the analytically fully 



extended de Sitter manifold (3.11). From (3.26), real static (t, r) coordinates cover only the quarter of 
the de Sitter manifold where W > and both W ± T > 0. This quarter is represented as the rightmost 
wedge of the Carter-Penrose conformal diagram of de Sitter space in Fig. [3j 
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FIG. 3: The Carter- Penrose conformal diagram for de Sitter space. Future and past infinity are at I±. Only the 



quarter of the diagram labeled as the static region are covered by the static coordinates of (3.27). The orbits of 



the static time Killing field d/dt are shown. The angular coordinates 9, <j> are again suppressed. 

The Regge- Wheeler radial coordinate r* can be denned in the static frame by 

l + Hr s 



1 
2H 
1 

H 



In 



1 



tanh -1 (Hr) , 



dr* 



1-Hr J H 
tanh(-ffr*) so that 

\/l — H 2 r 2 = sech(Hr*) , 



dr 



1 - H 2 r 2 ' 



(3.28a) 
(3.28b) 
(3.28c) 



and 



ds 2 



sech 2 ( Hr*) i-dt z + dr* 2 ) + — , tanh^Hr*) dQ, 2 

H A 

sech 2 (Hr*) \-H 2 dt 2 + H 2 dr* 2 + sinh 2 (Hr*) dfl 2 } 
H 2 1 J 



(3.29) 



Note that the horizon at r = if 1 = r H is mapped to r* = oo in these coordinates, and that the spatial 



part of the line element (3.29) in square brackets (the "optical" metric) is 



ds 2 opt = H 2 dr* 2 + smh 2 (Hr*) dtf = 4 ^ ^ 



(3.30) 



where the second form is obtained by defining 



= yn. 



tanh 



Hr* 



Hr 



Hl + y 2 ' 



H 



In 



H 2 r 2 

i + y 
l -v 



so that 



(3.31a) 
(3.31b) 

(3.31c) 
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Eq. (3.30) is a standard form of the line element of three-dimensional Lobachewsky (hyperbolic or 



Euclidean anti-deSitter) space H 3 . Thus, one expects conformal field theory (CFT) behavior at the 



horizon boundary, y = 1, r = r H . This is reflected also in the wave eq. analogous to (2.30) in which the 



corresponding potential Vi of (3.44) vanishes at the horizon, which we consider in the next subsection. 

Since these various coordinatizations of the de Sitter manifold are quite different globally, and involve 
singular coordinate transformations at the horizon, much of the discussion of the Schwarzschild geometry 
have their analogs in de Sitter spacetime. The Carter-Penrose conformal diagram for the analytically 
extended de Sitter hyperboloid, Fig. [3] is similar to the corresponding diagram in the Schwarzschild 
case, Fig. [TJ In each case the horizon is bifurcate (i.e. has two distinct parts) and the region covered 



by the static coordinates (3.27) is duplicated by a second region centered on the antipodal point of S 



where the sense of increasing static time t is reversed. Thus the static Killing vector K a (2.47) becomes 
null on either horizon and spacelike in the upper and lower quarter wedges of Fig. [3j 

C. Quantum Effects in de Sitter Spacetime 

Quantum fluctuations and their backreaction effects in de Sitter spacetime were considered in |47l I60T - 
163] . These studies indicate that fluctuations at the horizon scale r H are responsible for important 
backreaction effects on the classical de Sitter expansion that could potentially relax the effective cos- 
mological vacuum energy to zero. The authors of Refs. [6l] have performed a perturbative analysis of 
long wavelength gravitational fluctuations in non de Sitter invariant initial states up to two-loop order. 
This work indicates the presence of secular terms in the quantum stress tensor of fluctuations about 
de Sitter space, tending to decrease the effective vacuum energy density, consistent with the earlier 
considerations in Refs. |47 [ I60ff5"3"] . The authors of Ref. [S3] have studied the stress tensor for long 
wavelength cosmological perturbations in inflationary models as well, and also found a backreaction 
effect of the right sign to slow inflation. See also Refs. [66], [67], and a discussion of these various 
approaches in Ref. [68 1. 

One of the first indications of non-trivial quantum effects in de Sitter spacetime is particle creation 
|60j . A space or time dependent background field generally creates particles. Schwinger first studied 
this effect in QED in a series of classic papers [5T]. The rate of spontaneous decay of the electric field 
into charged particles is 

(eE) 2 ( m 2 c 3 \ 

e*P "T^F (3-32) 



ch 2 ir 2 V eEh 



in the limit eE -C m 2 c 3 /h. From this point of view the exponential de Sitter expansion ( 3.20 )-( 3.22) 
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provides a time dependent background field which can create particle pairs from the "vacuum," con- 
verting the energy of the classical gravitational background into that of particle modes. The rate of this 
spontaneous creation of matter in de Sitter space can be calculated in analogy to the Schwinger effect 
in an electric field, with a similar result for the decay rate per unit volume |60j . 



16tf 4 / 27rm\ 
r=-^exp^- w j (3.33) 

for a scalar massive field with arbitrary curvature coupling (i.e. m 2 = p? + £R = p 2 + 12^H 2 ). 
In the case of a spatially uniform electric field (and no magnetic fields) the Maxwell equation, 

M=-® < 3 ' 34 > 

indicates that the creation of these particle pairs tends (at least initially) to decrease the strength of 
the electric field. In cosmology the Friedmann equation, 

d\ 2 tt2 8ttG 



, H 2 = —p, (3.35) 
together with the equation of covariant energy conservation, 

p + 3H(p + p) = 0, (3.36) 
(with over dots here denoting differentiation with respect to proper time r) imply that 

H = - 4 ^(p + p). (3.37) 
In both cases there is a classical static background that solves the equation trivially, namely H or E a 



constant with zero source terms on the right hand side of (3.34) or (3.37). In the case of (3.37) this is 
de Sitter spacetime with pa + pa = 0. Any creation of matter with p + p > will tend to decrease the 
strength of the classical background field H. The particle creation process has been studied in greater 
detail in the electric field [69] and de Sitter [70J cases, and backreaction effects were also taken into 
account in the QED case the large N limit [69j . In this limit self-interactions between the created charge 
particles are neglected, but such scatterings are essential to the final decay of the coherent classical field 
into particles. In the gravitational case it has so far not been feasible to carry out a full dynamical 
calculation of backreaction and particle scattering which would establish p+p > and final decay of the 
background value of H . Unlike electromagnetism in gravity there are massless particles which couple to 



the background field, and from (3.33) these will give the largest effect. However for minimally coupled 



massless scalar particles and for gravitons themselves, the definition of the "vacuum" state becomes 
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more subtle, so that computing the infrared effects of massless particles and their interactions requires 
a different approach. 

One possible way to define a preferred state of a system is by its symmetries. Since the symmetry 
group of the maximally analytically extended de Sitter hyperboloid is 0(4,1), a maximally 0(4, 1) 
symmetric state can be defined, for most field theories, including massive and massless conformally 
invariant fields, such as the photon. Since under analytic continuation of T — > iT the de Sitter manifold 
becomes § 4 , which is compact, one can define a state by the requirement of maximal 0(5) invariance, 
regularity on S 4 , whose re-point functions are all analytic under the continuation. This is often referred 
to as the Bunch-Davies (BD) state in the literature [71]. Since the Euclidean S 4 radius is r H , the 
Euclidean Green's functions are periodic in imaginary time with period 2irr H /c. This is again the KMS 
property of a thermal Green's function, and hence the BD state defined by this analytic continuation 
from S 4 is a thermal state with temperature |421 14"4"] , 

r « = ^' < 3 - 38 > 

the Hawking temperature of de Sitter space. Hence the BD state in de Sitter space is not a "vacuum" 
state, but is a state thermally populated with quanta at temperature T„. In this respect it is rather 
like the Hartle-Hawking state in the Schwarzschild background. In these fully time symmetric states 
there is no net particle creation on average and no decay rate. Because of the full 0(4, 1) symmetry 
of the BD state, the expectation value of the energy-momentum tensor of any matter fields in this 
state must itself be of the form of a quantum correction to the cosmological term with p = —p ~ H 4 
constant. The BD state usually assumed, tacitly or explicitly in models of inflation. It is interesting 
that because of infrared divergences the BD state does not exist for massless, minimally coupled scalar 
fields and gravitons themselves [72] , which implies that these fields must be in a state with less than the 
full 0(4, 1) de Sitter symmetry. Once the de Sitter symmetry is broken, however, it is not clear what 
the residual symmetry and final state of the system is, or how to go about determining it. 

Even the existence of a maximally symmetric state does not guarantee its stability against small 
fluctuations. For example, it can be shown that in the uniform constant electric field background there 
are exactly 10 isometries, the same number as with zero external field [73]. Thus one can construct 
mathematically a state which respects all these isometries, including a discrete symmetry corresponding 
to time reversal invariance, and no net pair creation, without necessarily guaranteeing the stability of 
the vacuum against pair creation in an external electric field. A subtle point, easily overlooked in 
Schwinger's elegant effective action method, is that either time asymmetric boundary /initial conditions 
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must be specified in such a decay problem, or it must be shown that the time symmetric state is 
dynamically unstable to quantum fluctuations, spontaneously breaking the larger isometry group to a 
smaller subgroup. 

One may construct an argument analogous to that in the Schwarzschild case that the de Sitter 
invariant BD should be quantum mechanically unstable, under fluctuations in its Hawking temperature. 
Note that the vacuum energy within a spherical volume of radius r H is 

4irr 3 c 4 1 



which is inversely proportional to (3.38). Therefore 



dE„ E„ Trc 3 k„ 

and the heat capacity of the region of de Sitter space within one horizon volume is apparently negative. 
This indicates that if the region can exchange energy with its environment external to the cosmological 
horizon, the BD thermal state will be unstable to such energy exchanges, analogously to the black hole 
case. The problem of negative heat capacity is also similar. 

It has also been suggested that the Bekenstein-Hawking 1/4 area formula for the entropy 

S — k Ah - ^ K (3 411 

be taken over to the de Sitter case [74j . However what degrees of freedom this "entropy" counts are 



even less clear than in the black hole case. Note also that from (3.39) 



dE H = d(p A V H ) = dp a V H + p\ dV H 
= dp A V H - pAdV H 

+ T H dS H - PA dV H . (3.42) 
In fact T H dS H and dp A V H have opposite signs. This implies that S H cannot be interpreted as an entropy 



and/or additional contributions, such as surface terms, are missing from (3.42). 

Like the Schwarzschild case, de Sitter spacetime admits a static Killing field which is timelike in 
one region of its maximal analytic extension. If one considers any thermal state with a temperature T 
different from T H , defined with respect to the static Killing time in one patch of de Sitter space, then 
the renormalized stress-energy tensor is 



(n) * ^ 9^ (^3 " T " ] (1 " FV ^ 2 diag ( " 3 ' h h 1} ' (3 - 43) 
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as r — > r H , analogous to (2.45) in the Schwarzschild case. Thus any finite deviation of the temperature 



from the Hawking-de Sitter temperature (3.38) in the region interior to the horizon will produce a 



arbitrarily large stress-energy on the horizon. The previous discussion in Sec. [XT] about the dependence 
of {T a b ) R on the gauge invariant but non-local norm of the static Killing vector field, the non-commutivity 
of the limits H — )■ and r — > r H , and the breakdown of the analytic continuation hypothesis through 
coordinate or conical singularities apply to the de Sitter case as well. 

The wave equation of a quantum field propagating in de Sitter space can also be separated in the 



static coordinates (3.27). By following the analogous steps as used in (2.28)-(2.32) in the Schwarzschild 



case, one arrives at the identical form of the radial scattering eq. (2.30) in terms of the r* coordinate 



(3.28), with the scattering potential in de Sitter space given by 

£ + 1) 



dS 



H 



+ // + 2^(6£-l) 



1 



£{£ + 1) csch^ (Hr*) + ( - v 2 ] sech 2 (Hr*) 



where — £/2 is the i?<£ 2 coupling in the scalar field Lagrangian, and 



v 



I'- 
ll-' 



12f . 



(3.44) 



(3.45) 



A significant difference of the effective one dimensional scattering problem of spherical waves in the de 
Sitter case is that the coordinate singularity at the origin r = is in the physical range, compared to 
the Schwarzschild case in which the r* coordinate range (— oo, oo) covers only the exterior Schwarzschild 
geometry. Thus in order to avoid a singularity at the origin we must require that the scattering solutions 



to (2.30) with (3.44) satisfy 



Tu)i\ 



as 



as 



(3.46) 



thereby excluding a possible singular r * behavior for £ > 1. This means that in de Sitter space only the 



particular linear combination of ingoing and outgoing solutions of the corresponding (2.30) vanishing 



at the origin should be used in the quantization of the field and not the general linear combination 



in ( 2.33 )-( 2.35) appropriate in the Schwarzschild background. As a corollary this implies that there can 



be no analog of the Unruh state in de Sitter spacetime, since (3.46) is equivalent to the requirement 



of no net flux into or out of the origin at r = in static de Sitter coordinates (3.27). Otherwise the 



quantization of a scalar field in de Sitter space in these coordinates is analogous to the Schwarzschild 
case, and one can again easily find states regular at the origin which have diverging stress tensors on 



the horizon, as in (3.43). The appearance of such states and large stress tensors would necessarily mean 
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large quantum backreaction effects in the vicinity of the horizon r = r H , and the breakdown of 0(4, 1) 
de Sitter symmetry down to 0(3), with perhaps a very different global geometry than the analytic 



continuation implied by extending the coordinates (3.17) globally. As in the Schwarzschild case, this is 
ultimately a question of physics, not simply mathematical analytic continuation of coordinates. 

The earlier work on the behavior of the graviton propagator in de Sitter space |75H77] also indicates 
the existence of infrared divergent contributions to correlations at large distances. This violates cluster 
decomposition and makes the definition of a graviton scattering matrix in global de Sitter space prob- 
lematic. The existence of long distance correlations and infrared divergences is a signal of the breakdown 
of the global BD state. The relevance of matter self-interactions, first studied in Refs. \78\ 179] has been 
taken up again recently in Refs. [80] . where the phenomenon of runaway stimulated emission in de 
Sitter space is explored. This adds to the by now substantial literature on non-trivial quantum infrared 
behavior in de Sitter space [68]. Since infrared effects are the common feature of all these studies, one 
might well suppose that there should be some relevant operator(s) in the low energy effective theory 
of gravity which describes in a general and universal way (i.e. independently of specific matter self- 
interactions) the quantum effects of matter /radiation fields on the macroscopic scales of both the black 
hole and cosmological horizons, and which points to a possible resolution of the physics in both cases. 
This is what a predictive effective field theory (EFT) approach should provide. 



IV. EFFECTIVE FIELD THEORY: THE ROLE OF ANOMALIES 

The important role of quantum anomalies in the low energy effective field theory approach to the 
strong interactions is reviewed in this section, and their relevance to the EFT of gravity discussed. The 
first and best known example of this is the axial anomaly in QED and QCD. Since the two-particle 
correlations that give rise to non-trivial infrared effects can be seen already in this flat space case, we 
digress and review that here before returning to the curved spacetime application in gravity. 



A. The Axial Anomaly in QED and QCD 

Although the triangle axial anomaly in both QED has been known for some time |81H83] , the general 
behavior of the amplitude off mass shell, its spectral representation, the appearance of a massless 
pseudoscalar pole in the limit of zero fermion mass, and the infrared aspects of the anomaly generally 
have received only limited attention [84, 85J. It is this generally less emphasized infrared character of the 
axial anomaly which will be important for our EFT considerations in gravity, so we begin by reviewing 



45 



this aspect of the axial anomaly in QED in some detail. 
The vector and axial currents in QED are defined by 1 

JP(x) = i^{x)Y^{x) . (4.1a) 

Jg(x) = ^(x)^ 5 rp(x) . (4.1b) 

The Dirac eq., 

-i^{d IM -ieA^)i; + mtp = 0. (4.2) 
implies that the vector current is conserved, 

V = 0, (4.3) 

while the axial current apparently obeys 

J§ = 2im ijjj ijj (classically). (4-4) 

In the limit of vanishing fermion mass m — > 0, the classical Lagrangian has a U c h(l) global symmetry 
under tp — > e iai5 ip, in addition to £7(1) local gauge invariance, and J§ is the Noether current correspond- 
ing to this chiral symmetry. As is well known, both symmetries cannot be maintained simultaneously 
at the quantum level. Let us denote by (J§(z)) A the expectation value of the chiral current in the 



presence of a background electromagnetic potential A^. Enforcing U(l) gauge invariance (4.3) on the 
full quantum theory leads necessarily to a finite axial current anomaly, 



e 2 _ _ e 2 



d ^ ,„o = IS? ^" F " F " = 2^ E - B - <"' 5 » 

in an external electromagnetic field. Varying this expression twice with respect to the external A field 
we see that anomaly must appear in the amplitude, 

V^( v a ) = -i f d A x ( d 4 v e ip ' x+iq ' y ^ J ^°^ A 
^ q) - J J V 5A a (x)5A p {y) 



A=0 



ie 2 



J d 4 x J d 4 ye^ + ^ (T4(0)J a (x)J?(y))\ A=0 . (4.6) 



At the lowest one-loop order it is given by the triangle diagram of Fig. [4j plus the Bose symmetrized 
diagram with the photon legs interchanged. 



1 We use the conventions that {7^,7"} = ~ 2 g 1 *" = 2diag (H ), so that 7 = (7°)*, and -y 5 = i^j 1 ^ 3 = (7 5 ) 1 " are 

hermitian, and tr(7 5 7 A1 7 I/ 7 P 7 <T ) = — 4ie A1 " p<T , where e^ uprT — — e M „ p(T is the fully anti-symmetric Levi-Civita tensor, with 
£0123 = +1- 
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FIG. 4: The triangle diagram contributing to the axial current anomaly. The four-momentum of integration is I. 
Elementary power counting indicates that the triangle diagram of Fig. [4] is superficially linearly 



divergent. The formal reason why (4.3) and (4.4) cannot both be maintained at the quantum level 



is that verifying them requires the ability to shift the loop momentum integration variable I in the 
triangle amplitude. Because the diagram is linearly divergent, such shifts are inherently ambiguous, 
and can generate finite extra terms. It turns out that there is no choice for removing the ambiguity 
which satisfies both the vector and chiral Ward identities simultaneously, and one is forced to choose 
between them. Thus although the ambiguity results in a well-defined finite term, the axial anomaly has 
most often been presented as inherently a problem of regularization of an apparently ultraviolet linearly 
divergent loop integral [8TH83] . 

There is an alternative derivation of the axial anomaly that emphasizes instead its infrared character. 
The idea of this approach is to use the tensor structure of the triangle amplitude to extract its well- 
defined ultraviolet finite parts, which are homogeneous of degree three in the external momenta p and 
q. Then the remaining parts of the full amplitude may be determined by the joint requirements of 
Lorentz covariance, Bose symmetry under interchange of the two photon legs, and electromagnetic 
current conservation, 



P. 



a r^(p,q) = = q^(p,q), (4.7) 



at the two vector vertices. By this method the full one-loop triangle contribution to F IJ,a ^(p, q), becomes 
completely determined in terms of well-defined ultraviolet finite integrals which require no further 
regularization [86, 87J. The divergence of the axial current may then be computed unambiguously, 



and one obtains (4.5) in the limit of vanishing fermion mass. It is this latter method which makes 



clear that the anomaly is a consequence of symmetries of the low energy theory, no matter how its 
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UV behavior is tamed, provided only that the regularization respects these symmetries. There is of 
course no contradiction between these two points of view, since it is the same Ward identities which are 
imposed in either method, and in any case in the conformal limit of vanishing fermion mass the infrared 
and ultraviolet behavior of the triangle amplitude are one and the same. 

The details of the calculation by this method may be found in Refs. [86-88J. One first uses the 
Poincare invariance of the vacuum to assert that Y^{p,q) can be expanded in the set of all three- 
index tensors constructible from the p and q, with the correct Lorentz transformation properties. There 
are exactly eight such tensors, two of which are linear in p or q, namely e^ a "^p\ and e^ a ^q\, while the 
remaining six are homogeneous of degree three in the external momenta. However only certain linear 



combinations of these eight tensors satisfy (4.7). Define first the two index tensor, 

v a ^{p,q)=e a ^ a p p q ai (4. 

which satisfies 

v a(3 (p,q) = v^{q,p) } 
p a v^(p,q)= =q^(p,q). 

Then the six third rank tensors, if* (jp, q), i = 1, . . . , 6 which satisfy the conditions (4.7), 

PaT? aP (p, q) = = if °P(p, q)qp = 0, i = l,...,6, 

are given in Table [I| 



(4.9a) 
(4.9b) 

(4.10) 



i 




1 


-p ■ q e» a P x px - p 13 v^ a (p, q) 


2 


pV Q/3 V +p a v^(p,q) 


3 




4 


p . qe ^pX qx+q a v nPfa g) 


5 




6 


ql* v a ^{p,q) 



TABLE I: The 6 third rank (pseudo)tensors obeying (4.10) 



Hence we may express the amplitude (4.6) satisfying (4.7) as a linear combination, 

6 



(4.11) 



1=1 
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where fi = fi(k 2 ;p 2 ,q 2 ) are dimension —2 scalar functions of the three invariants, p 2 ,q 2 , and k 2 . Note 
from Table[I]that the two tensors of dimension one which could potentially have logarithmically divergent 
scalar coefficient functions occur only in linear combination with dimension three tensors. Hence the 



coefficient functions fi of these linear combinations obeying (4.7) are all ultraviolet finite. These finite 
contributions can be obtained unambiguously from the imaginary part of the triangle graph of Fig. [5j 
which are finite a priori and then determining the real part from the imaginary part by an unsubtracted, 
i.e. UV finite dispersion relation. This leads to the finite coefficients given in the literature |86H88j . viz. 



dy^ (4.12a) 
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x(l — x) 
D 

y(i - y) 

D 



(4.12b) 

(4.12c) 
(4.12d) 



where the denominator of the Feynman parameter integral is 

D = p 2 x(l — x) + q 2 y(l — y) + 2p ■ q xy + m 2 = (p 2 x + q 2 y)(l — x — y) + xy k 2 + m 2 , (4.13) 

strictly positive for m 2 > 0, and spacelike momenta, k 2 ,p 2 ,q 2 > 0. 
Thus, the full amplitude Y^ a "{p,q) satisfying 

(i) Lorentz invariance of the vacuum , 

(ii) Bose symmetry under interchange of photon lines , 



(hi) vector current conservation (4.7) 



(iv) unsubtracted dispersion relation of real and imaginary parts , 
with the finite imaginary parts determined by the cut triangle diagram of Fig. [5j is given by (4.11 ) and 



(4.12), without any need of regularization of ultraviolet divergent loop integrals at any step. 

Since this amplitude is fully determined by (i)-(iv), its contraction with k^, and the divergence of 
the axial vector current is determined as well. There is no further freedom to demand the naive axial 



Ward identity corresponding to (4.4). Instead we find 

k ll T^(p,q)=Av a P(p,q), (4.14) 
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FIG. 5: The discontinuous or imaginary part of the triangle diagram of Fig. |4]with respect to k 2 . The propagators 
which are cut by the dashed line are replaced by their corresponding on-shell delta functions. 

where 



A(k 2 ;p 2 ,q 2 ) = 2p-qf 1+ p 2 f 2 +q 2 f 5 



1 , f l ~ x , D-m 2 

. dx dy — 

V 2 Jo Jo D 



e 

2^2 



7T- 



m 



dx 



l-x 



dy 



D 



(4.15) 



The second term proportional to m 2 is what would be expected from the naive axial vector divergence 



(4.4) jSTj. The first term in (4.15) in which the denominator D is cancelled in the numerator is 

„2 t-l r-l-x „2 



tt2 , 
7T .l 



dx 



dy 



o 



e 

2^' 



(4.16) 



and which remains finite and non-zero in the limit m — > is the axial anomaly. 

Thus the finite anomalous term is unambiguously determined by our four requirements above, and 
may be clearly identified even for finite m, when the chiral symmetry is broken. This construction 
of the amplitude from only symmetry principles and its finite parts may be regarded as a proof that 
the same finite axial anomaly must arise in any regularization of the original triangle amplitude which 
respects these symmetries and leaves the finite parts unchanged. Explicit calculations in dimensional 
regularization and Pauli-Villars regularization schemes, which respect these symmetries confirm this 

Eg. 



Now the important consequence of the anomaly in (4.15) in the divergence is that when the photons 



q 2 = are on shell, the amplitude T fia ^(p, q) develops a simple pole at in k 2 as the electron mass 



m 



0. This follows from the first line of (4.15) which shows that when p 2 = q 2 = 0, 2p ■ q = k 2 so A 
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is explicitly proportional to k 2 . Since A is finite when m — > 0, the amplitude function fi(k 2 ; 0, 0) must 
develop a pole in k 2 . Indeed by explicit calculation from (4.12)-(4.13), we have 



h(k 2 ;p 2 = 0,q 2 = 0)\ m=0 = £- 2 l 2 . (4.17) 

The corresponding imaginary part (discontinuity in k 2 ) becomes a S(k 2 ) in the same limit. Moreover, 
even for arbitrary p 2 ,q 2 ,m 2 > 0, the spectral function obtained from this imaginary part obeys an 
ultraviolet finite sum rule [571 EH] • 



The appearance of a massless pseudoscalar pole (4.17) in the triangle anomaly amplitude in the 
massless fermion limit suggests that this can be described as the propagator of a pseudoscalar field 
which couples to the axial current. Indeed it is not difficult to find the field description of the pole. To 
do so let us note first that the axial current expectation value ( Jg ) A can be obtained from an extended 
action principle in which we introduce an axial vector field source, B^ into the Dirac Lagrangian, 

itp^ (jd^ -ieA^ Tp-mtp^ (d^ -ieA^ - igj 5 B^j ip - mtpip (4.18) 

so that the variation of the corresponding action with respect to B^ gives 

' vS g(4) A • (4.i9) 



Henceforth we shall set the axial vector coupling g = 1. Next let us decompose the axial vector B^ into 
its transverse and longitudinal parts, 

B„ = Bi + dpB (4.20) 
with d^Bj^ = and B a pseudoscalar. Then, by an integration by parts in the action corresponding to 



(4.18), we have 



d,(4) A = ~%- (4-21) 



Thus the axial anomaly (4.5) implies that there is a term in the one-loop effective action in a background 



A^ and B^ field, linear in B of the form, 



e 2 



Seff = ~ 16^2 / * x ^ paF ^ F P* B > ( 4 - 22 ) 



or since d x B\ = DB, 

Seff = j d 4 x j d 4 y , [e^F, v F pa ] x n~t [d x B x ] y , (4.23) 
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where is the Green's function for the massless scalar wave operator □ = d^d^. Thus from (4.19), 
this non-local action gives [90] 

= 4^ d " n ~ leaPpc7FapFp ° ' (424) 



which exhibts the massless scalar pole in the massless limit of (4.17), and which agrees with the explicit 
calculation of the physical (0\J§\p,q) triangle amplitude to two photons for p 2 = q 2 = m 2 = 0. The 
existence of this pole or 5 function in the two-particle intermediate e + e~ intermediate state is the 
first indication that a phenomenon similar to Cooper pairing in condensed matter systems can exist in 
relativistic quantum field theory, and this effect is connected to anomalies. 



In Ref. |88| a local action corresponding to (4.61) was found by introducing two pseudoscalar fields. 



The introduction of two fields is necessary if the only term in the effective action with this massless 



pole is the term (4.23) which is off-diagonal in FF and d ■ B. If on the other hand the effective action 
contains the perfect square, 

_! -...a a 



d 4 x [d x B x + n^F^F^j 

then this effective action can be represented by a single pseudoscalar field. The possible existence of the 
additional terms in massless QED necessary to complete the square is currently under investigation, 

Since it contains kinetic terms for the additional pseudoscalar field(s), the effective action of Ref. 
[88] describes additional pseudoscalar degree(s) of freedom. These degrees of freedom are two-particle 
CP correlated e + e~ states, composite bilinears of ip and ip, which appear in anomalous amplitudes as 
massless poles. In condensed matter physics, or electrodynamics at finite temperature or in polarizable 
media, where Lorentz invariance is broken, it is a familiar circumstance that there are low energy col- 
lective modes of the many body theory, which are not part of the single particle constituent spectrum. 
This occurs also in vacuo in the two dimensional massless Schwinger model, whose anomaly and longi- 
tudinal "photon" can be described by the introduction of an effective scalar field composed of an e + e~ 
pair [91]. In 3 + 1 dimensions, relativistic kinematics and symmetries severely limit the possibilities for 
the appearance of such composite massless scalars, with the triangle anomaly the only known exam- 
ple [S2]. The fact that the e + e~ pair becomes collinear in the massless limit shows that this effectively 
reduces the dimensionality back to 1 + 1. In the well studied 1 + 1 dimensional case, the commutation 
relations of fermion bilinear currents and Jg , which create the composite e + e~ massless state are 
due to the anomaly [93] . A similar phenomenon occurs in the triangle amplitude in 3 + 1 dimensions. 
The axial anomaly thus implies additional long range correlations and collective degrees of freedom in 
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the many body quantum theory, similar to Cooper pairs in a superconductor, which are not present in 
the classical or first quantized single particle theory. 

In real QED these infrared effects are suppressed by the non-zero physical electron mass m > 0, and 
the additional fact that macroscopic chirality violating sources for J§ which would be sensitive to the 
anomaly are difficult to create. In QCD the situation is complicated both by the strong interactions in 
the infrared and the spontaneous breaking of chiral symmetry. The neutral member of the isotriplet of 
pseudoscalar Goldstone bosons in the low energy EFT is the 7r°, whose decay to two photons, tt° — > 27 
is correctly given by the triangle amplitude |94| 195] . In fact, it was the experimental agreement between 
the measured decay rate to that predicted by the axial anomaly computed in the UV theory of 3 colors 
of fractionally charged quarks that gave one of the strongest early confirmations of QCD. The fact 
that this amplitude is non- vanishing in the chiral limit, yet cannot be described by a local operator of 
dimension 4 or less in the chiral Lagrangian, violates naive decoupling and illustrates how an anomaly 
couples UV to low energy physics. It is the fact that the anomaly may be computed in the UV theory 
of QCD but gives rise to a low energy amplitude of meson decay, tt° — > 27 that led to the principle of 
anomaly matching [96 . 



The apparent massless pseudoscalar anomaly pole of (4.17) in the isosinglet channel in the chiral limit 
of QCD is even more interesting. The 0~ state described by this pole in the isosinglet channel mixes 
with the pseudoscalar axial gluon density Q(x) = G'l lu (x)G a ^ u (x), and gives rise to a non-vanishing 
susceptibility of axial gluon densities, 

X (k 2 ) = J d 4 xe lk - x (Q(x)Q(0)), (4.25) 

as k 2 — > 0, despite the fact that Q is a total derivative and therefore one would naively expect x(k 2 ) 
to be proportional to k 2 and vanish in this limit. The fact that the susceptibility is non-vanishing 
is a direct effect of the massless anomaly pole |97j . The degree of freedom this infrared pole represents 
combines with a non-dynamical but gauge invariant Q 2 term in the effective action of QCD to yield 
finally one propagating massive isosinglet psuedoscalar state which can be identified with the rj 1 meson, 
solving the U(l) problem in QCD [97] • Thus there is no doubt that the pseudoscalar 0~ state which 
appears in the isosinglet anomaly channel in perturbation theory is physical and propagating in the 
final S-matrix of the theory, but it becomes massive by a topological variant of the Higgs mechanism 

The lesson to be taken away from this QCD example is that anomalies are a unique window which 
the low energy EFT provides to short distance physics. Two particle correlated pair states appear in 
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anomalous amplitudes, which can be described as propagating massless fields that have consequences 
for low energy physics. The anomalous Ward identities and the long distance effects they generate must 
be taken into account by explicitly adding the IR relevant terms they induce in the low energy effective 
action [MIES]. 



B. The (TJJ) Triangle Amplitude in QED 

In this section we consider the amplitude for the trace anomaly in flat space that most closely 
corresponds to the triangle amplitude for the axial current anomaly reviewed in the previous section, 
and give a complete calculation of the full (T^ 11 J a jP) amplitude for all values of the mass and the 
off-shell kinematic invariants. Although the tensor structure of this amplitude is more involved than 
the axial vector case, the kinematics is essentially the same, and the appearance of the massless pole 
very much analogous to the axial case. 

Classical fields satisfying wave equations with zero mass, which are invariant under conformal trans- 
formations of the spacetime metric, g^ u — > e 2a g a b have stress tensors with zero classical trace, = 0. 
In quantum theory the stress tensor Tt becomes an operator with fluctuations about its mean value. 



The mean value itself (Tt) is formally UV divergent, due to its zero point fluctuations, as in (3.3), 
and requires a careful renormalization procedure. The result of this renormalization consistent with 
covariant conservation in curved spacetime is that classical conformal invariance cannot be maintained 
at the quantum level. The trace of the stress tensor is generally non-zero when h ^ 0, in non-trivial 
background fields, provided one preserves the covariant conservation of T% (a necessary requirement 
of any theory respecting general coordinate invariance and consistent with the Equivalence Principle) 
yields an expectation value of the quantum stress tensor with a non-zero trace. 

The fundamental quantity of interest for us now is the expectation value of the energy-momentum 
tensor bilinear in the fermion fields in an external electromagnetic potential A^, 

(Tn A = (T%J A + (TZ) A (4.26) 

where 

T free = -^H^H + SH^T^A^ " "#V>) (4.27a) 

Tj% = - eJ^A^ + eg» u J x A x , (4.27b) 



are the contributions to the stress tensor of the free and interaction terms of the Dirac Lagrangian (4.18). 
The notations, t^ v > = it^ v + t u ^)/2 and di_i = (<9^ — <9^)/2, for symmetrization and ant-symmetrization 
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have been used. The expectation value (T'*") A satisfies the partial conservation equation, 

d v (Tn A = eF^{J u ) A , 



(4.28) 



upon formal use of the Dirac eq. of motion (4.2). Just as in the chiral case, the relation is formal 



because of the a priori ill-defined nature of the bilinear product of Dirac field operators at the same 



spacetime point in (4.27). Energy-momentum conservation in full QED (i.e. when the electromagnetic 



field is also quantized) requires adding to the fermionic T^ v of (4.27) the electromagnetic Maxwell 
stress tensor, 



1 Max A ^9 r r. 



X P 



(4.29) 



which satisfies d u T^ ax = —F^ U J U . This cancels (4.28) at the operator level, so that the full stress 



tensor of QED is conserved upon using Maxwell's eqs., d u F^ = J^. Since in our present treatment 
is an arbitrary external potential, rather than a dynamical field, we consider only the fermionic parts 



of the stress tensor (4.27) whose expectation value satisfies (4.28) instead 



At the classical level, i.e. again formally, upon use of (4.2), the trace of the fermionic stress tensor 
obeys 



-mipijj (classically) , 



(4.30) 



analogous to the classical relation for the axial current (4.4). From this it would appear that (T^ 



will become traceless in the massless limit m — > 0, corresponding to the global dilation symmetry of the 
classical theory with zero mass. However, as in the case of the classical chiral symmetry, this symmetry 
under global scale transformations cannot be maintained at the quantum level, without violating the 



conservation law satisfied by a related current, in this case the partial conservation law (4.28), implied 



by general coordinate invariance. Requiring that (4.28) is preserved at the quantum level necessarily 



leads to a well-defined anomaly in the trace [99-101J, 



III A \ m =Q 



24vr2 



(4.31) 



analogous to (4.5). It is the infrared consequences of this modified, anomalous trace identity and the 



appearance of massless scalar degrees of freedom for vanishing electron mass m = 0, analogous to those 
found in the axial case that we will study. 



The one-loop triangle amplitude analogous to (4.6) of the axial anomaly must satisfy vector current 
conservation, 



p a T^(p, q )=qpT^(p, q )=0. 



(4.32) 
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and the (partial) conservation law of the fermion stress-tensor (4.28) gives the Ward identity, 

KT^ip, q) = isTp,, - O^)n^(g) + (g"V - 8^W v {p) , (4.33) 
or since the polarizarion U^ u (p) is also a correlator of conserved currents, 

U a(3 (p) = (p 2 g al3 - p a /)n(p 2 ) . (4.34) 

we obtain 

K T^ al3 (p, q) = (gVV - q»g al3 p 2 + g^q a p 2 - g^p a p ■ g) U(p 2 ) 

+ (p fl q a qP - p^g^q 2 + ^ W - g^q^p ■ q\ U(q 2 ) . (4.35) 

These relations are still formal since the one-loop expressions are formally divergent. However, analo- 
gously to the axial case, the joint requirements of: 

(i) Lorentz invariance of the vacuum, 

(ii) Bose symmetry, Y^{p,q) = T^ ul3a (q,p) 



(hi) vector current conservation (4.32), 



(iv) unsubtracted dispersion relation of real and imaginary parts, and 



(v) energy- momentum tensor conservation (4.35), 



are sufficient to determine the full amplitude r^" Q ^(p, q) in terms of its explicitly finite pieces, and yield 
a well-defined finite trace anomaly. As in the axial anomaly case conisdered previously, this method of 
constructing the full T fiual3 (p,q) may be regarded as a proof that the same finite trace anomaly must 
be obtained in any regularization scheme that respects (i)-(v) above. It is particularly important to 
recognize that the last condition (v) is necessary to obtain a covariantly conserved stress tensor and 
avoid any gravitational anomaly (i.e. breaking of general coordinate invariance) from arising at the 
quantum level. If this fifth condition is not applied, and the naive conformal Ward identity arising from 



(4.30) is used instead, there is a gravitational anomaly, general coordinate invariance is broken, but 
scale invariance is preserved and the /3 function of the electromagnetic or QCD coupling would vanish 
|88|, I99|. in contradiction with experiment [102J. 

The tensor analysis in this case is somewhat more involved than in the axial case and is given detail 
in Ref. [88]. Lorentz invariance of the vacuum (i) is again assumed first by expanding the amplitude 
in terms of all the possible tensors with four indices depending on p a , q@ and the flat spacetime metric 
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gap _ _ j) e £ ne th e two- index tensors, 

u^(p,q) = {p-q)g aP -q a ^, (4.36a) 

w a(3 (p, q) = p 2 q 2 g al3 + (p • g)j*V - - p 2 ?V , (4.36b) 

each of which satisfies the conditions of Bose symmetry, 

u^(p,q)=u^(q,p), (4.37a) 

w afS (p,q)=w^(q,p), (4.37b) 

and vector current conservation, 

P*u a P(p, q)=0 = qpu a P{p, q) , (4.38a) 

p a w a P(p, q)=0 = qpw a P(p, q) . (4.38b) 

Making use of u a ^(p,q) and w a ^(p,q), one finds that there are exactly 13 linearly independent four- 
tensors tj Ua ^(p, q), i = 1, . . . , 13, which satisfy 

PaC^fo l)=0 = Q^ip, q), i = 1, . . . , 13 . (4.39) 

These 13 tensors are catalogued in Table [II} 

Only the first two of the thirteen tensors possess a non-zero trace, 

V C /3 (p,g) = 3fc 2 ^(p,< ? ) I (4.40a) 

9^ vaf3 (p,q) = Zk 2 w a P(p, q ), (4.40b) 

while the remaining eleven tensors are traceless, 

9^C a0 M = O, i = 3,...,13. (4.41) 

In the limit of zero fermion mass, the entire trace anomaly will reside only in the first amplitude function, 
F^e-p^q 2 ). 

To proceed, one expands Y^ ua "{p, q) in terms of these 13 tesnors with scalar coefficient functions Fi 



of the invariants k 2 ,p 2 ,q 2 analogous to (4.11), and then fixes as many of the 13 scalar functions Fi as 
possible by examining the finite terms in the formal expressions for the triangle amplitude. This turns 
out to give enough information to fix 10 linear combinations of the 13 scalar functions. The information 
needed to fix the remaining three functions comes from our fifth and final requirement on the amplitude 



namely the Ward identity (4.35). In this way, after some algebra one finds that all the Fi are completely 
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i 




1 


(k 2 g^ -k^k v )u af3 {p.q) 


2 


(k 2 g^ u - k^k u ) w al3 (p.q) 


3 


(p 2 g^ -Ap^p v ) u af) (p.q) 


4 


(p 2 g^ v -Ap^p v ) w af3 {p.q) 


5 


(q 2 g^ -Aq^q v )u af3 {p.q) 


6 


(q 2 g^ -4q>*q») w afi {p.q) 


7 


[p ■ qg^ - 2(ql J -p v + p»q v )] u af3 {p.q) 


8 


[p ■ qg>" v - 2{q^p v + p»q»)} w al3 (p.q) 


9 


(p-qp a - p 2 q a ) [p P {q»p v + p^q") -p-q (gP u pf + g'^p")} 


10 


(p-qq - q V) [q a {q^v v + p'V ) -p-q + g at "q u )] 


11 


(p ■ qp a - p 2 q a ) [2 - q 2 {g fiu q l * + g^ l q v )} 


12 


(p-qq p - q 2 p p ) [2p a p^p u - p 2 {g au p^ + g a ^p v )} 


13 


(p»q u +p v q^)g af} +p- q {g av g & ^ + g a »g0 v ) - g» v u afi 
-(/V + g Pt *P u )q a - {g ai/ q^ + g atl q v )p p 



TABLE II: The 13 fourth rank tensors satisfying (4.39) 



determined and hence there is no remaining freedom in the trace part F\ . The result of this calculation 
is [88] 

.2 r i rl-x (D-m 2 ) 



3k 2 F x 



e 

6^2 



2 



dx 
o Jo 



e 2 m 2 



dy (1 - Axy) 

-X-x 



2vr 2 J dX J - D 
The second term which vanishes in the limit m — > is the non-anomalous part which one would have ex- 



dy 



D 

(1 - Axy) 



(4.42) 



pected by naive application of the tree level Ward identity (4.30). The first term where the denominator 
D appears also in the numerator of the integrand and cancels, giving a well-defined contribution to the 
trace independent of m is the anomaly. As in the axial case it again implies the existence of a pole in k 2 
for F% in the full amplitude in the massless limit (and when the photons are on-shell: p 2 = q 2 = 0). In 
the imaginary part of the (TJJ) triangle amplitude a 5 function develops in the corresponding spectral 
density when p 2 = q 2 = and m — > + . Otherwise there is again an exact ultraviolet finite sum rule for 
this spectral density |88j . showing that the state persists even away from the conformal on-shell limit. 

The kinematics of the state appearing in the imaginary part and spectral function in this limit is 
essentially 1 + 1 dimensional, and can be represented as the two-particle collinear e + e~ pair in Fig. 
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[6j This is the only configuration possible for one particle with four-momentum k^ converting to two 
particles of zero mass, p 2 = q 2 = as k 2 — > as well. Although this special collinear kinematics 
is a set of vanishing measure in the two particle phase space, the 6(s) in the spectral function and 
finiteness of the anomaly itself shows that this pair state couples to on shell photons on the one hand, 
and gravitational metric perturbations on the other hand, with finite amplitude. 




FIG. 6: The two particle intermediate state of a collinear e + e pair responsible for the 1/fc 2 pole and <5-fh. in 
the k 2 discontinuity of the (TJJ) triangle amplitude. 

The two-particle correlated state of e + e~ behaves like an effective massless scalar exchange, and can 
be so described in an effective field theory approach. Because this effective field theory contains massless 
fields, it can have long distance, macroscopic effects. Before giving the general form of this effective 
theory in four dimensional curved space, let us consider the somewhat simpler case of two dimensions. 



C. The Trace Anomaly and Quantum Gravity in Two Dimensions 

In two dimensional curved space the trace anomaly takes the simple form, [103] . 

( Ta a) = 7^R, (d = 2) (4.43) 

where A = N$ + Np is the total number of massless fields, either scalar (Ag) or (complex) fermionic 
(Np). The fact that the anomalous trace is independent of the quantum state of the matter field(s), and 
dependent only on the geometry through the local Ricci scalar R suggests that it should be regarded as 
a geometric effect. However, no local coordinate invariant action exists whose metric variation leads to 



(4.43). This is important because it shows immediately that understanding the anomalous contributions 
to the stress tensor will bring in some non-local physics or boundary conditions on the quantum state 
at large distance scales. 



A non-local action corresponding to (4.43) can be found by introducing the conformal parameteri- 
zation of the metric, 

9ab = e 2a g ab , (4.44) 
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and noticing that the scalar curvature densities of the two metrics g a b and g a b are related by 



R- 



-g = R* 



-gUa . 



(d = 2) 



(4.45) 



a linear relation in a in two (and only two) dimensions. Multiplying (4.43) by y/—g, using (4.45) and 



noting that ^—g{T a a ) defines the conformal variation, 5T^ /5a of an effective action T^ 2 \ we conclude 
that the a dependence of can be at most quadratic in a. Hence the Wess-Zumino effective action 
[104| in two dimensions, T\J 7 is 



,(2) 
WZ 



N I' — 

] = 24W d2x V^{-^aa + Ra) 



(4.46) 



Mathematically the fact that this action functional of the base metric g a b and the Weyl shift parameter 



a cannot be reduced to a single local functional of the full metric (4.44) means that the local Weyl group 



(2) 

of conformal transformations has a non-trivial cohomology, and ^ wz is a one-form representative of this 
cohomology [106, 1107] . This is just a formal mathematical statement of the fact that a effective action 
that incorporates the trace anomaly in a covariant EFT consistent with the Equivalence Principle must 
exist but that this S anom [g] is necessarily non-local. 

(2) 

It is straightforward in fact to find a non-local scalar functional Sanom[g] such that |105| 



r&fo a] = Sj» om [g = e 2 °g] - S®, 



(4.47) 



By solving (4.45) formally for a, and using the fact that W—gU 



-g\3 is conformally invariant in 



(2) _ 

two dimensions, we find that T)J 7 can be written as a Weyl shift (4.47) with 



wz 

2 



9(2) r i = 9. 

167T 



d 2 x- 



-g I d 2 x' 



-g'R(x) U' 1 {x,x')R(x l ) , 



(4.48) 



and □ _1 (x,x / ) denoting the Green's function inverse of the scalar differential operator □. The param- 
eter Q 2 is —N/6 if only matter fields in a fixed spacetime metric are considered. It becomes (25 — N)/6 
if account is taken of the contributions of the metric fluctuations themselves in addition to those of 
the N matter fields, thus effectively replacing N by N — 25 [108]. In the general case, the coefficient 
Q 2 is arbitrary, related to the matter central charge, and can be treated as simply an additional free 
parameter of the low energy effective action, to be determined. 



The anomalous effective action (4.48) is a scalar under coordinate transformations and therefore fully 



covariant and geometric in character, as required by the Equivalence Principle. However since it involves 
the Green's function □~ 1 (x,x / ), which requires boundary conditions for its unique specification, it is 
quite non-local, and dependent upon more than just the local curvature invariants of spacetime. In 
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this important respect it is quite different from the classical terms in the action, and describes rather 
different physics. In order to expose that physics it is most convenient to recast the non-local and 

(2) 

non-single valued functional of the metric, S an0 m into a local form by introducing auxiliary fields. In 



the case of (4.48) a single scalar auxiliary field, <p satisfying 



Dip = R 



(4.49) 



is sufficient. Then varying 

S a 2 nom[9'i ^] 



167T 



d z x ( g ab V b y> - 2R <p 



(4.50) 



with respect to (p gives the eq. of motion (4.49) for the auxiliary field, which when solved formally by 
ip = — U~ 1 R and substituted back into S a 2 n 0m [g; ip] returns the non-local form of the anomalous action 



(4.48), up to a surface term. The non-local information in addition to the local geometry which was 
previously contained in the specification of the Green's function D (x,x') now resides in the local 



auxiliary field <p(x), and the freedom to add to it homogeneous solutions of (4.49). 

;ric yi 

2 5Sanom[g;p] 



The variation of (4.50) with respect to the metric yields a stress-energy tensor, 

7 (2) 



-9 



Sg 



ab 



Q 2 

4-7T 



-V a V b (p + g ab Dp- ^(V a <p)(V b tp) + ^g ab (V c <p)(V c <p) 



(4.51) 



which is covariantly conserved, by use of (4.49) and the vanishing of the Einstein tensor, G ab = R, 



ab 



Rdab/2 = in two (and only two) dimensions. The classical trace of the stress tensor (4.51) is 

Q 2 



9 ab T^[g;p} = ^Cs 



Q 



4tt 



R 



(4.52) 



which reproduces the quantum trace anomaly in a general classical background (with Q 2 proportional 



to h). Hence (4.50) is exactly the local auxiliary field form of the effective action which should be added 



to the action for two dimensional gravity to take the trace anomaly of massless quantum fields into 
account. 

Since the integral of R is a topological invariant in two dimensions, the classical Einstein-Hilbert 

(2) 

action contains no propagating degrees of freedom whatsoever in d = 2, and it is Sanom which contains 



the only kinetic terms of the low energy EFT. In the local auxiliary field form (4.50), it is clear that 
5/1 



7 (2) 



describes an additional scalar degree of freedom ip, not contained in the classical action S\ . 
Once the anomalous term is treated in the effective action on a par with the classical terms, its effects 
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become non-perturbative and do not rely on fluctuations from a given classical background to remain 
small. 



Extensive study of the stress tensor (4.52) and its correlators, arising from this effective action estab- 
lished that the two dimensional trace anomaly gives rise to a modification or gravitational "dressing" 
of critical exponents in conformal field theories at second order critical points [108]. Since critical expo- 
nents in a second order phase transition depend only upon fluctuations at the largest allowed infrared 
scale, this dressing is clearly an infrared effect, independent of any ultraviolet cutoff. These dressed 
exponents and shift of the central term from N — 26 to N — 25 are evidence of the infrared fluctuations 
of the additional scalar degree of freedom tp which are quite absent in the classical action. The scaling 
dimensions of correlation functions so obtained are clearly non-perturbative in the sense that they are 
not obtained by considering perturbatively small fluctuations around flat space, or controlled by a uni- 
form expansion in A <C 1. The appearance of the gravitational dressing exponents and the anomalous 



effective action (4.48) itself have been confirmed in the large volume scaling limit of two dimensional 
simplicial lattice simulations in the dynamical triangulation approach |109[ IllUj . Hence there can be 



little doubt that the anomalous effective action (4.50) correctly accounts for the infrared fluctuations of 
two dimensional geometries. 

The importance of this two dimensional example is the lessons it allows us to draw about the role 
of the quantum trace anomaly in the low energy EFT of gravity, and in particular the new dynamics 
it contains in the conformal factor of the metric. The effective action generated by the anomaly in two 
dimensions contains a new scalar degree of freedom, relevant for infrared physics, beyond the purely 



local classical action. It is noteworthy that the new scalar degree of freedom in (4.49) is massless, 
and hence fluctuates at all scales, including the very largest allowed. In two dimensions its propagator 
x') is logarithmic, and hence is completely unsuppressed at large distances. This is the precise 
analog of the massless pole found in flat space triangle amplitudes in four dimensions. In d = 2 this 
pole appears already in two-point amplitudes of current correlators in the Schwinger model [91J, and 
in correlators of the energy momentum tensor of conformal fields [93]. Physically this means that 
the quantum correlations at large distances require additional long wavelength information such as 
macroscopic boundary conditions on the quantum state. 



The action (4.50) due to the anomaly is exactly the missing relevant term in the low energy EFT 
of two dimensional gravity responsible for non-perturbative fluctuations at the largest distance scales. 
This modification of the classical theory is required by general covariance and quantum theory, and 
essentially unique within the EFT framework. 
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D. The General Form of the Trace Anomaly in Four Dimensions 



The line of reasoning in d = 2 dimensions just sketched to find the conformal anomaly and construct 
the effective action may be followed also in four dimensions. In d = 4 the trace anomaly takes the 
somewhat more complicated form, 



(T a a ) =bF + b' (e-^Dr\ +b"nR + J2PiHi, 

^ ' i 

in a general four dimensional curved spacetime, where we employ the notation, 



(4.53) 



E =*R abcd *R abcd = RabcdR abcd ~ 4R ab R ab + R 1 , and (4.54a) 
with R a bcd the Riemann curvature tensor, *R a b c d = ^abefR e ^ c d/^ its dual, and C a b c d the Weyl conformal 



tensor. Note that E is the four dimensional Gauss-Bonnet combination whose integral gives the Euler 
number of the manifold, analogous to the Ricci scalar R in d = 2. The coefficients b, b' and b" are 



dimensionless parameters multiplied by K. Additional terms denoted by the sum (3iHi in (4.53) may 
also appear in the general form of the trace anomaly, if the massless conformal field in question couples to 
additional long range gauge fields. Thus in the case of massless fermions coupled to a background gauge 



field, the invariant H =tr(F a i,F ab ) appears in (4.53) with a coefficient (3 determined by the anomalous 
dimension of the relevant gauge coupling. 



As in d = 2 the form of (4.53) and the coefficients b and b' are independent of the state in which the 
expectation value of the stress tensor is computed, nor do they depend on any ultraviolet short distance 
cutoff. Instead their values are determined only by the number of massless fields [103, lllj, 

b = ]_20(47r) 2 {Ns + 6Nf + UNv) ' (455a) 
6' = -^ F (iV S + ^ + 62iV y) , (4.55b) 

with (Ns, Np, Ny) the number of fields of spin (0, |, 1) respectively and we have taken h = 1. Notice 
also that b > while b' < for all fields of lower spin for which they have been computed. Hence 
the trace anomaly can lead to stress tensors of either sign. The anomaly terms can be utilized to 
generate an effective positive cosmological term if none is present initially. Such anomaly driven inflation 



models |112j require curvatures comparable to the Planck scale, unless the numbers of fields in (4.55) is 
extremely large. It is clear that conformally flat cosmological models of this kind, in which the effects 
of the anomaly can be reduced to a purely local higher derivative stress tensor, are of no relevance to 
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the very small cosmological term (3.6) we observe in the acceleration of the Hubble expansion today. 



Instead it is the essentially non-local effects of the anomaly on the horizon scale, much larger than Lpi 



which should come into play. This requires a covariant action functional analogous to (4.50) for a proper 
treatment. This is what we now turn to computing. 

Three local fourth order curvature invariants E, F and dR appear in the trace of the stress tensor 



(4.53), but only the first two (the b and b') terms of (4.53) cannot be derived from a local effective 
action of the metric alone. If these terms could be derived from a local gravitational action we would 
simply make the necessary finite redefinition of the corresponding local counterterms to remove them 
from the trace, in which case the trace would no longer be non-zero or anomalous. This redefinition of a 
local counterterm (namely, the R 2 term in the effective action) is possible only with respect to the third 
b" coefficient in ( |4.53| ) , which is therefore regularization dependent and not part of the true anomaly. 



Only the non-local effective action corresponding to the b and b' terms in (4.53) are independent of 



the UV regulator and lead to effects that can extend over arbitrarily large, macroscopic distances. The 
distinction of the two kinds of terms in the effective action is emphasized in the cohomological approach 
to the trace anomaly [107 . 



To find the WZ effective action corresponding to the b and b' terms in (4.53), introduce as in two 



dimensions the conformal parameterization (4.44), and compute 



-gF = 
~g (e 



-gF 



■DR 



E 



DR +4V 3 ^A 4 cr, 



(4.56a) 
(4.56b) 



whose a dependence is no more than linear. The fourth order differential operator appearing in this 
expression is [ffiZl E31 033 



A 4 = D 2 + 2i? a6 V a V 6 - ^RO + ^(V a R)V a , 
which is the unique fourth order scalar operator that is conformally covariant, viz. 



(4.57) 



(4.58) 



for arbitrary smooth o~(x) in four (and only four) dimensions. Thus multiplying (4.53) by yj—g and 
recognizing that the result is the a variation of an effective action T\yz, we find immediately that this 
quadratic effective action is 

(4.59) 



T W z[g;o-\ = b \ (Fxyf^F 'a + b' \ tfxyf-t) \ \E- ^DR ) cr + 2o-A 4 cr 
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up to terms independent of a. This Wess-Zumino action is a one-form representative of the non- 
trivial cohomology of the local Weyl group in four dimensions which now contains two distinct cocycles, 



corresponding to the two independent terms multiplying b and b' |107j. By solving (4.56b) formally for 

(4.60) 



a, using (4.58), and substituting the result in (4.59) we obtain 

l [g = e 2a g]-S l 



Twz[g;&] = S a 



anom Li/ J i 



with the non-local anomalous action 



S a 



d x^fg \ d x \J g' 



E 
2~ 



UR 



A 4 1 (x,x / ) 



bF + b' 



E 
2~ 



UR 



(4.61) 



and A 4 1 (x,x / ) denotes the formal Green's function inverse of the fourth order differential operator 



defined by (4.57). From the foregoing construction it is clear that if there are additional Weyl invariant 
terms in the anomaly (4.53) they should be included in the S anom by making the replacement bF —> 



bF + in the last square bracket of (4.61). The case of the stress-energy of charged fermions 



coupled to photons in QED of Sec. |IV B| produces just such an additional term in the anomaly action, 
the effective action for which is the flat space limit of the general anomaly action. 



E. Anomaly Effective Action and Massless Scalar Fields 



As detailed in Ref. |115j we may render the non-local anomaly action (4.61 ) into a local form, by the 
introduction of two scalar auxiliary fields (p and ip which satisfy fourth order differential eqs., 

^ = \{e-\ur 



(4.62a) 

., - •• 'in- • .y^ - •■ ■ (4.62b) 

where we have added the last term in ( 4.62b[ ) to take account of the background gauge field. For the 
case of Dirac fermions, b = 1/3207T 2 , b' = — 11/5760-7T 2 , and c = — e 2 /24-/r 2 . The local effective action 



corresponding to (4.61) in a general curved space is given by 

= b'S^L + bSZL + \J ' d 4 x^-g F ap F 



S t , 



a/3 



If. 



(4.63) 



where 



SffinF \ J dV=s{-( □ ^) 2 +2 (rT- f <r)(V^)(V^)+ (e - 2 -URj ^ 



□ ^)(Dr)+2(^-|^)(V /( ^)(V,/r) 



(4.64) 
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The free variation of the local action (4.63)-(4.64) with respect to ip and ip yields the eqs. of motion 



(4.62). Each of these terms when varied with respect to the background metric gives a stress-energy 



tensor in terms of the auxiliary fields satisfying eqs. (4.62) 



If we are interested in only the first variation of the action with respect to g^, around flat spacetime, 



in order to compare to our calculation of the (TJJ) amplitude in Sec. IV B we may drop all terms in 



(4.63) which are second order or higher in the metric deviations from flat space. Since 5R is first order 



in variation around flat space, we may assume from (4.62a) that ip is as well. Then the entire b'S^ 



contribution to (4.63) is at least second order in this variation from flat space and cannot contribute 



to (TJJ). From (4.62b) the field ip has a contribution from F /1U F' J ' 1/ even in flat space, and a potential 
second order pole, LI -2 — > k~ 4 in its stress tensor. However, retaining only terms in (4.63) and (4.64) 



which contribute at first order in variation of the metric from flat space, we find that the second order 
pole is partially cancelled, and we obtain the simpler first order pole form, 



(4.65) 



which is valid to first order in metric variations around flat space, or 

R 



Sanom[g, A; if, Ip'] = j d^X 

its local equivalent, where 



-fj 



+ -F^F^<p 



(4.66) 



i// = bDip, 



□ p 



R 



(4.67a) 
(4.67b) 

(4.67c) 



Then after variation we may set <p = in flat space, and the only terms which remain in the stress 



tensor derived from (4.63) are those linear in ip', viz. 

THV>'0)]= 2 dS " 



-g bg^{ z ) 



(4.68) 



flat,(p=0 

which is independent of b and b' , and contain only second order differential operators, after the definition 



(4.67a). Solving (4.67b) formally for ip' and substituting in (4.68), we find 



(4.69) 



a result that may be derived directly from (4.65) as well 
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By varying (4.69) again with respect to the background gauge potentials, and Fourier transforming, 
we obtain 



r£E&(p,?) = / d A x J d 



4 ip-x+iq-y " '■"/">»' 

y & 



<J -LanomyJ) 



18vr 2 k 2 

which gives the full trace for massless fermions, 



5A a (x)A (y) 
1 (g^k 2 -m u ) u al3 (p,q), 



„ pV _ p pa/3 1 p pa0 

ilpiv-L anom ~ cr ap r — 247T 2 



(4.70) 



(4.71) 



in agreement with (4.53). The tree amplitude of the effective action (4.66) which reproduces the pole 



in the trace part of the (TJJ) triangle amplitude computed in Sec. IV B is illustrated in Fig. [7} 



T^[i/j f 
WvW-- 



FIG. 7: Tree Diagram of the Effective Action (4.66), which reproduces the trace of the triangle anomaly. The 
dashed line denotes the propagator G^i v = of the scalar intermediate state, while the jagged line denotes 

the gravitational metric field variation = Sg^. Compare to Fig. [6] 

The massless degrees of freedom ip and ip' are a necessary consequence of the trace anomaly, required 
by imposition of all the other symmetries. In this case these are scalar rather than pseudoscalar degrees 
of freedom. An important physical difference with the axial case is that the introduction in QED of a 
chiral current J§ and axial vector source corresponding to it appear rather artificial, and difficult 
to realize in Nature, whereas the trace of the stress tensor obtained by a conformal variation of the 
effective action is simply a particular metric variation already present in the QED Lagrangian in curved 
space, required by general coordinate invariance and the Equivalence Principle, without any additional 
couplings or extraneous fields. Since the stress-energy tensor couples to the universal force of gravity, we 
should expect that physical processes can excite the scalar ip and tp' scalar degrees of freedom required 
by the trace anomaly with a gravitational coupling strength, which can produce effects of arbitrarily 
long range. 
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F. The Scalar Anomaly Pole and Gravitational Scattering Amplitudes 

In order to verify the existence of the massless scalar pole in a physical process, consider the simple 
tree diagram of gravitational exchange between an arbitrary conserved stress-energy source T /fJ " u and 
photons illustrated in Fig. [8} 




FIG. 8: Tree Level Gravitational Scattering Amplitude 
This process is described by the scattering amplitude [116J, 




(4.72) 



The relative factor of — ^ between the two terms is dictated by the requirement that there be no scalar 
or ghost state exchanged between the two sources, and is exactly the prediction of General Relativity, 
linearized about flat space. That only a spin-2 propagating degree of freedom is exchanged between the 
two sources in Fig. [8] can be verified by introducing the following 3 + 1 decomposition for each of the 
conserved stress tensors, 

T 00 = Too , (4.73a) 
T 0i = _y±i _&}_foo, (4.73b) 




(4.73c) 



where diV^ 1 = 0, d{F^ % i = T\ % = 0, and V -2 denotes the static Green's function of the Laplacian 
operator, V 2 = d l di in flat space. This parameterization assumes only the conservation of the stress- 
tensor source(s), i.e. d^T^ = 0, so that there remain six independent components of T^ v which must 
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be specified, and we have chosen these six to be Tqq, V ±t , T^ 1 ^ and the total trace T^, which is a 



spacetime scalar. Substituting the decomposition (|4.73|) into (|4.72|) gives 
M = 8itG I d 4 x' I <I l .v 



rpfX 



2VJ 



1 



v 2 ^ 



1 



2 uu (v^,r 



□ T 00 + *T' 



00 



/ x.x v y x . 



'- 00 



which becomes the five terms 



M 



-8vrG 



T 



1 



1 



k 2 11 11 



(4.74) 



(4.75) 



in momentum space. These expressions show that only the spatially transverse and tracefree components 
of the stress tensor, Zy exchange a physical propagating helicity ±2 graviton in the intermediate state, 
characterized by a Feynman (or for classical interactions, a retarded) massless propagator — — > k~ 2 



pole in the first term of (4.74) or (4.75). All the other four terms in either expression contain only an 
instantaneous Coulomb-like interaction — V~ 2 — > k~ 2 or V~ 4 — > k~ 4 between the sources, in which 
no propagating physical particle appears in the intermediate state of the cut diagram. This is the 
gravitational analog of the decomposition, 



P: 



r = r 



of the conserved electromagnetic current and corresponding tree level scattering amplitude, 



dV [d*xJ ,fl (x')(— ) J„ 

J V n Jx'x 



(*) 



- V^— J 

k 2 



J, 



7± 



1 

k 2 



Jt+P' 



/I 

k 2 



(4.76a) 
(4.76b) 



(4.77) 



which shows that only a helicity ±1 photon is exchanged between the transverse components of the cur- 



rent, the last term in (4.77) being the instantaneous Coulomb interaction between the charge densities. 



We now replace one of the stress tensor sources by the matrix element of the one-loop anomalous 
amplitude, considering first the trace term with the anomaly pole in F%. This corresponds to the 
diagram in Fig. [9j We find for this term, 

(OlToob.g), = -k 2 F l (k 2 )u^(p 1 q)A a {p)A p {q) 



(OnMi = 3k 2 F 1 (k^u^(p,q)A a (p)A^q) 



\ T ii \P, 



0. 



Hence the scattering amplitude (4.75) becomes simply, 



Mi = ^GT\F 1 {k 2 )u a f j {p,q)A a {p)A^q) 



4vrG 



" k 2 



(4.78a) 
(4.78b) 
(4.78c) 



(4.79) 



69 



where (4.42) has been used for m = 0. Thus for massless fermions the pole in the anomaly amplitude 



becomes a scalar pole in the gravitational scattering amplitude, appearing in the intermediate state 
as a massless scalar exchange between the traces of the energy-momentum tensors on each side. The 
standard gravitational interaction with the source has produced an effective interaction between the 



scalar auxiliary field and the trace T ^ with a well defined gravitational coupling. Thus we may 
equally well represent the scattering as Fig. [9] involving the fermion triangle, or as the tree level diagram 
Fig. [To] of the effective theory, with a massless scalar exchange. 




FIG. 9: Gravitational Scattering of photons from the source T' Ml/ via the triangle amplitude 





FIG. 10: Gravitational scattering of photons from the trace of a source T ,fl via massless scalar exchange in the 



effective theory of (4.80) 



This tree diagram is generated by the effective action in flat space modified from (|4.66|) to 
Seff[g,A;(p,ip'} = / d A x^g 



(4.80) 



to include the coupling to the trace of the energy-momentum tensor of any matter T ^ source. Corre- 



spondingly the eq. (4.67c) for tp becomes 



(4.81) 
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instead of (4.67c). The eq. of motion for ip' remains (4.67b). We note that if the source T'^ v generates 



the curvature R by Einstein's eqs., then R = SirGT^, so that (4.80) and (4.81) are equivalent to 



(4.66) and (4.67c) at leading order in G. 

We conclude that in the conformal limit of massless electrons, the pole in the trace sector of the 
(TJJ) anomaly amplitude contributes to gravitational scattering amplitudes as would a scalar field 
coupled to the trace of the energy-momentum tensor of classical sources. The gravitationally coupled 
intermediate scalar can be understood as arising from collinear e + e~ correlated pairs in a total spin + 
state in Fig. [6j Although the result appears similar in some respects to a Brans-Dicke scalar |117j . and 
indeed ( 4.81[ ) is identical in form to that of a Brans-Dicke scalar with vanishing Brans-Dicke coupling 
(u) = 0) , because of the unusual off-diagonal kinetic term dictated by the structure of the trace anomaly, 
( |4.80 ) is certainly not a Brans-Dicke theory. Thus, although (4.81) tells us that tp is sourced by the 
trace of matter stress tensor with gravitational strength, ip cannot react back on matter except through 



i/)' and hence F a /3F al3 , reproducing (4.65), from whence it was derived. Each of the two scalar fields 
each couples to a different source, with an off-diagonal propagator, G^/^. There is no direct coupling of 
the trace of matter stress tensors T ^ to T u u via a scalar exchange as there would be in a classical scalar 



tensor theory of the Jordan-Brans-Dicke kind. Hence the phenomenology of (4.80) will be quite different, 
and the observational limits on a Jordan-Brans-Dicke scalar [5] do not apply. Note also that the matrix 
elements of F a pF a ^ = 2(E 2 — B 2 ) vanish for monochromatic photons on shell. Thus to leading order the 
scattering diagram in Fig. [9] also does not contribute to photon scattering on shell. The contribution 
of the massless scalar pole to higher order or off shell physical processes and the prospects for detecting 
its effects experimentally are important and interesting questions currently under investigation. 



G. The Effective Action of Low Energy Gravity 

With the foregoing detailed consideration of anomalies, massless poles and their long distance ef- 
fects, we consider finally the EFT of four dimensional macroscopic gravity. This gravitational EFT is 
determined by the same general principles as in other contexts |118j , namely by an expansion in powers 
of derivatives of local terms consistent with symmetry. Short distance effects are parameterized by 
the coefficients of local operators in the effective action, with higher order terms suppressed by inverse 
powers of an ultraviolet cutoff scale M uv . The effective theory need not be renormalizable, as indeed 
Einstein's theory is not, but is expected nonetheless to be quite insensitive to the details of the under- 
lying microscopic degrees of freedom, because of decoupling |118j . It is the decoupling of short distance 
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degrees of freedom from the macroscopic physics that makes EFT techniques so widely applicable, and 
which we assume applies also to gravity. 

As a covariant metric theory with a symmetry dictated by the Equivalence Principle, General Rela- 
tivity may be regarded as just such a local EFT, truncated at second order in derivatives of the metric 
field g a b( x ) [H9J • When quantum matter is considered, the stress tensor T a b becomes an operator. 
Because the stress tensor has mass dimension four, containing up to quartic divergences, the proper 
covariant renormalization of this operator requires fourth order terms in derivatives of the metric. How- 
ever the effects of such higher derivative local terms in the gravitational effective action are suppressed 
at distance scales L 3> Lpi in the low energy EFT limit. Hence surveying only local curvature terms, 
it is often tacitly assumed that Einstein's theory contains all the low energy macroscopic degrees of 
freedom of gravity, and that General Relativity cannot be modified at macroscopic distance scales, 
much greater than Lpi, without violating general coordinate invariance and/or EFT principles. As we 
have argued previously in two dimensions, this presumption should be re-examined in the presence of 
quantum anomalies. 

When a classical symmetry is broken by a quantum anomaly, the naive decoupling of short and 
long distance physics assumed by an expansion in local operators with ascending inverse powers of 
M uv fails. In this situation even the low energy symmetries of the effective theory are changed by 
the presence of the anomaly, and some remnant of the ultraviolet physics survives in the low energy 
description. An anomaly can have significant effects in the low energy EFT because it is not suppressed 
by any large energy cutoff scale, surviving even in the limit M vv — > oo. Any explicit breaking of 
the symmetry in the classical Lagrangian serves only to mask the effects of the anomaly, but in the 
right circumstances the effects of the non-local anomaly may still dominate the local terms. The axial 



anomaly in QCD, discussed in Sec. IV A has low energy effects, unsuppressed by the EFT ultraviolet 
cutoff scale, M uv ~ Aqcd in that case. Although the quark masses are non-zero, and chiral symmetry 
is only approximate in Nature, the chiral anomaly gives the dominant contribution to the low energy 
decay amplitude of ir° — > 2j in the standard model |81| [M] , a contribution that is missed entirely by 
a local EFT expansion in pion fields. Instead the existence of the chiral anomaly requires the explicit 
addition to the local effective action of a non-local term in four physical dimensions to account for its 
effects |104L 1118] . Although when an anomaly is present, naive decoupling between the short and long 
distance degrees of freedom fails, it does so in a well-defined way, with a coefficient that depends only 
on the quantum numbers of the underlying microscopic theory. 

The low energy effective action for gravity in four dimensions contains first of all, the local terms 
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constructed from the Riemann curvature tensor and its derivatives and contractions up to and including 
dimension four. This includes the usual Einstein-Hilbert action of General Relativity, 

Seh[9\ = / d ' x V=9(R- 2A) (4.82) 



16vrG . 

as well as the spacetime integrals of the fourth order curvature invariants, 

j( 4 ) _ I / ,. — 7, (,,/■• . .r 



Sfocai M = 2 J d x ^9 {aC abcd C abcd + (3R Z ) , (4.83) 
with arbitrary dimensionless coefficients a and (3. There are two additional fourth order invariants, 



namely E =*R a i, cc i*R and \3R, which could be added to (4.83) as well, but as they are total 
derivatives yielding only a surface term and no local variation, we omit them. All the possible local 
terms in the effective action may be written as the sum, 

1 f 00 
S local [9] = Y^n d ' X ^ 9{R - 2A) + S ^ + £ S iocli ■ ( 4 - 84 ) 

J n=3 

with the terms in the sum with n > 3 composed of integrals of local curvature invariants with dimension 
2n > 6, and suppressed by M~^ n+4 at energies much less than M uv . Here M uv is the ultraviolet cutoff 
scale of the low energy effective theory which we may take to be of order M v \. The higher derivative 
terms with n > 3 are irrelevant operators in the infrared, scaling with negative powers under global 
rescalings of the metric, and may be neglected at macroscopic distance scales, at least with respect to 
the classical scaling dimensions. On the other hand the two terms in the Einstein-Hilbert action n = 0, 1 



scale positively, and are clearly relevant in the infrared. The fourth order terms in (4.83) are neutral 
under such global rescalings. 

The exact quantum effective action also contains non-local terms in general. All possible terms 
in the effective action (local or not) can be classified according to how they respond to global Weyl 
rescalings of the metric, i.e. a = <jq = const. If the non-local terms are non-invariant under global 



rescalings, then they scale either positively or negatively under (4.44). If mT 1 is some fixed length 
scale associated with the non-locality, arising for example by the integrating out of fluctuations of fields 
with mass to, then at much larger macroscopic distances (mL S> 1) the non-local terms in the effective 
action become approximately local. The terms which scale with positive powers of e a ° are constrained 



by general covariance to be of the same form as the n = 0, 1 Einstein-Hilbert terms in Si oca i, (4.82). 
Terms which scale negatively with e a ° become negligibly small as mL 3> 1 and are infrared irrelevant 
at macroscopic distances. This is the expected decoupling of short distance degrees of freedom in an 
effective field theory description, which are verified in detailed calculations of loops in massive field 
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theories in curved space. The only possibility for contributions to the effective field theory of gravity at 



macroscopic distances, which are not contained in the local expansion of (4.84) arise from fluctuations 



not associated with any finite length scale, i.e. m = 0. These are the non-local contributions to the low 
energy EFT which include those associated with the anomaly. 

The anomaly effective action is associated with non-trivial cohomology of the Weyl group in the space 



of metrics. Since Twz (4.59) from which the effective action of the anomaly was derived satisfies WZ 



consistency, i.e. is closed but not exact under the Weyl group, it is unique up to an arbitrary admixture 
of local trivial cocycles, which in physical terms are either trivial because they are completely Weyl 
invariant effective actions obeying 



(4.85) 



and drop out of difference (4.60 ), or they are purely local terms easily catalogued by ascending powers of 



the Riemann curvature tensor, its covariant derivatives and contractions in (4.84). Thus the classification 



of terms according to their global Weyl scaling properties tells us that the exact effective action of any 
covariant theory must be of the form [107 , 



S exact [9] — Sl oca l{g] + Sini;[<7] + S anom [g] , 



(4.86) 



with Si oca i given by the expansion (4.84), Si nv the (generally non-local) Weyl invaraint terms satisfying 



(4.85), and S anom the anomaly action given by (4.63)-(4.64). The higher dimension local terms in (4.84) 



are strictly irrelevant in the IR, since they scale to zero with negative powers of e a ° and may be neglected 
for physics far below the Planck scale, while the lower dimension local terms are nothing but the terms 



of the usual Einstein-Hilbert classical action (4.82). These classical terms grow as positive powers of 



e a ° under global dilations and are clearly IR relevant terms. Indeed the naive classical scaling of these 
terms are positive powers (L and L 2 ) under rescaling of distance, and are clearly relevant operators of 
the low energy description. 

The local dimension four terms involving the Weyl tensor squared C 2 is fully locally Weyl invariant 
(while that involving R 2 is invariant only under global Weyl rescalings) and among the many terms that 
can appear in Si nv . Because both of these are neutral under global dilations, scaling like L° we expect 
them to be marginally irrelevant in the IR (as well as in the UV). All the higher dimension local terms 



in the sum in (4.84) for n > 3 scale to zero as <tq — > 00 and are clearly strictly IR irrelevant. However 



among all the possible terms that can be generated by quantum loops in the exact effective action of 
gravity, the anomaly effective action is unique in scaling logarithmically under the global Weyl group. 
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Indeed we note that in the form (4.64) the simple shift of the auxiliary field (p by a spacetime constant, 



tf + 2<7 



(4.87) 



corresponding to a global logarithmic variation of length scales yields the entire dependence of S a 



on the global Weyl rescalings (4.44), viz 



mi [g; <p, i>] ->■ Sanom 

[e ao g;ip + 2a ,ij} 

Sanom 

[g; (f, if)] + (Tq / (fx 



-fj 



bF + b' ( E - -UR 



owing to the strict invariance of the terms quadratic in the auxiliary fields under (4.44) and eqs. (4.56). 
Hence S anom scales logarithmically (~ log L) with distance under Weyl rescalings. 

Because of this infrared sensitivity to global rescalings, unlike local higher derivative terms in the 
effective action, which are either neutral or scale with negative powers of L, the anomalous terms 
should not be discarded in the low energy, large distance limit. Ordinarily, i.e. absent anomalies, the 
Wilson effective action should contain only local infrared relevant terms consistent with symmetry |120j . 
However, like the anomalous effective action generated by the chiral anomaly in QCD, the non-local 
Sanom must be included in the low energy EFT to account for the anomalous Ward identities, even 



in the zero momentum limit, and indeed logarithmic scaling with distance (4.88) indicates that S, 



is an infrared relevant term. Even if no massless matter fields are assumed, the quantum fluctuations 
of the metric itself will generate a term of the same form as S anom in the infrared [106]. The scalar 
fields of the local form (4.63)-(4.64) of S anom describe massless scalar degrees of freedom of low energy 



gravity, not contained in classical General Relativity. As we have seen, these massless scalars may be 
understood as correlated two-particle states of the underlying anomalous QFT, and show up as poles 
in gauge invariant physical scattering amplitudes in the EFT where the original quantum fields appear 
only in internal quantum loops. Thus the effective action of the anomaly S anom should be retained in 
the EFT of low energy gravity, which is specified then by the first two strictly relevant local terms of 

,, i.e. 



the classical Einstein-Hilbert action (4.82), and the logarithmic S, 



S eff[g] = S EH [g] + Sanom[g; <f, 1p] 



(4.89) 



with S anom is given by (4.61) or in its local form by ( 4.63 )-( 4.64). 



The complete classification of the terms in the exact effective action (4.86) into just three categories 
means that all possible infrared relevant terms in the low energy EFT, which are not contained in Si oca i 



of (4.84) must fall into S anom , i.e. they must correspond to non-trivial co-cycles of the local Weyl group 
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\107\ 1115] . The Weyl invariant terms in the exact effective action (4.85) are by definition insensitive 
to rescaling of the metric at large distances. Hence the (generally quite non-local) terms in Si nv do 
not give rise to infrared relevant terms in the Wilson effective action for low energy gravity. By this 
classification of terms (local or non-local) according to their behavior under global Weyl rescalings, the 



Wilson effective action (4.89) contains all the infrared relevant terms in low energy gravity for energies 
much less than M p i. 



Note also that it would be inconsistent with the semi-classical Einstein eqs. (2.42) to have as their 
source a stress tensor which is not covariantly conserved. Thus in a clash of symmetries which a 
quantum anomaly presents, it is necessary to choose the option that conformal invariance is broken, not 
general coordinate invariance. The action S an om is invariant under general coordinate transformations 
and under no conditions (i.e. in the presence of an horizon or not) is there a gravitational anomaly 
|89j . This is not a state dependent condition dependent on whether an horizon exists or not, but a 



condition of consistency of the semi-classical Einstein eqs. (2.42). Under the defining assumptions of 
general covariance and the EFT hypothesis of decoupling of physics associated with massive degrees 
of freedom, any infrared modifications of Einstein's theory generated by quantum effects is tightly 



constrained and the effective action (4.89) becomes essentially unique. The addition of the anomaly 
term(s) and the scalar degrees of freedom ip and ip they contain to the low energy effective action of 
gravity amounts to a non-trivial infrared modification of General Relativity, required by the existence 
of the trace anomaly, fully consistent with both quantum theory and the Equivalence Principle. 

V. MACROSCOPIC EFFECTS OF THE TRACE ANOMALY 

Having a fully covariant effective action the most straightforward application is to compute the 
covariantly conserved stress-energy tensor corresponding to it, and study its effects in particular back- 
grounds. If the effects are significant it will be necessary then to include the anomalous stress tensor 
as a source for Einstein's eqs. to obtain new solutions, but at first one can evauate the stress tensor 
in certain fixed backgrounds, such as the Schwarzschild black hole geometry discussed in Sec. IH] and 



the de Sitter geometry of Sec. III. This will provide the first evidence of the relevance of the anomaly 
effective action and the scalars tp and ip for macroscopic gravity in the presence of horizons. 

After these studies in fixed classical backgrounds one can consider next dynamical effects of the 
fluctuations associated with the anomaly scalars, their role in cosmology and their relevance to both 
the problem of gravitational collapse and cosmological dark energy. 
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A. Anomaly Stress Tensor in Schwarzschild Spacetime 



From (4.63) and (4.64) the stress tensor of the anomalous effective action consists of two independent 
terms, 

2 fiSanom i 



namely 



and 



p(anom) 
l ab 



-9 $9' 



ab 



: [g;<p,i>] = b'E ab + bF ab 



E a b 



2 (y, a<P ) (v 6) □ <p) + 2 v c [(y c<P ) (v a v b ip)] - - v a v b [( v c¥ >) (v V)] 



+ ; 



(V^)(VV) -4i? c (a (V b) ^)(V c vp) + -R(V a ip)(V b <p) 



+ 



^g ab { -3 ( □ if) 2 + □ [(Vc¥>) ( V V)] + 2 (3i? cd - i?<f d ) (V c¥ >) (V d <p) } 
- - V a V b Dip -4 C a \ d V c V d y? - 4 V 6) V c ^ + ^R ab □ y> + - R V a V b <p 
{V {a R) V b) ip + ^g ab {2D V + 6i? cd V c V d¥ > - 4i?n^ + (V C «)V C ^} 



F afe = -2(V (a ^)(V 6) n^)-2(V (a ^)(V b) n<^) + 2V c [(V c ^)(V a V^) + (V c ^)(V a V^)] 
■^V.V, [(V C ^)(V C ^)] + i 



-^V a V 6 [(V^)(V C ^)] + ^ a6 (V c ^)(V c ^) -4ir ( J(V b) ^)(V c V0 + (V^)(V c ^)] 



+ ^R(V (a cp)(V b) iP) + l -g ab [ - 3(D^)(DV) + □ [(V C ^)(V C V)] 
+2 (sf^ - i? 5 cd ) (V^)(V d V)} - 4 V C V, (c (a c b) <V) - 2C a c ^ c V 



Each of these two tensors are individually conserved and they have the local traces, 



(5.1) 



(5.2) 



(5.3) 



E\ = 2A 4 ip = E--aR, 
F a a = 2A 4 V = F = C abcd C 



abed 



(5.4a) 
(5.4b) 



corresponding to the two terms respectively in the trace anomaly in four dimensions in (4.53) (with 
Pi = 0). 

In the four dimensional Schwarzschild geometry the full contraction of the Riemann tensor is given 



by (2.6), and since R ab = the invariants F\ = E\ s =48 M jr . (where we temporarily set GM = M 
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to simplify the expressions below). A particular solution of either of the inhomogeneous Eqs. (4.62) is 
given then by <p(r), with 



dip 
dr 



AM r v \ 



1 

2M 



1 



AM 
r 



(5.5) 



and f(r) 



1 



2M 



The general solution of (4.62) for <p = <p(r) away from the singular points r 



(0, 2M, 00) is easily found and may be expressed in the form 

q-2 



dip 
dr 



dip 
s dr 
q — 2 ( r 



2Mc„ 
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(5.6) 



in terms of the three dimensionless constants of integration, c H , Cqo, and q. This expression has the 
limits, 



dip 
dr 



dip 
dr 



s 



r - 2M 

2 Cqo 



+ 



q - 2 / r 
2Af n V2M 



1 



1 

2M 



3c 



+ 



q Coo 

' 4M 2 ' AM r 3r 2 "* ~ V2iW7 

Hence c H controls the leading behavior as r approaches the horizon, while Coo controls the leading 
behavior as r — > 00, which is the same as in flat space. The leading behavior at the horizon is determined 



2M / r \ 2M 



q In 



+ 



7 

18 



(9-2) 



+ ... 



2M; (5.7a) 
^00. (5.7b) 



by the homogeneous solution to (4.62), c H In / = c H \n{—K a K a )^ where K = dt is the timelike Killing 



field of the Schwarzschild geometry for r > 2M . Because of these singular behaviors in (5.7b), (5.6) is 



clearly a solution of (5.4) in the distributional sense, i.e. containing possible 5 functions or derivatives 



thereof at the origin, at the horizon and at infinity. 



To the general spherically symmetric static solution (5.6) we may add also a term linear in t, i.e. we 
replace ip(r) by 



<p(r,t) = <p(r) + j M ' t - 



(5.8) 



Linear time dependence in the auxiliary fields is the only allowed time dependence that leads to a time- 
independent stress-energy. The conformal transformation to flat space near the horizon corresponds to 
the particular choice c H = ±p = 1, leaving the subdominant terms in (5.7a) parameterized by q and Coo 
undetermined. 

Since the equation for the second auxiliary field ip is identical to that for 93, its solution for tp = tp(r, t) 



is of the same form as (5.6) and (5.8) with four new integration constants, d H ,d oa , q' and p' replacing 
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c H , c oo, Q, and p in cp(r, t). Adding terms with any higher powers of t or more complicated t dependence 



produces a time dependent stress-energy tensor. Inspection of the stress tensor terms in ( 5.1 ) also shows 
that it does not depend on either a constant <po or ipo but only the derivatives of both auxiliary fields 
in Ricci flat metrics such as Schwarzschild spacetime. For that reason we do not need an additional 



integration constant for either of the fourth order differential equations (4.62). 



With the general spherically symmetric solution for (p(r,t) and ^(r,t), we can proceed to compute 



the stress-energy tensor (5.1) in a stationary, spherically symmetric quantum state. For example the 
Boulware state is characterized as that state which approaches the flat space vacuum as rapidly as 
possible as r — > oo [35]. In the flat space limit this means that the allowed r 2 and r behavior in the 
auxiliary fields (r and constant behavior in their first derivatives) must be set to zero. Inspection of the 



asymptotic form (5.7b) shows that this is achieved by requiring 



Coo = doo = and (5.9a) 
q = q = . (Boulware) (5.9b) 

If we set p = p' = as well, in order to have a static ansatz for the Boulware state, then the remaining 
two constants c H and d H are free parameters of the auxiliary fields, which lead to a stress-energy 
which diverges as r — > 2M on the horizon. Matching the leading divergence of the stress-energy in the 



Boulware state (2.41) by adjusting c„ and d H appropriately, one then has a one parameter fit to the 



numerical data of |121j . The results of this fit of (5.1) with c H and d H treated as free parameters are 



illustrated in Figs. [IT]to[13| for all three non-zero components of the stress tensor expectation value of 
a massless, conformally coupled scalar field in the Boulware state. The (approximate) best fit values 
plotted were obtained with c H = — ^ and d H = ||. 

For comparison purposes, we have plotted also the analytic approximation of Page, Brown, and 



Ottewill [122T I124] (dashed curves in Figs. 11 to 13). We observe that the two parameter fit with the 



anomalous stress tensor in terms of the auxiliary ip and t/j fields is more accurate than the approximation 
of refs. j!22H124j for the Boulware state. 

The next important point to emphasize is that the stress-energy diverges on the horizon in an en- 
tire family of states for generic values of the eight auxiliary field parameters (c H , q, Coo,p; d H , q' , doo,p'), 
in addition to the Boulware state. Hence in the general allowed parameter space of spherically sym- 
metric macroscopic states, horizon divergences of the stress-energy are quite generic, and not resticted 
to the Boulware state. On the other hand, the condition that the stress-energy on the horizon be 
finite gives four conditions on these eight parameters, in order to cancel the four possible divergences 
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FIG. 11: The expectation value (T t *) of a conformal scalar field in the Boulware state in Schwarzschild spacetime, 
as a function of s = r -^- in units of 7r 2 T^/90. The solid curve is Eq. (^-IpTl with ( pTHj ), dsTg] ) , ( jHTo] ) and 
c H = —^Q,d H — || . The dashed curve is the analytic approximation of |122j , and the points are the numerical 
results of [T2"T] . 
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-2 



-3 



-4 




FIG. 12: The radial pressure (T r r ) of a conformal scalar field in the Boulware state in Schwarzschild spacetime. 
The axes and solid and dashed curves and points are as in Fig. [TTJ Note the better agreement of the anomaly 
stress tensor (solid curve) with the numerical results of |121| . compared to the analytic approximation of [122] . 
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FIG. 13: The tangential pressure (T s ) of a conformal scalar field in the Boulware state in Schwarzschild spacetime. 
The axes and solid and dashed curves and points are as in Fig. [TT] 

s~ 2 , s _1 , In 2 s, In s. The simplest possibility with the minimal number of conditions on the auxiliary 
field parameters is via 



(26 + b')c 2 + p(2bp + b'p) = 



(In s and Ins) 



(5.10a) 
(5.10b) 
(5.10c) 



It is clear that these finiteness conditions on the horizon are incompatible with the conditions for fall 



off at infinity (5.9), in their values of q and q' . Thus the effective action and behavior of the anomaly 



scalar fields and stress tensor illustrate at a glance the general topological obstruction discussed in Sec. 
[I!]for vacuum behavior of (T a b ) R at both the horizon and infinity simultaneously. With conditions (|5.9l) 
the anomaly scalar fields give a very weak Casimir-like long-range interaction between massive bodies, 
that falls off very rapidly with distance (at least as fast as r~ 7 ). 

By taking different values of the parameters (c H ,q, Coo,p; d H ,q' , d^^p') of the homogeneous solutions 



to (4.62), corresponding to different states of the underlying quantum theory, it is possible to approxi- 



mate the stress tensor of the Hartle-Hawking and Unruh states as well |115|, I125j . In states which are 
regular on the horizon there is no particular reason to neglect the Weyl invariant terms Si nv in the exact 
effective action, and the stress-energy tensor it produces would be expected to be comparable in magni- 
tude to that from S anom . In this case of bounded stress tensors, the contributions from both S anom and 
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Si nv are both of order M 4 and negligibly small on macroscopic scales in any case. However in states 
such as the previous Boulware example, the diverging behavior of the stress tensor near the horizon is 



captured accurately by the terms in (5.1) arising from the anomaly, which have the same generically 
diverging behaviors as the quantum field theory expectation value (T a b ) R . This general behavior has 
been found in the Reissner-Nordstr0m case of electrically charged black holes as well [126J. We find no 
a priori justification for excluding generic states with stress-energy tensors that grow without bound 
as the horizon is approached. In any such states the backreaction of the stress-energy on the geometry 
will be substantial in this region and lead to large backreaction effects near the horizon. 



B. Anomaly Stress Tensor in de Sitter Spacetime 



In conformally flat spacetimes with g a b = e 2a rj a b, one can choose tp = 2a and ip = to obtain 
the stress tensor of the state conformally transformed from the Minkowski vacuum. In this state F a b 
vanishes, and tp can be eliminated completely in terms of the Ricci tensor with the result |107| 1127] . 

2 _ „ „ d 5 



2 c 14 

E a b = -VaVfci? + 2R a c R bc — —RR a b + g a b ( —-OR - R c dR ca + g-R' 



(5.11) 



Thus all non-local dependence on boundary conditions of the auxiliary fields tp and tp drops out in 
conformally flat spacetimes for the state conformally mapped from the Minkowski vacuum. In the 
special case of maximally 0(4, 1) symmetric de Sitter spacetime, R a b = 3H 2 g a h with R = 12H 2 a 
constant, 



E, 



ab 



dS 



SH^Qab. 



(5.12) 



Hence we obtain immediately the expectation value of the stress tensor of a massless conformal field of 
any spin in the Bunch-Davies (BD) state in de Sitter spacetime, 



Tab 



BD,dS 



66' 



'H 4 g a b = ~^2 9ab (N s + llN f + 62N V ) 



H 



(5.13) 



which is determined completely by the trace anomaly. 

States in de Sitter space which are not maximally 0(4, 1) symmetric are easily studied by choosing 



diiferent solutions of (4.62) for tp and ip. For example, if de Sitter spacetime is expressed in the spatially 



flat coordinates (3.20) and tp = tp(r), we obtain from (5.1) the stress-energy in spatially homogeneous, 
isotropic states. Since in de Sitter spacetime the A4 operator factorizes, 



dS 



'□ -2H 2 ) (Dtp) = 12H 4 



(5.14) 
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it is straightforward to show that the general solution to this equation with <p = (f(r) is 



ip(r) = 2Ht + c + c_ie Hr + c_ 2 e ' itlT + c_ 3 e 



-2Ht 



-3Ht 



(5.15) 



When the full solution for of (5.15) is substituted into (5.2) we obtain additional terms in the 
stress tensor which are not de Sitter invariant, but which fall off at large r, as e~ 4HT . The stress tensor 
of this time behavior is traceless and corresponds to the redshift of massless modes with the equation 
of state, p = p/3. 

States of lower symmetry in de Sitter spacetime may be found by considering static coordinates 



(3.27), in terms of which the operator A4 again separates. Inserting the ansatz if = (p(r) in (5.14), the 



general 0(3) spherically symmetric solution regular at the origin is easily found: 

l-Hr\ 2c„-2-o, /1-Hr 



ip(r) 



dS 



In (1 - H 2 r 2 ) + c 



In i t- I + ~" g „„ 2 g ln 



(5.16) 



2"'\l + HrJ 2Hr \l + Hr / 

A possible homogeneous solution proportional to 1/r has been discarded, since it is singular at the origin. 
An arbitrary linear time dependence 2Hpt could also be added to 9?(r), i.e. (p(r) — > (p(r, t) = ip(r)+2Hpt. 
The particular solution, 



tp BD ( r , t) = In (1 - H 2 r 2 ) + 2Ht = 2Hr . 



(5.17) 



is simply the previous solution for the Bunch-Davies state we found in the homogeneous flat coordinates 



(3.20). 



In the general 0(3) symmetric state centered about the origin of the coordinates r = 0, the stress 
tensor is generally dependent on r. In fact, it generally diverges as the observer horizon r = H v is 



approached, much as in the Schwarzschild case considered previously. From (5.16) we observe that 

1 - Hr 



dS 



Ch + K-1-|) (1-Hr) + ... 



hi 



+ 0(1 -Hr) 



(5.18) 



as Hr — > 1, so that the integration constant c H controls the most singular behavior at the observer 
horizon r = H -1 . 

The second auxiliary field ifi satisfies the homogeneous equation, 



A 4 V> = , 

which has the general spherically symmetric solution linear in t 

Hr 



i/j(r,t) 



d Q + 2Hp't+\\n(\—- 
2 \ 1 + Hr 



(5.19) 



2d H ~ Q 1 
H — Mn 

2Hr 



1 - Hr 
1 + Hr 



(5.20) 
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Note that the constant d H enters this expression differently than c H enters the corresponding Eq. (5.16 ), 



due to the inhomogeneous term in (5.14), which is absent from the ifi equation. Since the anomalous 



stress tensor is independent of Co and do, it depends on the six parameters (c H ,d H ,q,q' ,p,p') in the 
general stationary 0(3) invariant state. The simplest way to insure no r~ 2 or r~ l singularity of the 
stress tensor at the origin is to choose q = q' = 0. With q = q' = there are no sources or sinks at the 
origin and the zero flux condition, 

4# 2 



(bpq' + bp q + b'pq) = 



(5.21) 



is satisfied automatically, for any p and p' . Because of the subleading logarithmic behavior (5.18) of ip 
on the observer horizon at r = H^ 1 , there will ln 2 (l — Hr) and ln(l — Hr) divergences in the other 



components of the stress tensor (5.1) in general. These divergences are removed when q = q' = and 



1 . All the divergences of the stress-energy at both the origin and the observer horizon are cancelled 



if the four conditions, 



2b (d H +pp') = b' (1 -p 2 ) 
ch = 1 and 
a = a' = 



(5.22a) 
(5.22b) 
(5.22c) 



are satisfied. These are satisfied by the Bunch-Davies state with p = ±1 and d = p' = 0. In fact, with 



conditions (5.22) the tensor F a b of (5.3) vanishes identically, the ip field drops out entirely, and the full 



anomalous stress tensor is given by (5.2), which takes the Bunch-Davies form (5.13). 



If the first two of the four conditions (5.22) are relaxed, then the stress-energy remains finite at the 



origin, but becomes divergent at r = H~ l . This is the generic case. It includes in particular the analog 
of the static Boulware vacuum [35] in de Sitter space. There is no analog of the Unruh state |38| in 
de Sitter spacetime, since by continuity a flux through the future or past observer horizon at r = H^ 1 



would require a source or sink of flux at the origin r = 0, a possibility we have excluded by (5.21 ) above. 



The important conclusion from these studies of the stress tensor (5.1) obtained from the anomaly 



effective action is that rather than being a pathology of a single state the divergent behavior of the stress 
tensor on the horizon is generic. The general spherically symmetric solutions of the linear scalar field 



eqs. (4.62) are easy found, which allows an overview of a variety of quantum states in the underlying 



field theory which would normally require a laborious effort to study individually, and for each spin 



field separately. This overview of possible states shows that at least four conditions (5.10) on the set 



of eight parameters is necessary to eliminate divergences on the future horizon of a Schwarzschild black 
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hole, and a similar set of conditions (5.22) are necessary to remove divergences on the de Sitter horizon. 
This suggests that in the generic state of the gravitational collapse problem the backreaction effects will 
be large, and alter significantly the classical picture of a black hole horizon, which is the source of the 
paradoxes discussed in Sec. |IT| and that the black hole and cosmological horizon singularities may be 
related. 

C. Conformal Phase of 4D Gravity and Infrared Running of A 

In order to understand the dynamical effects of the kinetic terms in the anomaly effective action, one 



can consider simplest case of the quantization of the conformal factor in the Wess-Zumino action (4.59) 
in the case that the fiducial metric is flat, i.e. g a b = e 2r7 r] a b. Then the Wess Zumino effective action 



(4.59) simplifies to 

r W z [m *] = - / <^ x ( D a ) 2 ' ( 5 - 23 ) 

where 

Q 2 = -32n 2 b' . (5.24) 

This action quadratic in a is the action of a free scalar field, albeit with a kinetic term that is fourth 
order in derivatives. The propagator for this kinetic term is (p 2 )~ 2 in momentum space, which is a 
logarithm in position space, 

G a {x,x') = -~\n[ l i 2 {x-x') 2 ] , (5.25) 

the same behavior as that for a massless field in two dimensions. Of course this is no accident but 
rather a direct consequence of the association with the anomaly of a conformally invariant differential 
operator, □ in two dimensions and A4 in four dimensions, a pattern which continues in all higher even 
dimensions. Because of this logarithmic propagator a similar sort of infrared fluctuations, conformal 
fixed point and dressing exponents as those obtained in two dimensional gravity |114| . 
The classical Einstein-Hilbert action for a conformally flat metric g a b = e 2a r] a b is 

-!- f d 4 x [3e 2 °(d a a) 2 -Ae 4 °] , (5.26) 

07TCx J 

which has derivative and exponential self-interactions in a. It is remarkable that these complicated 



interactions can be treated systematically using the the fourth order kinetic term of (5.23). In fact, these 
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interaction terms are renormalizable and their anomalous scaling dimensions due to the fluctuations of 
a can be computed in closed form |114l 11281 H29j . Direct calculation of the conformal weight of the 
Einstein curvature term shows that it acquires an anomalous dimension fii given by the quadratic 
relation, 



2 + 



(5.27) 



2Q 2 • 

In the limit Q 2 — > oo the fluctuations of a are suppressed and we recover the classical scale dimension 



of the coupling G 1 with mass dimension 2. Likewise the cosmological term in (5.26) corresponding to 
the four volume acquires an anomalous dimension given by 

Pi 



A) = 4 + 



(5.28) 



2Q 2 " 

Again as Q 2 — > oo the effect of the fluctuations of the conformal factor are suppressed and we recover 



the classical scale dimension of A./G, namely four. The solution of the quadratic relations (5.27) and 



(5.28) determine the scaling dimensions of these couplings at the conformal fixed point at other values 
of Q 2 . This can be extended to local operators of any non- negative integer mass dimension p, with 



associated couplings of mass dimension 4 — p, by 



(3 p = 4-p + 



2Q 2 



In order to obtain the classical scale dimension 4 — p in the limit Q 2 
is determined so that 



P P = Q 2 



l 



(8 - 2p) 

Q 2 



(5.29) 

oo the sign of the square root 

(5.30) 



valid for Q 2 > 8 — 2p for all p > 0, and thus Q 2 > 8. These scaling dimensions were computed both by 
covariant and canonical operator methods. In the canonical method we also showed that the anomalous 
action for the conformal factor does not have unphysical ghost or tachyon modes in its spectrum of 
physical states [130J. 



In the framework of statistical mechanics and critical phenomena the quadratic action ( 5.23 ) describes 
a Gaussian conformal fixed point, where there are no scales and conformal invariance is exact. The 
positive corrections of order 1/Q 2 (for Q 2 > 0) in (5.27) and (5.28) show that this fixed point is stable 
in the infrared, that is, both couplings G^ 1 and A/G flow to zero at very large distances. Because both 
of these couplings are separately dimensionful, at a conformal fixed point one should properly speak 
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only of the dimensionless combination hGA/c 5 = A. By normalizing to a fixed four volume V = J d 4 x 
one can show that the finite volume renormalization of A is controlled by the anomalous dimension, 



28 - 1 = 2 



1 + a/1 



< 0. 



(5.31) 



This is the anomalous dimension that enters the infrared renormalization group volume scaling relation 

d 



V-X = H25-1)X. 



(5.32) 
0.414 at 



The anomalous scaling dimension (5.31) is negative for all Q 2 > 8, starting at 1 — y/2 - 
Q 2 = 8 and approaching zero as —1/Q 2 as Q 2 — > oo. This implies that the dimensionless cosmological 
term A has an infrared fixed point at zero as V — > oo. Thus the cosmological term is dynamically driven 



to zero as V — > oo by infrared fluctuations of the conformal part of the metric described by (5.23). 



We emphasize that no fine tuning is involved here and no free parameters enter except Q 2 , which 
is determined by the trace anomaly coefficient b' by (5.24). Once Q 2 is assumed to be positive, then 



25 — 1 is negative, and A is driven to zero at large distances by the conformal fluctuations of the metric, 
with no additional assumptions. 



The result (5.32) does rely on the use of (5.23) or its curved space generalization (4.61) as the free 



kinetic term in the effective action for gravity, treating the usual Einstein-Hilbert terms as interactions or 
marginal deformations of the conformal fixed point. This conformal fixed point represents a new phase 
of gravity, non-perturbative in any expansion about flat space. In this phase conformal invariance is 
restored and the mechanism of screening A due to quantum effects proposed in |114j is realized. 

Identifying the fluctuations responsible for driving A to zero within a framework based on quantum 
field theory and the Equivalence Principle, free of ad hoc assumptions or fine tuning is an important first 
step towards a full solution of the cosmological constant problem. However, this is not yet a complete 
or testable cosmological model. Near the conformal fixed point the inverse Newtonian constant G^ 1 is 
also driven to zero when compared to some fixed mass scale m |129j . This is clearly different from the 



situation we observe in our local neighborhood. Under what conditions can (4.61) play a decisive role 



in realistic cosmological models? How would the conformal fixed point behavior of the scaling relations 
such as (5.32) be detectable by cosmological observations? 
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VI. COSMOLOGICAL CONSEQUENCES OF THE ANOMALY EFFECTIVE ACTION 



Absent a full cosmological model, the simplest possibility is to consider the essentially kinematic ef- 
fects of conformal invariance in cosmology and predictions for observables such as the Cosmic Microwave 
Background (CMB), without specifying the dynamics or spacetime history of the universe which gave 
rise to that conformal invariance. Later in this section we take a first step to a cosmological model by 
considering the semi-classical linear response of quantum fluctuations generated by the anomaly terms 
in the effective action around de Sitter space. 



A. Conformal Invariance in de Sitter Space 

To illustrate how conformal invariance on flat spatial sections of FLRW cosmology can be a natural 
result of a de Sitter phase in the evolution of the universe, consider the 0(4, 1) de Sitter symmetry 



group and its ten Killing vectors, £a . In coordinates (|3. 20h . the Killing eq., 



Va^ a) +V^i Q) = 0, a = l,...,10, (6.1) 

becomes 

d T Ci a) =0, (6.2a) 

d T 6* ] + dS a) - 2H^ a) = , (6.2b) 

ditjf + d^\ a) - 2Ha\ j ^ = . (6.2c) 

For £ T = we have the three translations, a = Tj, 

6 Tj) =0, ^=a 2 5{, j = 1,2,3, (6.3) 

and the three rotations, a = 

e^)=0, g Rt) =a 2 e l£m x m , £ = 1,2,3. (6.4) 
This accounts for 6 of the 10 de Sitter isometries which are self-evident in the flat FLRW coordinates 



with £ r = 0. The 4 additional solutions of (6.2) have £ T 7^ 0. They are the three special conformal 
transformations of M 3 , a = C n , 

f (C n ) = _2Hx n , eP n) = H 2 a 2 {8 n i 5 jk x j x k - 25 ij x j x n ) - 8% n = 1,2, 3, (6.5) 



S8 



and the dilation, a = D, 



6 D) = h £ D) n<r«, r i J - 

This last dilational Killing vector is the infinitesimal form of the finite dilational symmetry, 

x — > Ax, 
77 — >• A?7 

a(r) ->■ A~ 1 a(r) , 
r — )• r — H^ 1 In A 



(6.6) 



(6.7a) 
(6.7b) 
(6.7c) 
(6.7d) 



of de Sitter space. The existence of this symmetry explains why Fourier modes of different \k\ leave 
the de Sitter horizon at a shifted FLRW time r, so in an eternal de Sitter space, in which there is no 
preferred r, one expects a scale invariant spectrum. 



The existence of the three conformal modes of R (6.5) implies in addition that any 0(4, 1) de 



Sitter invariant correlation function must decompose into representations of the conformal group of 
three dimensional flat space. Fundamentally this is because the de Sitter group 0(4, 1) is the conformal 



group of flat Euclidean M , as eqs. (6.2)-(6.6) shows explicitly. Moreover, because of the exponential 



expansion in de Sitter space, the decomposition into representations of the conformal group become 
simple at distances large compared to the horizon scale 1/H. 

If we consider a massive scalar field in the Bunch-Davies state in de Sitter space (with m 2 = ii 2 +£,R = 
fj? + 12^H 2 ), its propagator function G(z(x,x');m 2 ) satisfies 



(-□ +m 2 ) G(z(x,x');m 2 ) 



d 2 , . d m 2 

^-^ + 2{1 - 2 ^-jp 



G{z;m 2 ) = 5 4 (x,x'), 



(6.8) 



where z(x, x') is the de Sitter invariant distance function given by 

l-z(x,x') = ^ rf [-( V - V ') 2 + (x-x') 2 } 
H 2 

= —ra{r])a(i]') [-(77 - rj) 2 + (x - x) 2 ] , 



(6.9) 



in the conformally flat coordinates of (3.24). Since for x 7^ x' , (6.8) is a standard form of the hyperge 
ometric eq., it has the solution, 

H 2 



G(z;m 2 



16vr 2 



r(a)r(/3)F(a,/3;2;z) 



.10) 
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in terms of the hypergeometric function 2F1 = F, with parameters 



(6.11a) 
(6.11b) 



where v is defined by (3.45). The particular solution (6.10) and its normalization is dictated by the 



requirement that (— □ + m 2 )G(x, x') yield a delta function of unit magnitude at x = x' in (6.8) and 
possess no other singularities. 



To explicitly see the conformal behavior in the massive case consider (6.10) in the limit a(r), a(r') 3> 



1, that is, after many e-foldings of the exponential de Sitter expansion [131] . From (6.9) in this limit, 



H 2 

1 — z{x,x) — > -^-a(T)a(T ! ) (x — x) 2 — > 00 



(6.12) 



for fixed x,x'. Using the asymptotic form of the hypergeometric function in (6.10) for z — > —00, we 
obtain [131J 



G(z(x,x');M 2 ) ->• A + \x - x 



-t ^ 1-3+21/ 



+ A- \x — x 



-3-2v 



(6.13) 



a sum of simple power law scaling behaviors in the distance \x — x'\, with coefficients 



1 



H 



±2v 



r(§=Fi/)r(±2i/) 



exp 



±u)H(t + t') 



(6.14) 



2Htt 2 v 2 ) T(i ± v) 

Thus even massive non-conformal fields in de Sitter space exhibit conformal power law scaling behavior 
on the flat M 3 spatial sections after sufficiently long exponential expansion in a de Sitter inflationary 



phase. The corrections to (6.13) are integer power law terms which are a Taylor series in z ^1 given 



by (6.12). Moreover if m 2 < 9H 2 /4 then only the leading power law term \x — x'| _3+2l/ survives in the 
limit. This would identify the CFT scaling exponent to be 3/2 — v if the fluctuations responsible for 
the CMB were generated in the late time de Sitter expansion by a field of mass m such that u defined 



in (3.45) is real. The conformal behavior of de Sitter invariant correlation functions is an example of 



a kind of dS/CFT correspondence [107] , and the result of the mathematical isomophism between the 
conformal group of the flat M 3 sections and the four-dimensional de Sitter group SO(4, 1) expressed by 
the solutions of the Killing equation. 

Note that in this picture the fluctuations of a massless, minimally coupled scalar field with m = 
or gravitons themselves have v = 3/2 and A = u — 3/2 = 0, as should be expected for the fluctuations 
of a dimensionless variable such as the metric tensor, and a nearly conformally invariant spectrum is 
generated on the FLRW flat spatial sections. This presumes that the de Sitter isometry group is broken 
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in a specific way, to a standard spatially homogeneous, isotropic FLRW cosmology. In slow roll inflation 
this breaking of de Sitter invariance is achieved by the addition of a scalar field, the inflaton, with a 
nearly flat potential [132[ 1133] . The fluctuations of this field generate fluctuations in its stress tensor 
which act as a source for scalar metric fluctuations through Einstein's eqs. as the system rolls slowly 
down the potential and out of de Sitter space everywhere in space uniformly, maintaining overall spatial 



homogeneity and isotropy and the form (3.20) of the FLRW line element. Since gravity is treated purely 
classically, in the limit of a strictly constant inflaton potential or in the absence of any inflaton field at 
all there would be no coupling of the scalar inflaton modes to the metric and hence the amplitude of the 
CMB power spectrum would vanish. For this reason also in this picture the non-Gaussian bi-spectrum 
is expected to have a very small amplitude, higher order in the slow roll parameters [134J. 

Now the effective action of the trace anomaly has the potential to modify this picture in at least two 
fundamental ways. First, the scalar degrees of freedom (p and ip are present in the effective action and 
will generate scalar fluctuations in a de Sitter epoch, independently of any inflaton field or potential. 
Thus, the deviations from a strict Harrison- Zel'dovich spectrum and primordial non-Gaussianities need 



not be controlled by any small slow roll parameters. Second, as discussed in Sec. VC the anomaly 
effective action gives rise to dynamical dark energy, and the possibility of a departure from global FLRW 
kinematics usually assumed in cosmological models. 



B. Conformal Invariance and the CMB 



Let us consider first the consequences of conformal invariance on flat FLRW sections, independently 
of any detailed model. Our studies of fluctuations in de Sitter space suggest that the fluctuations 
responsible for the screening of A take place at the horizon scale. In that case then the microwave 
photons in the CMB reaching us from their surface of last scattering should retain some imprint of the 
effects of these fluctuations. It then becomes natural to extend the classical notion of scale invariant 
cosmological perturbations, pioneered by Harrison and Zel'dovich |135j to full conformal invariance. 
In that case the classical spectral index of the perturbations should receive corrections due to the 
anomalous scaling dimensions at the conformal phase |136j . In addition to the spectrum, the statistics 
of the CMB should reflect the non-Gaussian correlations characteristic of conformal invariance. 

Scale invariance was introduced into physics in early attempts to describe the apparently universal 
behavior observed in turbulence and second order phase transitions, which are independent of the 
particular short distance dynamical details of the system. The gradual refinement and development of 
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this simple idea of universality led to the modern theory of critical phenomena, one of whose hallmarks 
is well-defined logarithmic deviations from naive scaling relations |120| . A second general feature of 
the theory is the specification of higher point correlation functions of fluctuations according to the 
requirements of conformal invariance at the critical point [137J . 

In the language of critical phenomena, the observation of Harrison and Zel'dovich that the primordial 
density fluctuations should be characterized by a spectral index n = 1 is equivalent to the statement 
that the observable giving rise to these fluctuations has engineering or naive scaling dimension p = 2. 
This is because the density fluctuations bp are related to the metric fluctuations by Einstein's equations, 
5R ~ G5p, which is second order in derivatives of the metric. Hence, the two-point spatial correlations 
(5p(x)5p(y)} ~ (5R(x)5R(y)} should behave like \x — y\~ 4 , or \k\ x in Fourier space, according to simple 
dimensional analysis. 

One of the principal lessons of the modern theory of critical phenomena is that the transformation 
properties of observables under conformal transformations at the fixed point is not given by naive 
dimesnional analysis. Rather one should expect to find well-defined logarithmic deviations from naive 
scaling, corresponding to a (generally non-integer) dimension A. The deviation from naive scaling is 
the "anomalous" dimension of the observable due to critical fluctuations, which may be quantum or 
statistical in origin. Once A is fixed for a given observable, the requirement of conformal invariance 
determines the form of its two- and three-point correlation functions up to an arbitrary amplitude, 
without reliance on any particular dynamical model. 

Consider first the two-point function of any observable O a with dimension A. Conformal invariance 
requires [HQl [[37] 

(O a ( Xi )O a (x 2 )) ~ \xi - x 2 \~ 2A (6.15) 

at equal times in three dimensional flat spatial coordinates. In Fourier space this gives 

G 2 {k) = (6 A (k)G A {-k)) ~ |£f A - 3 . (6.16) 

Thus, we define the spectral index of this observable by 

n = 2A-3. (6.17) 

In the case that the observable is the primordial density fluctuation Sp, and in the classical limit where 
its anomalous dimension vanishes, A — > 2, we recover the Harrison- Zel'dovich spectral index of n = 1. 

In order to convert the power spectrum of primordial density fluctuations to the spectrum of fluc- 
tuations in the CMB at large angular separations we follow the standard treatment [138J relating the 
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temperature deviation to the Newtonian gravitational potential <p at the last scattering surface, ^ ~ 5ip, 
which is related to the density perturbation in turn by 

V 2 5<p = AirG5p . (6.18) 

Hence, in Fourier space, 

r~^F7' (fU9) 

and the two-point function of CMB temperature fluctuations is determined by the conformal dimension 
A to be 

<h{0) = (^(h) 6 -^(h) 

d'k (j^j 2 G 2 (k)e ik -™ ~ T(2 - A)(r 2 2 ) 2 ~ A , (6.20) 

where r\i = (f\ — r 2 )r is the vector difference between the two positions from which the CMB photons 
originate. They are at equal distance r from the observer by the assumption that the photons were 
emitted at the last scattering surface at equal cosmic time. Since r\ 2 = 2(1 — cos#)r 2 , we find then 

C 2 (0)~r(2-A)(l-cos6O 2 - A (6.21) 

for arbitrary scaling dimension A. 

Expanding the function C2(0) in multipole moments, 



C 2 (6) = ^]T(2^+l)4 2) (A)P £ (cosfl) , (6.22) 



cf (A) ~ r(2 - A) sin [vr(2 - A)] ^ + ^ + , (6.23) 
shows that the pole singularity at A = 2 appears only in the £ = monopole moment. This singularity 



is just the reflection of the fact that the Laplacian in (6.18) cannot be inverted on constant functions, 
which should be excluded. Since the CMB anisotropy is defined by removing the isotropic monopole 
moment (as well as the dipole moment), the 1 = term does not appear in the sum, and the higher 
moments of the anisotropic two-point correlation function are well-defined for A near 2. Normalizing 

(2) 

to the quadrupole moment c 2 (A), we find 

(A) - (A) T{6 ~ A) r ^~ 2 + A ) (6 24) 

93 



f2l (2"! 

which is a standard result |138j . Indeed, if A is replaced by p = 2 we obtain £(£ + l)c\ >(p) = Qc\ >(p), 
which is the well-known predicted behavior of the lower moments (I < 30) of the CMB anisotropy where 
the Sachs- Wolfe effect should dominate. 

Turning now from the two-point function of CMB fluctuations to higher point correlators, we find a 
second characteristic prediction of conformal invariance, namely non-Gaussian statistics for the CMB. 
The first correlator sensitive to this departure from gaussian statistics is the three-point function of the 
observable 0a, which takes the form |137j 



{Oa(x 1 )Oa(x 2 )Oa{x3)} ~ \xi - x 2 \ A \x 2 -x 3 \ A \x 3 -x 1 \ A , 



(6.25) 



or in Fourier space, 2 



d 6 p |p| A ~ 3 \p + h l \*-' i \p-k 2 



A-3 



|A-3 



r(3-M) 



du I dv 



[u(l-u)v] 2 (1-v) 



-1+4 



o Jo [u{l-u){l-v)k 2 1 +v{l-u)kl + uv{k 1 + k 2 ) 2 \ 



3.26) 



This three-point function of primordial density fluctuations gives rise to three-point correlations in the 



CMB by reasoning precisely analogous as that leading from Eqns. (6.16) to (6.20). That is, 

JT JT 
— (ri) — (r 2 ) — (r 3 ) 



C^{0 12 1 823, 031) 

d 3 h d 3 k 2 



k\ k\ (k x + k 2 ) 2 



G 3 (ki, k 2 ) e 



(6.27) 



where = (f« — fj)r and rf- = 2(1 — cos#j 



In the general case of three different angles, this expression for the non-Gaussian three-point corre- 

c form analogous 
9, it follows from 



lation function (6.27) is quite complicated, although it can be rewritten in parametric form analogous 
to (6.26) to facilitate numerical evaluation. In the special case of equal angles Oij 



its global scaling behavior that the three-point correlator is 

C 3 {9) ~ (1 -cosfl)i (2 - A) . 



(6.28) 



Expanding the function C 3 {6) in multiple moments as in (6.22) with coefficients c\ , and normalizing 
to the quadrupole moment, we find 



C^(A) - c (3) f M r ( 4 +j( 2 - A )) r ^-K 2 - A )) 

1 1 ' 2 { j r(2-|(2-A))r(£ + 2 + |(2-A)) 



(6.29) 



Note that (6.261 corrects two minor typographical errors in eq. (16) of Ref. [1361 
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In the limit A — > 2, we obtain 1(1 + l)cl = Qc^\ which is the same result as for the moments dp 
of the two-point correlator but with a different quadrupole amplitude. The value of this quadrupole 

(3) 

normalization c 2 (A) cannot be determined by conformal symmetry considerations alone, and requires 
more detailed dynamical information about the origin of conformal invariance in the spectrum. 

For higher point correlations, conformal invariance does not determine the total angular dependence. 
Already the four-point function takes the form, 

(0 A (x 1 )OA(x 2 )0 A (x 3 )0 A {x i )} ^— , (6.30) 

FT . . r ■ ■ 

i ii<j ij 

where the amplitude is an arbitrary function of the two cross-ratios, ^i3?*24/ r i2 r 34 an d r i4 r 23/ r i2 r 34- 
Analogous expressions hold for higher p-point functions. However in the equilateral case 9™ = 6, the 
coefficient amplitudes A p become constants and the angular dependence is again completely determined, 
with the result, 

C p {0) ~ (l-cos#)i (2 - A) . (6.31) 
The expansion in multiple moments yields coefficients cf 1 of the same form as in Eqn. (6.29) with 3/2 



replaced by p/2. In the limit A = 2, we obtain the universal ^-dependence lit + V)cf = Qcf ■ 

Again it bears emphasizing that these results depend upon the hypothesis of conformal invariance 
on the flat spatial sections of FLRW geometries, but otherwise makes no dynamical assumptions, such 
as in scalar field inflaton models. Because of the possibilities of large backreaction on the horizon in 
de Sitter space, this assumption of global spatial homogeneity and isotropy implied by the FLRW line 



element (3.20) may not in fact be realized. We study these horizon backreaction effects to linear order 



in the fluctuations next. 



C. Linear Response in de Sitter Space 



The first step to a cosmological model based on the effective action (4.89 ) going beyond the essentially 
purely kinematic considerations of the previous subsection is to study the dynamical effects of (4.89) 
in de Sitter space. This one can do by performing a linear response analysis of perturbations about de 
Sitter space with the anomaly scalars ip and i[>, in place of any ad hoc inflaton field. 



Because of the 0(4, 1) maximal symmetry, with 10 Killing generators (6.2), the maximally symmetric 
Bunch-Davies (BD) state in de Sitter spacetime is the natural one about which to consider perturbations. 
This is the state usually considered in cosmological perturbation theory. The linear response approach 
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requires a self-consistent solution of equations (2.42), around which we perturb the metric and stress 
tensor together. In the BD state 0(4, 1) de Sitter symmetry implies that (TV) is also proportional to 
5 a b , which therefore guarantees that de Sitter space is a self-consistent solution of the semi-classical eqs. 



(2.42). The self-consistent value of the scalar curvature R, including the quantum contribution from 



(TV) is determined by the trace of (2.42), i.e. 



-R + Ak = 8irGb'E 



4vrG6' 



R 
12 



A 



1 



R 2 

3 

16vrG6'A 



or 



+ ... 



(6.32a) 
(6.32b) 



3 V 3 

in an expansion around the classical de Sitter solution with GA|6'| <C 1. Thus, in this limit the stress 



tensor source of the semi-classical Einstein's equations (2.42) in the BD state (5.13) gives a small finite 



correction to the classical cosmological term in the self-consistent de Sitter solution. The trace of (5.13) 
is exactly b'E = b'R 2 /6 used in determining the self-consistent scalar curvature including the quantum 



BD corrections in (6.32). The solutions for the anomaly scalar fields, 



tp = 21no = 2Ht, 
^ = 0, 



(6.33a) 
(6.33b) 



correspond exactly to the BD state, and it is consistent to expand the semi-classical Einstein eqs. (2.42) 



around de Sitter space using the effective action and stress tensor of the anomaly scalar fields with 



background values (6.33) 



With the self-consistent BD de Sitter solution (5.13), (6.32), and (6.33), one may consider the linear 



response variation of the semi-classical Einstein eqs. (2.42) [ 139} 1140] 

s{R a b -^8 a b + A6%)=8irG6(T%) R , 



(6.34) 



with the source the anomaly generated stress tensor (5.1). Considering the general form of the exact 



quantum effective action (4.86), its variation includes three kinds of terms, corresponding to the local 



terms up to fourth order in derivatives of the metric, Weyl invariant non-local terms, and those terms 



coming from the anomaly generated effective action ( 4.63 )-( 4.64). Parameterizing the local fourth order 



geometric terms in (4.83) with finite coefficients a R and j3 R , (6.34) with (4.86) gives then 

SR 



SRV 



5\ = 8ttG 5 



(T a b ) loc + (T a b ) mv + (T a b ) 



8ttG [-a R 5A a b - (3 R 5B\ + (T%) inv + b' 5E a b + b 8F% 



6 ' 



(6.35) 
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where all terms to linear order in 5g a b = h a b and in the variations of the auxiliary fields, 



5(p — <p — (p = 4> 
5ip = ip — if) = if> . 



(6.36a) 
(6.36b) 



are to be retained. All indices are raised and lowered by the background de Sitter metric g a b at linear 
order in the perturbations. Since (T a b ) mv is the variation of a Weyl invariant action, it has zero trace. 

Since the terms coming from the fourth order local invariants C a b c dC abcd and R 2 , namely SA a b and 
SB a b are higher order in derivatives than the Einstein-Hilbert terms, they are important only in the 
extreme ultraviolet regime at energies of order Mpi, where in any case one should not trust the semi- 
classical effective theory. These local terms were explicitly analyzed in Ref. |140| and do not affect any 



physics on the cosmological horizon scale H <C Mpi, In particular, from the trace of (6.35) one obtains 
the purely local eq., 



3/3 R + - □ + 



1 



+ 2b' H 2 



5R = 0. 



3.37) 



3 J " ' \16irG 

Because of the G~ l = M pl term, the only non-trivial solutions of (6.37) are on Plank scale, and outside 
the range of applicability of the EFT approach. Hence we restrict attention to only the remaining 



solutions of (6.35) and (6.37) that satisfy 



SR = 0. 



(6.38) 



Notice that by so doing we are excluding the trace anomaly driven inflationary solutions studied by 
Starobinsky |112j . 

In general the non-local Weyl invariant terms in the exact quantum effective action S{ nv are difficult 
to calculate and are not known. However in the case of conformally invariant matter /radiation fields 
in the conformaly flat de Sitter geometry, it is possible to keep track of these non-local terms, and 
calculate the linearized stress tensor they provide 5(T a b ) mv completely from the corresponding quantity 
in the conformally related flat space |141| 1142] . By keeping track of these terms in the de Sitter linear 
response analysis, we will then be able to provide a non-trivial check on our general argument based on 



the classification of terms in (4.86) according to their response under global Weyl transformations that 



these terms play no role in the infrared, and the Wilson effective action for low energy gravity consists 



of the classical Einstein-Hilbert and anomaly terms only in (4.89) 



The last two terms in the last line of (6.35) arise from the anomaly action and give rise to new 
physical effects in cosmology. The variation of the stress tensors E a b and F a b depend on both the 
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variations of the metric and the variations of the auxiliary fields. The variation of the auxiliary field 



equations (4.62) gives 



R 



6(n 2 <p)-j6(nip) 
«y(nV)-^*(n^) 





- □ + 2H 2 ) 5{- □ <p) = -2(V a V b £)5R a b , 
-□ +2H 2 ) (-n)tp = 0. 



(6.39a) 
(6.39b) 



The first of these two equations shows that at linear order there is mixing between the fluctuations of 
the anomaly scalar field ip and the metric perturbation 5g ab = h ab around de Sitter space. This mixing 
is algebraically simplest to study by making a suitable gauge choice, although of course the results in 
the end must be gauge invariant. Here we present only the results of the gauge invariant analysis |140j . 

The metric perturbations which are scalar with respect to the background three-metric gij = a 2 r]ij 
can be parameterized in terms of four functions, (A,B,C,£) in the form [1431 1144] , 



htt 
h tj 
hij 



-2A 



adjB — > iakjB 



2a z 



V 3 



(6.40a) 
(6.40b) 
(6.40c) 



in momentum space. As is well known, only two linear combinations of these four functions are gauge 
invariant (for k ^ 0). The two gauge invariant metric perturbation variables may be taken to be 



T A = A + d T {aB) - d T {a 2 d T £) , 
V 2 



T C ^C 



£ + (d T a)B — a(d T a)d T £ , 



(6.41a) 
(6.41b) 



which correspond to the gauge invariant Bardeen-Stewart potentials denoted by $,4 and <&c m Refs. 
[143, 1144] . For k = there is only one gauge invariant metric combination, namely d T T c — HT A . 
The variation for 5R can be written in the form 



5R 



-ah + v a v b h ab -R ao h ab 



jab 1 



6(d 2 T c - Hd T T A ) + 24H(d T T c - HT A 



2 9 4 n 

V T V T 

a 2 A a 2 c ' 



(6.42) 



where V 2 is the Laplacian on the flat M 3 FLRW sections. Hence condition (6.38), 5R = 0, is gauge 
invariant, and equivalent to 



8 



HI. 



V 2 



T, 



2V Z 



X. =0. 



(6.43) 



dr ' ""J \8t "~ a J 3a 2 A 3a 2 ^ c 
which provides one constraint between the two gauge invariant potentials T, and T c . This means that 
there remains only one gauge invariant metric function to be determined by linear response in this scalar 



!i> 



sector. As proven in [140J the information about this remaining metric degree of freedom is contained 



completely in the rr component of the linear response equations (6.35). Thus, recalling the condition 



(6.38), one can define the gauge invariant quantity 

2a 2 



H 2 



5G 7 



2a 2 
H 2 



5R 7 



12a 2 fdT 



H 



dr 



HT. 



H 2 



V 2 T C , 



(6.44) 



which appears in the rr component of the linear response eq. (6.35), and which contains the only 



remaining gauge invariant information in the sector of scalar metric perturbations (6.40) after condition 



(6.38) has been imposed. 



It is easily checked and verified in Ref. [140] that the quantity 

* = 4> + 2{d T a)B - 2a(d T a){d T £) 



(6.45) 



is gauge invariant. This is similar to the gauge invariant variable that can be constructed from the 
scalar field in scalar infiaton models of slow roll inflation [132] . The second anomaly scalar field tp is 
already gauge invariant. Defining also the explicitly gauge invariant quantities, 



D $ + 



2 dT A 



H dr H dr 



and 



1 



d V 5 



1>: 



where the differential operator T> n is defined for arbitrary integer n by 

1 / r) 2 r) X7 2 

' + {2n + l)H^- + n{n + l)H 2 



h 2 \d 



dr 



(6.46) 



(6.47) 



(6.48) 



one finds that the linear response equations around de Sitter space can be written in terms of the gauge 
invariant variables u, v, 5R, and q. With the condition 5R = imposed, the equations for u and v are 



found from (6.39) and become simply 



V 2 u = 


H 2 


( d 2 

\dr 2 




+ 6H 2 - 


V 2 

a 2 


V 2 v = 


IP 


[dr 2 




+ 6H 2 - 


V 2 
a 2 



(6.49a) 
(6.49b) 



while that for q is the somewhat more complicated, because of the non-local Weyl invariant terms in 
5(T a b ) mv . For the case of a conformal matter /radiation field in de Sitter space these terms can be 
computed exactly j!41l 1142] , with the result |14QJ 

e k 2 



be\ k 2 
1 - R - — J q = eHr (V iq ) - a R (Dig) + e' u - -j^v 



e 
2^ 



dr\ K(rj -rf;k; //) [a 2 (Dig)],,/ 



(6.50) 



'/o 
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in Fourier space where —V 2 — > k 2 , and [a 2 (T>iq)] in the integrand is evaluated at rj . The u and v terms 



coming from the stress tensor of the anomaly and obeying the eqs. (6.49) are particular realizations of 



the possible state dependent terms mentioned in Ref. |103| . The last term in (6.50) is the non-local 
Weyl invariant term from 5(T a b y nv that involves the kernel K is defined in conformal time r/ by 



K{r]-rf;k;ij) = 



2^ 



-uj 2 + k 2 — ie sgnu> 
J 2 ^ 



(6.51) 



It depends on an arbitrary renormlization scale [i, and derives from the Weyl invariant term Si nv in 
the exact quantum one- loop action around de Sitter space. The small (order K) a R ,e,e' parameters are 
defined by 



Ctr 



p. 



R — 
£ = 

e' = 



16irGH 2 a R , 
l927tGH 2 f3 R 
32irGH 2 b , 



32vr 



GH 2 b' . 



(6.52a) 
(6.52b) 
(6.52c) 
(6.52d) 



Let us make some observations about the exact linear response eqs. ( 6.49 )-( 6.50) for conformal quantum 



fields about de Sitter space. First it is important to point out that although the anomaly scalar eqs. 



(4.62) are fourth order in time derivatives, in fact only the gauge invariant quantities u and v satisfying 



the second order eqs. (6.49) couple to the linearized Einstein eq. (6.50) to linear order. In other words, 



half of the solutions of the fourth order eqs. (6.35) are annihilated by the operator T>q in ( 6.46 )-( 6.47) 



and decouple entirely. Only the other half of the solutions with non-zero (u,v) satisfying (6.49) couple 



to the physical metric perturbations. This seems to be the analog of the elimination of the ghost degrees 



of freedom of the fourth order action T\yz of (4.59) by the diffeomorphism constraints found previously 
in its quantization on IR (g> §> 3 [130 . 

Secondly, although the variations of the fourth order invaraints in Si oca i are fourth order in time 



derivatives, they also appear in (6.50) only with the second order differential operator (6.48) on the 



gauge invariant quantity q. This is a consequence of the fact that the classical Einstein eqs. contain no 



dynamical metric perturbations in the scalar sector (6.40) at all, and would be completely constrained 



(i.e. q = 0) if it were not for the quantum corrections on the right side of (6.50) which vanish in the 
limit h — > 0. When h ^ 0, the character of the eqs. changes discontinuously in that there are new 



non-trivial gauge invariant solutions of the second order equations ( 6.49 )-( 6.50) 



Next we showed in Ref. [140J that most of the complications of the exact linear response eq. (6.50) 



are irrelevant for macroscopic physics. First, because all the dimensionless parameters defined in (6.52) 
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are very small compared to unity for H « Mpi, the corrections in the parentheses on the left side of 



(6.50) may be neglected. Then we analyzed the general homogeneous solutions of (6.49) and (6.50), i.e. 
with u = v = but including the effect of the non-local term involving K. If u = v = the non-trivial 



homogeneous solutions of (6.50) involve oscillations on the Planck scale Mpi, and hence lie outside the 
range of applicability of the effective action of low energy gravity. Hence all these homogeneous solutions 



of (6.50) may be consistently neglected in the low energy theory for k <C Mpi, just as the non-trivial 



solutions of (6.37) with 5R ^ lie outside the range of applicability of the semi-classical effective theory 



and should be discarded as unreliable. With the homogeneous solutions of (6.50) thereby discarded for 
the same reason, the only non-trivial solutions remaining are the inhomogeneous ones for which either 



u 7^ or v 7^ (or both). For these solutions all the terms on the right side of (6.50) linear in q are 
suppressed by powers of H/Mpi <C 1 or kp^yg/Mpi = k/aMpi -C 1 compared to the u, v terms. 

Thus, as far as the predictions valid for macroscopic gravity are concerned with physical wavelengths 
of perturbations much greater than the Planck length, one may replace the complicated, non-local exact 



(6.50) by the much simpler and fully determined 
5R\ 



2a 2 



16ttGH 2 V 2 . L , 

T [b u + bv) 



(6.53) 



so that 5R T T if and only if driven by the non-trivial solutions of the additional degrees of freedom 
provided by the anomaly scalar field equations ( 6.49| ). Hence we have verified explicitly by this analysis 
of the full linear response of conformal matter /radiation field fluctuations around de Sitter space that the 
^lolai an d Sinv terms in the exact quantum effective action do not influence low energy or macroscopic 



physics, consistent with the general classification and Weyl scaling arguments of Sec. IV Instead 



we could have started with only the low energy effective action (4.89) and obtained (6.53) much more 
directly, which is correct at scales k p h ys <C Mpi . Conversely, without inclusion of the anomaly generated 
terms one would miss entirely the physics associated with the degrees of freedom u and v and the 



solutions (6.49) in the scalar sector of cosmological perturbation theory not present in the purely classical 
theory for which 5R T T = in this scalar sector. 



The general solution of (6.49) for either u or v in FLWR coordinates is easily found in Fourier space, 
namely 



1 



u± = v± 



exp ± 



ik 
~H~a 



Ak-x 



2 2 ^fikr)+ik-x 



H z r] z e 



(6.54) 



Thus, the auxiliary fields of the anomaly action yield the non-trivial gauge invariant solutions (6.53) 



with (6.54) for the linearized Ricci tensor perturbations 5R 7 



-5R 1 : in the scalar sector, which clearly 
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are not present in the purely classical Einstein theory without the anomaly action. Since with (6.38) 

2V 2 



<)R\ <)R', -Gil ( - HT A ] + =^-T L ■ 



(6.55) 



using the solution for 5R T T (6.53) in (6.55), differentiating and substituting in (6.43), one can solve for 



Y c , to obtain 



-T c = 8itGH 2 (h^- + 4H 2 + ^) {b'u + bv) , 



(6.56) 



where we have canceled a common factor of V 2 from both sides (valid for k ^ 0). This can be rewritten 

as 

72 • 



V 2 T C = 8irGH 2 (h^- + 2H 2 + ^ [a 2 (b'u + bo)] . 



Since from (6.54) 



a 2 (b'u + bv) = [c + e~ ikr) + c_e 
with c± constants (depending on b, b'), at late times rj — > 0, 



V 2 T C -> 167rGff 4 (c+ + c_)e^ as r -> oo . 



(6.57) 



(6.58) 



(6.59) 



Thus the Bardeen-Stewart potential T c describing the gauge invariant linearized perturbations of the 
de Sitter geometry, generated by the stress tensor of the conformal anomaly in the scalar sector remains 
non- vanishing at late times for every k (while and 5R T T falls off with a~ 2 ). 



Being solutions of (6.49) which itself is independent of the Planck scale, these new solutions due to 



the effective action of the anomaly vary instead on arbitrary scales determined by the wavevector k, and 
are therefore genuine modes of the semi-classical effective theory. This is a non-trivial result since these 
modes appear in the tracefree sector of the semi-classical Einstein equations, with 5R = 0, and hence 
cannot be deduced directly from the local form of the trace anomaly itself, but only with the help of the 



covariant action functional (4.61) and the additional scalar degrees of freedom which the local form of 



this action implies. These additional modes, which couple to the scalar sector of metric perturbations 
in a gauge invariant way are due to a quantum effect because the auxiliary scalar fields from which they 
arise are part of the one-loop effective action for conformally invariant quantum fields, rather than a 
classical action for an ad hoc scalar inflaton field usually considered in inflationary models [132t 1133] . 
This demonstrates the relevance of the anomaly action for describing physical infrared fluctuations in 
the effective semi-classical theory of gravity, on macroscopic or cosmological scales unrelated to the 
Planck scale. 
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The Newtonian gravitational constant G and the Planck scale enter Eq. ( 6.53 ) only through the small 



coupling parameters e and e' between the anomaly scalar fields and the perturbation of the geometry 



5R T T , which are defined in (6.52). Thus in the limit of either flat space, or arbitrarily weak coupling 
GH 2 — > the modes due to the anomaly scalars decouple from the metric perturbations at linear order. 



As in the case of the gravitational scattering of Sec. IV this explains why the anomaly scalars are so 



weakly coupled and difficult to detect directly in the flat space limit. 



The effective action (4.89) expanded to quadratic order about the self-consistent de Sitter solution 



(6.32) in terms of the gauge invariant variables u, v, Ta, and Tc is 

H A u 2 H 2 u5R 



S 



(2) I 
effldS 



>7 



Sg + b' I d 3 x dr a 



+ 



+b I d 3 x dr a 3 



-H A uv + 



H 2 vSR 4 In a 



+ 



3 a 4 



V 2 (T, 



T 



c 



(6.60) 



where 



S G 



8ttG 



d 3 x dr a 3 



dr 



c 



Or 



+ 6HT A 



8T 



c 



+\ (VT A ) ■ (VT C ) + 



(VT C ) 2 



dr 
3H 2 T 2 A 



(6.61) 



is the Einstein-Hilbert part of the action, and 5R is given by (6.38). Varying (6.60) with respect to 



$,ip, Ta and Tc, and setting 5R = yields the low energy form of the gauge invariant linear response 



equations (6.49) and (6.53). The 5R terms in (6.60) cannot be set to zero until after the variations 



of S^g are performed since they generate the coupling between the u,v and metric variables T^,T 



c- 



Note that the variations of (6.60) must be performed with respect to the original set of gauge invariant 



variables <!>, ip, Ta, Tc, in order to obtain eqs. (6.46)-(6.49) and (6.53). The linearized equations for 



the anomaly scalar fields so obtained are equivalent to (6.39), when it is recognized that 



A 4 ^ = (- □ + 2H 2 )(- D)V> = H A V 2 (V i(>) , 



(6.62) 



i.e. that the fourth order conformal differential operator A4 associated with the Euler- Gauss-Bonnet 
invariant of the trace anomaly factorizes into two second order differential operators in two different 
ways in de Sitter space. As a practical matter this means that any solution of the fourth order eq. 



(5.19), for which T>^ 7^ 0, automatically is a solution of (6.49) for u or v. The disappearance from the 



action (6.60) of any solution of the scalar field eqs. (6.35) that is in the kernel of T>q is again noteworthy. 
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D. Cosmological Horizon Modes 



The solutions ( 6.53 )-( 6.54) corresponds to a linearized stress tensor perturbation of the form 



Thus each mode of fixed k redshifts like mode of a classical conformal radiation field. 

Not apparent from this form of this stress tensor perturbation for a given k mode in FLRW co- 
ordinates is their possible relevance to the vacuum polarization effects near the horizon discussed in 
previous sections. To address this we need the other components of the stress tensor perturbation, and 



the behavior in the static coordinates of de Sitter space (3.27). 

The other components of the stress tensor perturbation (6.63) in FLRW coordinates can be found by 
a general tensor decomposition for scalar perturbations analogous to Eqs. (6.40) for the metric. That 
is, the general perturbation of the stress tensor 5{T ab ) R in the scalar sector can be expressed in terms 
of 5{T TT ) plus three additional functions. These three functions are determined by the conditions of 
covariant conservation, 



V b 6(T ab ) R = (6.64) 
for a = t and a = i (two conditions), plus the tracefree condition, 

6(T a a ) R = Q, (6.65) 



imposed as a result of the 5R = condition. Choosing the arbitrary proportionality constant in (6.63) 
to be unity, a straightforward calculation using these conditions and the Christoffel coefficients in the 
flat FLRW coordinates gives |140j 

S(T TT ) R =H 2 ^ e Tikri+ik-x = y2 ^ik^+ik-S > 

6{T Ti ) R = ±H 2 ^ e ^kv+ik-x = IP ^ikn+ik-S ^ (a66b) 

Ob Ob CJX CJT 
5{T ij ) R = H 2 — ^- e Tikv+ik-x = eT ik v +ik-x (6 _g 6c) 

a° a 3 ox t ox3 

for the other components of the stress tensor variation for these modes in the flat FLRW coordinates 
of de Sitter space. If one averages this form over the spatial direction of k, a spatially homogeneous, 
isotropic stress tensor is obtained with pressure p = p/3. In FLRW coordinates this averaging describes 
incoherent or mixed state thermal perturbations of the stress tensor which are just those of massless 
radiation. 
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Next we use (3.21) and (3.26) to read off the coordinate transformation from FLRW flat coordinates 



r and g = \x\ to static t and r coordinates, given by 



r = a \x\ = ag = g e 
1 
2H 



t = T 



Ht 

2 2 2Ht\ 



In (1 - H 2 g 2 e 2Hr ) . 



(6.67a) 
(6.67b) 



The inverse transformations are 



y g Ht 

g = \x\ = — , 

' ' VI - H 2 r 2 

r = t + ~Hl-H 2 r 2 ). 

The Jacobian matrix of this 2x2 transformation of bases is 

l 



(6, 
(6.68b) 



d(t,r) 



( dt 


dt \ 




dr 


dg 


H 


dr 


dr 




\ dr 


dg ) 





Hr 



1-H 2 r 2 g(l-H 2 r 2 ) 

Hr r/g 



(6.69) 



<9(t, q) 

Using these relations, one may express the action of the differential operators in Eq. (6.49) in terms of 



the static coordinates (3.27) instead. Since 



1 d 



g 2 dg \ dg 







+ 



L 



21 



(6.70) 



where — L is the scalar Laplacian on the sphere S , a straightforward calculation using (6.69) shows 



that 



V 2 

2 V 



(1-H 2 r 2 ) 2 



H 2 r 2 ^+H(3-H 2 r 2 : 



dt 2 



d_ 
dt 



2Hr d 2 v 1 d 



1 — H 2 r 2 dtdr r 2 dr 



,dv 
Or 



+ 



L 2 



(6.71) 



while 



H 2 v 



dr 2 dr 



v- 



tp = H 2 V ip 



1 



ay 

1 - H 2 r 2 \ ~W 



2H 



dip 
~di 



'1 - H z r 



2J2\ 



1 3 / 2 <9'0 \ L 



(6.72) 



and 



V 2 v = 



2 d 

+ 5H— + 6H 2 
dr z dr 



a 2 



(d\ 



r dv\ 



d\ 



1^2-2 + 2H Yt ) ~ (1 " fV) 9^ 



+2 3H 2 r 



l\ dv 
dr 



L 2 

+ 6H 2 v + -^v, 



(6.73) 
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thus converting these differential operators from flat FLRW to de Sitter static coordinates (3.27). 



In (5.16) and (5.20) the general solutions of the homogeneous eq. (5.19) as functions of the static r 



in de Sitter space are given. Of the four, 



1 /1-Hr 

In 

Hr \1+Hr 



In 



1 - Hr 



1 



(6.74) 



1 + Hr 

the second and third solutions are solutions of the second order eq. for a minimally coupled scalar field, 
i.e. they satisfy □'0 = 0, while the first and last solution satisfy the second order eq, for a conformally 
coupled scalar field, i.e. they satisfy (— □ + 2H 2 )t/j = 0. The fourth solution 1/r is singular at the 
origin and so was not considered in Sec. |Vj In any case this solution and the constant solution in (I6.74I) 



give vanishing contribution to v in (6.72) while the first and second solutions give for v, 

4 4 l + H 2 r 2 



(6.75) 



1 - H 2 r 2 ' Hr 1 - H 2 r 2 ' 
respectively. The second gives a singular contribution to v and the stress tensor at r = 0, so we consider 



only the first solution. Note that if we allow also for a term linear in t in the list (6.74) of homogeneous 



solutions to (5.19), then from (6.72) it would also produce the identical form for v as the first member 



of (6.75). Thus the general spherically symmetric static solution to the homogeneous eq. (5.19), which 



because of (6.62) is also a solution of (6.73), which is non-singular at the origin is 

2 1 

-{b'u + bv) oc 1 _ ff2r2 . 



(6.76) 



Choosing the arbitrary normalization constant of proportionality to be unity, then from (6.53), (6.71) 



and (6.76), we obtain 



5{T 7 



-H 2 



1 

1 - H 2 r 2 



H A /3 2H 2 r 2 
'I - H 2 r 2 ) 2 { 2 + 1 - H 2 r 2 



(6.77) 



This shows that a linear superposition of solutions (6.54) of the linear response equations in static 



coordinates can lead to gauge invariant perturbations which diverge on the de Sitter horizon. 



To see what stress tensor (6.77) corresponds to in the static (t, r) coordinates, we use the form of the 

n FL 

dt\ 



other components in (6.66) in FLRW coordinates, (6.69), and the transformation relation for tensors, 

2 



dt\ ( dt 



( dt 



dt \ 



r T 3 



(6.78) 



with (6.69), (6.71), to obtain 



-(l-H 2 r 2 )5(T a ) R 



QH 4 
1 - H 2 r 2 ) 2 ' 



(6.79) 
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Here use has also been made of the identities, 



dt dt dg dt Xi 
dx % dg dx % dg g 



(6.80) 



and 



d_ 

Xid 1 



dr d dt d 

1 = Hr — 

dr dr dr dt dr 



x 



d_ 
' dg 



d_ 

ir 
dr d 



+ 



dt d 



dg dr dg dt 
d 



dr 



XiXjd l d J = (x l di)(x J dj) — x l di = r 



d 



dr V dr 

valid when operating on functions of r only. Likewise we find 

2H 4 



dr 



(6.81a) 
(6.81b) 

(6.81c) 



5(T r r ) R = 5{T« ( 



5{T\ 



1 



(6.82) 



(l-# 2 r 2 ) 2 3 

corresponding to a perfect fluid with p = p/3, but now in static coordinates. As we know from the 



conservation eq. (2.56) in static coordinates, and (3.43), this eq. of state leads to a quadratic divergence 
in the metric factor (1 — H 2 r 2 )~ 2 as r — > r H . This is the exact analog of the large blueshift factor 



discussed previously for the Schwarzschild case, (2.41). Because of the buildup of the effects of the 
anomaly scalar fields on the cosmological horizon in de Sitter space, they may be called cosmological 
horizon modes. 



The form of the stress tensor (6.79), (6.82) is the form of a finite temperature fluctuation away from 
the Hawking-de Sitter temperature T H = H/2ir of the Bunch-Davies state in static coordinates 



Since the equation the solutions (6.74) satisfy is the same as (5.19), it follows that there exist linear 



combinations of the solutions (6.54) found in Sec. VIC which give (6.76) and the diverging behavior 
of the linearized stress tensor on the horizon, corresponding to this global fluctuation in temperature 
over the volume enclosed by the de Sitter horizon. Since the solution i/j = pt also generates a solution 



of (6.49) for v of the same form as the first member of (6.75), we can expect that perturbations of the 
quantum state described by a non-zero time derivative (in static time coordinates) within one causal 



horizon volume at some initial time t = will give rise to very large stress tensors of the form of (6.79) 



and (6.82) on tha horizon of that causal region at later times. 

Note that in static coordinates the stress tensor p = p/3 does not involve averaging over directions 



of k, but a particular coherent linear superposition over modes (6.54) with different k in order to obtain 



a particular isotropic but spatially inhomogeneous solution of (6.76). This selects a preferred origin and 



corresponding horizon in static coordinates (3.27). The fluctuations in Hawking-de Sitter temperature 
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thus preserve an 0(3) subgroup of the de Sitter isometry group 0(4, 1). Clearly the origin chosen at 
r = is completely arbitrary, and a particular 0(3) subgroup is chosen at random by the perturbation. 
This is similar to a bubble nucleation process in a homogeneous thermal ensemble: no origin for the 
bubble nucleation center is preferred over any other, but some definite particular origin around which 
the bubble nucleates is selected by the fluctuations. In other words, the de Sitter group 0(4, 1) is 
spontaneously broken to 0(3) by any random fluctuation of temperature away from its Hawking-de 
Sitter value, centered on some arbitrary but definite point of de Sitter space at an initial time t = 0. 
The free energy functional relevant to this spontaneous symmetry breaking of de Sitter symmetry is 



the quadratic effective action (6.60), in which the last In a term acts as a the negative of a potential term 



(for b > and In a > 0) favoring a non-zero and spatially inhomogeneous V 2 (T 4 — T c ) perturbation 
away from de Sitter space. The —b'u 2 term behaves similarly (for b' < 0). In FLRW coordinates the 
action is time dependent. However, in static coordinates the volume measure is time independent and 
it is clear that a solution for u of the form (6.76) makes the Euclidean continued action under t — > it 



arbitrarily negative from its divergent behavior (6.76) at the horizon. Thus perturbations of this kind 



should be energetically favored, and destabilize global de Sitter space in favor of a finite patch up to 
the horizon, centered on a fixed origin. To understand what becomes of the near horizon region clearly 
requires going beyond the linearized expansion around the de Sitter background. 



Referring back to the full anomalous effective action (4.89)-(4.64), one observes that the 
C a p^C a ^ u <p term acts negative potential term in general as well, favoring the breaking of Weyl 
invariance unless cp vanishes, which is impossible in de Sitter space because of the non-vanishing source 



term for tp in (4.62). Note also that the Euler-Gauss-Bonnet source term E is non- vanishing for any 



solution of the vacuum Einstein's eqs. with A / 0. Only A e yj = (and UR = 0) can give a vanishing 
source for ip, and eliminate the C 2 (p term in the effective action. Thus the effective action of the trace 
anomaly provides a possible way to distinguish vacua with different values of the vacuum energy in 
gravity, and in the absence of any boundary terms, selecting A e ff = 0. 



To follow the diverging behavior (6.79), (6.82) all the way to the horizon one would clearly require a 



linear combination of the solutions (6.54) with large Fourier components. However once 87rG times the 
perturbed stress tensor in (6.79) becomes of the same order as the classical background Ricci tensor H 2 , 



the linear response theory breaks down and non-linear backreaction effects must be taken into account. 



The perturbation becomes of the same order as the background at r 



Ar near the horizon, where 



Ar ~ Lpi , 



(6.83) 
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or because of the line element (3.27), at the proper distance from the horizon of 

£ ~ y/r H L Pl » Lpi . (6.84) 
This correspond to a maximum k p h ys ~ l/£ <C Mpj, where the semi-classical description may still be 



trusted. At the distance (6.84) from the horizon, the state dependent contribution to the stress-energy 
tensor becomes comparable to the classical de Sitter background curvature, the linear approximation 
breaks down, and non-linear backreaction effects may be expected. 

VII. GRAVITATIONAL CONDENSATE STARS 

The paradoxes of black hole physics and the possibility of large backreaction effects on the 
Schwarzschild horizon were discussed in Sec. Ill) The problem of dark energy and corresponding possible 



large backreaction effects on the de Sitter horizon were discussed in Sees. Ill and VI. A considerable 
technical machinery and intuition associated with quantum effects near horizons has been elaborated. 
The effective action of the trace anomaly contains massless degrees of freedom which can become sig- 
nificant and potentially lead to large backreaction effects on the classical geometry close to horizons. In 
addition the fluctuations of the massless scalar degrees of freedom associated with the conformal factor 
of the metric described by the anomaly action allow the vacuum energy to become dynamical in which 
case it may vary in both space and time, unlike in the classical Einstein theory. 

To bring all these ideas and results together it is attractive to consider the possibility of a resolution 
of both the black hole and cosmological dark energy problems at one stroke, by matching an interior de 
Sitter geometry with one positive value of the vacuum energy to an exterior Schwarzschild geometry with 
another smaller value of the vacuum energy (which we can take to be zero) . The matching should take 
place in a small but finite region at the position of the common de Sitter and Schwarzschild boundary 
r n — r s- A spherically symmetric vacuum "bubble" of this kind was suggested several years ago, and 
called a gravitational condensate star or gravastar |28| [29"1 1145] . Similar ideas in which the black hole 
event horizon is replaced by a quantum critical surface were proposed in Refs. |146} 1147] . 

In the gravastar model, we replace the effective action and the stress tensor it generates by an easier 
to treat effective fluid description, and reconsider the Einstein eqs. with a perfect fluid source as a 
simplified model of the quantum effects to be described finally in the full EFT. In an effective mean 



field treatment for a perfect fluid at rest in the coordinates (2.1), any static, spherically symmetric 



collapsed object must satisfy the Einstein eqs. (2.54), together with the conservation eq. (2.56) which 



ensures that the other components of the Einstein eqs. are satisfied. In the general spherically symmetric 
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situation the tangential pressure, p± = T d g = is not necessarily equal to the radial normal pressure 
p = T r r . However, the simplest possibility which illustrates the main idea is to take p± = p, (except 
possibly at the boundaries between layers). In that case, we have three first order eqs. for four unknown 
functions of r, viz. f, h, p, and p. The system becomes closed when an eq. of state for the fluid, relating 
p and p is specified. If we define the mass function m(r) by 

2Gm(r) 



h(r) 



(7.1) 



so that (2.54a) becomes 



dm 
dr 



Airr p . 



(7.2) 



and eliminate / between (2.54b) and (2.56), we obtain 

dp G(p + p) (m + 4-7rr 3 p) 
dr 



(7.3) 



r(r — 2Gm) 

which is the TOV equation of hydrostatic equilibrium [3]. 

Because of the considerations of the previous sections, we allow for three different regions with the 
three different eqs. of state, 



I. Interior : < r < r\ , p = —p , 

II. Thin Shell : r\ < r < r<i , p = +p , 

III. Exterior : r2<r, p = p = . 



(7.4) 



In the interior region p = —p is a constant from (2.56). This is an effective cosmological "constant" in 
the interior. Let us call the constant p v = 3Hq /8ttG. If we require that the origin is free of any mass 
singularity then the interior is determined to be a region of de Sitter spacetime in static coordinates, 
i.e. 



I. f(r) = Ch{r) = C{l-H$r z ), < r < n , 



(7.5) 



where C is an arbitrary constant, corresponding to the freedom to redefine the interior time coordinate. 
The unique solution in the exterior vacuum region which approaches flat spacetime as r — > oo is a 



region of Schwarzschild spacetime (2.2), viz. 

III. f(r) = h(r) = 1 



2GM 



r 2 < r . 



(7.6) 



The integration constant M is the total mass of the object. 
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The only non- vacuum region is region II. Let us define the dimensionless variable uu in this section 
by w = 8irGr 2 p, so that eqs. ( 2.54 )-( 2.56) with p = p may be recast in the form, 

dh 



dr 
r 



1 — w — h ' 



(7.7) 



dh 1 — w — h dw 
h 1 + w — 3h w 



(7.8) 



together with pf oc wf/r 2 a constant. Eq. (7.7) is equivalent to the definition of the (rescaled) Tolman 
mass function by h = 1 — 2m(r)/r and dm(r) = 4ttG pr 2 dr = wdr/2 within the shell. Eq. (7.8) can 



be solved only numerically in general. However, it is possible to obtain an analytic solution in the thin 



shell limit, < h <C 1, for in this limit we can set h to zero on the right side of (7.8) to leading order, 
and integrate it immediately to obtain 



, , 2Gm (1 + w) 2 1 

h = 1 ~ e — < 1 , 

r w 



(7.9) 



in region II, where e is an integration constant. Because of the condition h <C 1, we require e<l, with 



w of order unity. Making use of eqs. (7.7)-(7.8) and (7.9) we have 



dr 



-e dw 



:i + w) 



■io- 



(7.10) 



Because of the approximation e <C 1, the radius r hardly changes within region II, and dr is of order 
edw. The final unknown function / is given by / = (r /r\) 2 {wi/w)f{ri) ~ [w\/w)f{ri) for small e, 
showing that / is also of order e everywhere within region II and its boundaries. 

At each of the two interfaces at r = r% and r = r2 the induced three dimensional metric must be 
continuous. Hence r and f(r) are continuous at the interfaces, and 



f(r 2 ) ~ — fin) = (1 - iVi) = 1 

u>2 w 2 r 2 



To leading order in e 1 this relation implies that 

1 



ri ~ — ~ 2GM ~ r 2 



(7.11) 



(7.12) 



Thus the interfaces describing the phase boundaries at n and r 2 are very close to the classical event 
horizons of the interior de Sitter and exterior Schwarzschild geometries, while the full solution has no 
horizon at all. 

The significance of < e <ti 1 is that both / and h are of order e in region II, but are nowhere 
vanishing. Hence there is no event horizon, and t is a global time. A photon experiences a very large, 
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0(e 2) but finite blue shift in falling into the shell from infinity. The proper thickness of the shell 
between these interface boundaries is 



1 1 f Wl 3 1 



r 2 

I drh~ 2 ~ r s ea / dww~ 2 ~e 2 r s , (7.13) 



and very small for e — > 0. Because of (7.12) the constant vacuum energy density in the interior is just 
the total mass M divided by the volume, i.e. p v ~ 3M/47rr^, to leading order in e. The energy within 
the shell itself, 

E u = 4:ir pr 2 dr ~ eM — (1 + w) ~ eM (7.14) 
J n J W2 

is extremely small. 

We can estimate the size of e and I by consideration of the expectation value of the quantum stress 
tensor in the static exterior Schwarzschild spacetime. In the static vacuum state corresponding to no 
incoming or outgoing quanta at large distances from the object, i.e. the Boulware vacuum [35], the 
stress tensor near r = r s is the negative of the stress tensor of massless radiation at the blue shifted 
temperature, T/ oc = T H /^ /(r) and diverges as ~ f 2 (r) as r — > r s . The location of the outer 
interface occurs at an r where this local stress-energy oc M~ 4 e -2 , becomes large enough to affect the 
classical Schwarzschild curvature ~ M~ 2 , i.e. when 



M \ M r 

where Mpi is the Planck mass y^hc/G ~ 2 x 10 -5 gm. Thus e is indeed very small for a stellar mass 
object, justifying the approximation a posteriori. With this semi- classical estimate for e we find 

I ~ VL^s - 3 x 10- 14 (J£-) 2 cm. (7.16) 



This is the same estimate that was obtained in (6.84), for where the vacuum polarization effects described 



by the stress tensor of the anomaly scalar fields becomes of the same order as the interior de Sitter 
curvature. Although still microscopic, the thickness of the shell is very much larger than the Planck 
scale Lpi ~ 2 x 1CP 33 cm. The energy density and pressure in the shell are of order M~ 2 and far below 
Planckian for M 3> Mpi, so that the geometry can be described reliably by Einstein's equations in both 
regions I and II. 

One may think of £ as the analog of the skin depth of a metal arising from its finite conductivity that 
cuts off the divergence in the Casimir stress tensor near a curved boundary [56., [57] . or of the healing 
length of an inhomogeneous Bose- Einstein condensate [148] , or finally as the thickness of a boundary 
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layer or stationary shock front in hydrodynamics |149| . In all of these examples from other areas of 
physics a lowest order macroscopic description of a bulk medium must be supplemented with some first 
order information about microscopic interactions in order to describe a rapid spatial crossover between 
two regions in which the bulk macroscopic description is completely adequate. In an EFT language 
this means that certain higher derivative interaction terms in the mean field eqs. related to fluctuations 
about the mean which are negligible in the bulk medium must be taken into account in the surface 
crossover layer. Although these higher derivative terms are present only because of some underlying 
microscopic degrees of freedom, the scale at which they become important are typically much larger 
than the fundamental microscopic or atomic scale (here Lpi), so that an EFT approach is still possible. 

Although f(r) is continuous across the interfaces at n and r 2 , in our simple model the discontinuity 
in the eqs. of state does lead to discontinuities in h{r) and the first derivative of f(r) in general. 
Defining the outwardly directed unit normal vector to the interfaces, n b = 5 r b ^h(r), and the extrinsic 
curvature K a b = V a n b , the Israel junction conditions |150| determine the surface stress energy r\ and 
surface tension a on the interfaces to be given by the discontinuities in the extrinsic curvature through 

~Vhdf~ 



2f dr 



4ttG(?7 - 2a) , 



(7.17a) 
(7.17b) 



Since h and its discontinuities are of order e, the energy density in the surfaces r\ ~ ea, while the 
surface tensions are of order e~ 2 . The simplest possibility for matching the regions is to require that 



the surface energy densities on each interface vanish. From (7.17b) this condition implies that h{r) is 



also continuous across the interfaces, which yields the relations, 

(l + «;i) 2 



h( n ) 

h(r 2 ) 

IM 

h(r 2 ) 



Hlr\ ~ e 
2GM 



T2 

w\ f(ri) 
w 2 /i(r 2 ) 



W 2 



C 



1 + Wl 
1 + ^2 



(7.18a) 
(7.18b) 
(7.18c) 



From (7.10) dw/dr < 0, so that u>2 < w\ and C < 1. In this case of vanishing surface energies r\ = 

(3 + u>i) 



the surface tensions are determined by (7.17) to be 

1 



a 2 ^ 



32vrG 2 M (1 + wi) 



€ 



1 



W2 



32ttG 2 M (1 + w 2 ) V e 



(7.19a) 
(7.19b) 
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to leading order in e at n and r<i respectively. The negative surface tension at the inner interface 
is equivalent to a positive tangential pressure, which implies an outwardly directed force on the thin 
shell from the repulsive vacuum within. The positive surface tension on the outer interfacial boundary 
coresponds to the more familiar case of an inwardly directed force exerted on the thin shell from without. 

The entropy of the configuration may be obtained from the Gibbs relation, p + p = sT + np, if the 
chemical potential p is known in each region. In the interior region I, p + p = and the excitations are 
the usual transverse gravitational waves of the Einstein theory in de Sitter space. Hence the chemical 
potential p may be taken to vanish and the interior has zero entropy density s = 0, consistent with a 
single macroscopic condensate state, S = k B InW(E) = for W(E) = 1. In region II there are several 
alternatives depending upon the nature of the fundamental excitations there. The p = p eq. of state 
may come from thermal excitations with negligible p or it may come from a conserved number density 
n of gravitational quanta at zero temperature. Let us consider the limiting case of vanishing p first. 

If the chemical potential can be neglected in region II, then the entropy of the shell is obtained from 
the eq. of state, p = p = (a 2 /8irG)(k B T/h) 2 . The T 2 temperature dependence follows from the Gibbs 
relation with p = 0, together with the local form of the first law dp = Tds. The Newtonian constant 
G has been introduced for dimensional reasons and a is a dimensionless constant. Using the Gibbs 
relation again the local specific entropy density s{r) = a 2 k 2 B T(r) j '4irh 2 G = a{k B /K){p/2 , KG) 1 i for local 
temperature T(r). Converting to our previous variable w, we find s = (ak B /4irhGr) and the entropy 
of the fluid within the shell is 

S = Att r s r 2 drh-^^ie l 2 hx(^) , (7.20) 
J n KG \w 2 J 



to leading order in e. Using (7.13) and (7.16), this is 



3 

$- ak B ^~W 57 ak B (^-V < S [:u . (7.21) 



The maximum entropy of the shell and therefore of the entire configuration is some 20 orders of magni- 



tude smaller than the Bekenstein-Hawking entropy (2.24) for a solar mass object, and of the same order 



3 

of magnitude as a typical progenitor of a few solar masses. The scaling of (7.21) with Ma is also in 



agreement with our general estimate (2.26) for a relativistic star, and the same order of magnitude as 



that for a supermassive star with M > 50 M & , whose pressure is dominated by radiation pressure 
Thus the formation of the gravastar from either a solar mass or supermassive stellar progenitor does 
not require an enormous generation or removal of entropy. Since the entropy is of the same order as 
that of a typical stellar progenitor, there is no information paradox and no significant entropy shedding 
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needed to produce a cold gravitational vacuum or 'grava(c)star' remnant. Because of the absence of 
an event horizon, the gravastar does not emit Hawking radiation. Since w is of order unity in the shell 
while r ~ r s , the local temperature of the fluid within the shell is of order T H ~ h/k B GM. The strongly 
redshifted temperature observed at infinity is of order y/eT H , which is very small indeed. Hence the rate 
of any thermal emission from the shell is negligible. There is no negative specific heat and no instability 
since the total bulk rest mass energy of the configuration remains essentially constant as the thin shell 
cools. 

If we do allow for a positive chemical potential within the shell, p > 0, then the temperature and 
entropy estimates just given become upper bounds, and it is possible to approach a zero temperature 
ground state with zero entropy. As the shell does cool, the entropy decreases very slowly from its initial 



value (7.21) and eq. of state of the shell must be replaced by a more accurate quantum stress tensor, 
not accounted for in this simple model. The non-singular final state of ultimate gravitational collapse 
is then a cold, completely dark object sustained against any further collapse solely by quantum zero 
point pressure of the interior. 

Realizing this alternative requires that a quantum gravitational vacuum phase transition intervene 
before the classical event horizon can form. This is exactly what the fluctuations of the anomaly scalar 
fields (f and ip described in Sec. [V] can provide. In a more realistic model based on the EFT of the 



anomaly (4.89), with (5.1)-(5.3), the rapid change of these fields near the horizon and their stress tensor 
can provide the suitable boundary layer which replaces the phenomenological eq. of state p = p in 



region II, and the somewhat artificial sharp interface boundaries with the surface stresses (7.19). The 



entire surface layer will still be of order £ in (7.16) in physical thickness, and if treated as very thin, 



be completely consistent with the classical formula (2.15) with an actual surface tension carried by the 
physical interface boundary. This is now under investigation. 

Incidentally, because of the necessarily non-vanishing surface tensions at the interface, p± ^ p there 
and the gravastar solution also explicitly evades the Buchdahl lower bound on the radius of a compact 
object R > 9GM/4, since this bound is derived under the assumption that the radial and tangential 
pressures are everywhere equal (and positive) inside the object [7] 1151} H52j . 

Since the exterior spacetime is Schwarzschild until distances of order of the diameter of an atomic 
nucleus from r = r s , a gravastar cannot be distinguished from a black hole by present observations of X- 
ray bursts |153j . However, the shell with its maximally stiff eq. of state p = p, where the speed of sound 
is equal to the speed of light, could be expected to produce explosive outgoing shock fronts in the process 
of formation. Active dynamics of the shell may produce other effects that would distinguish gravastars 
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from black holes observationally, possibly providing a more efficient particle accelerator and central 
engine for energetic astrophysical sources. The spectrum of gravitational radiation from a gravastar 
should bear the imprint of its fundamental frequencies of vibration, and hence also be quite different 
from a classical black hole. 

The interior de Sitter region with p = —p may be interpreted also as a cosmological spacetime, with 
the horizon of the expanding universe replaced by a quantum phase interface. The possibility that the 
value of the vacuum energy density in the effective low energy theory can depend dynamically on the 
state of a gravitational condensate may provide a new paradigm for cosmological dark energy in the 
universe. The proposal that other parameters in the standard model of particle physics may depend 
on the vacuum energy density within a gravastar has been discussed by Bjorken [154J. A stable bubble 
of positive vacuum zero point energy also realizes Dirac's idea for an extensible model of the electron 
[155|, I156j , but in the case of the attractive gravitational force instead (so that the repulsive quantum 
zero point force can be balanced by rather than added to the classical self-force of the extended body) . 
Similar suggestions for removing the singularity in gravitational collapse were also made by Sakharov 
and Gliner [157J. 

In the original paper on gravastars [28], the stability of the configuration was studied in the same 



hydrodynamic approximation used to construct it. In each of the three regions (7.4), fluctuations are 
completely stable, for the clear physical reasons that the de Sitter and Schwarzschild geometries are 
stable to fluctuations and the fluid making up the this shell also has a physical eq. of state. In this 
hydrodynamic model the positions of the interfaces at r\ and r-i are fixed by the externally provided 



information from the estimate (6.84). Thus the exact positions of r\ and r2 are not fixed dynamically, 
and there are zero modes of neutral stability corresponding to shifts of r\ and r<i [28, 29J. This can only 
be remedied by a fully dynamical theory in which the gravastar solution is obtained by the variation of 
a well-defined action functional. This we did not possess in 2001, but do have now have in the effective 



theory of the Einstein- Hilbert term plus the covariant action functional (4.63) with (4.64) generated 
by the anomaly. The subsequent studies of the solutions of the equations of motion of the anomaly 



scalar field (p and ip in both the Schwarzschild and de Sitter geometries |115j . reviewed in Sec. IV, show 



that the stress tensor of the anomaly scalars (5.2)- (5.3) can become large near both the black hole and 
cosmological horizons and provide exactly the dynamical stresses that were hypothesized in [28} 129] . 
Work is now currently in progress to find a gravastar solution to the effective field equations including 
the anomaly terms. As an extremum of a well-defined action principle, the dynamical stability of the 
solution can then be studied without any additional assumptions. 
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Since the appearance of the papers [2H1 ES] , a number of other authors have addressed the issue 
of gravastar stability in a number of different ways [158-166], all requiring some phenomenological 
assumptions and parameters. Not surprisingly, the results obtained depend upon the values of these 
parameters, although generally a wide class of stable gravastar-like solutions have been found. The 
authors of Ref. [164J claimed to have found a more generic instability of the ergoregion of rotating 
non-black hole compact objects. However it was shown in Ref. |165j that this conclusion can be avoided 
by making somewhat less restrictive assumptions than those of [164J. 

More to the point, it is by no means certain that a rotating gravastar possesses an ergoregion at all. 
Our discussion of the stress tensor dependence upon the norm of the timelike Killing field (K a K a )z in 
Sec. |TI]and the global aspects of the trace anomaly in Sees. |IV] and [V] lead to the expectation that large 
effects in the expectation value of the stress tensor are to be expected where K a first becomes null. In 
the non-rotating Schwarzschild geometry this is the horizon at r = r s . However, in the Kerr geometry 
the unique timelike Killing vector at large r first becomes null at the outer boundary of the ergosphere. 
Hence if the region of large quantum effects is here, the interior ergosphere region will be changed and 
likely not exist at all, removing the question of stability of this region. 

Clearly this question and the entire subject of rotating gravastars, their possible magnetic fields and 
detailed interactions with matter, relevant for realistic collapsed stars and astrophysical observations 
require further study before definite answers can be given or predictions made. It has been pointed out 
in Ref. |153j that astrophysical observations or non-observations of X-ray bursts or other emisions from 
compact objects are not definitive for or against the existence of a physical surface. The interpretations 
of this data also depend heavily upon models of any surface and its possible interactions with matter. An 
interesting possibility for detecting emission from the surface based on a model of these interactions has 
been given in |167j . A more model independent prediction of a physical surface for a compact object is 
that it will have normal modes of vibration, and a characteristic discrete energy spectrum which should 
be discernible in gravitational wave observations. For this reason it is important to develop both the 
theory and the observational techniques for detecting a gravitational wave spectrum that may differ 
from that expected from pure classical GR |168j . 

VIII. SUMMARY AND OUTLOOK 

Although quantum effects can often be quite subtle, and it is often assumed that in gravity they 
play a role only at the Planck scale, the challenges presented to current theory by quantum effects in 
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black hole and cosmological spacetimes suggest otherwise. There are two macroscopic systems where 
quantum effects may in fact be crucial, namely the final state of gravitational collapse and in accounting 
for cosmological dark energy of the universe itself. These are the principal challenges presented to the 
reconcilation of General Relativity with quantum mechanics on macroscopic scales. 

In this article we have reviewed the status of the horizons in the classical Schwarzschild and de 
Sitter solutions of Einstein's equations. The paradoxes of black hole physics in particular argue for a 
careful analysis of the possible importance of large backreaction effects of quantum vacuum polarization 
in the vicinity of the event horizon. Large quantum effects of this kind do not violate the Principle 
of Equivalence, but simply serve to emphasize that quantum coherence and entanglement effects can 



depend on gauge invariant but non-local integrals of gauge potentials such as (2.46) in electrodynamics 



or (2.48) and (2.51) in gravity. Physically this is because quantum matter has wavelike properties and 
cannot be absolutely localized to a point. Mathematically, it is expressed by the fact that gauge theories 
including gravity are theories of gauge connections with generally non-trivial fiber bundle structure. 
Thus singular coordinate frame transformations of the kind often considered in classical GR ignore the 
possibility of the new degrees of freedom associated with these improper "gauge" transformations, and 
may be unwarranted in quantum theory. 

Based on a number of examples from both flat space QED and QCD, and gravity in both 2 and 
4 dimensions, the central role of anomalies was discussed. Anomalies play a special role in EFT's 
since they provide an essentially unique way to violate the naive decoupling of ultraviolet from infrared 
sectors of a Lorentz invariant or generally covariant theory. This is because they are at the same time 
local in the divergence of the classically conserved current, but also imply a non-local structure in the 
effective action. The non-locality is associated with the propagation of massless correlated pair states, 
analogous to Cooper pairs in a superconductor or collective modes in other media. Both the chiral 
anomalies of QED and QCD and the trace anomaly in a gravitational background lead generically to 
additional propagating degrees of freedom than is apparent from the original classical or single particle 
Lagrangian. The spectral density in these propagating anomaly channels obey ultraviolet sum rules 
and hence remain physical states, even away from the strictly conformal on-shell limit. The states in 
the trace anomaly channel couples only to fourth order curvature invariants and therefore only very 
weakly to matter, so that they would not have been detected so far in terrestrial or solar system tests 
of gravity. 

The + states in trace anomaly amplitudes can be described as local scalar fields in an effective field 
theory description. The logarithmic scaling with distance of the effective action of the trace anomaly 
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(4.63)-(4.64) implies that it remains a (marginally) relevant operator in the EFT of low energy gravity, 
and should be explicitly appended to the classical Einstein-Hilbert terms to take account of quantum 
effects of pair correlations due to the trace anomaly. This has the consequence that General Relativity 
receives quantum corrections relevant at macroscopic distances. 

Due to the topological character of the Euler-Gauss-Bonnet invariant in the trace anomaly, the scalar 
degrees of freedom in the anomaly effective action are sensitive to non-local or global aspects of the 
underlying quantum theory, and this is the fundamental reason why they can play a decisive role near 
classical event horizons, where they have macroscopically large effects. As the horizon is approached, 
the large blueshift of frequencies overwhelms all finite mass and energy scales, so that all fields become 
essentially massless, and the behavior is conformal, accurately described by the conformal anomaly 
effective action. This has been shown quantitatively by comparison of the stress tensor of the anomaly 



(5.1)-(5.3) in terms of the two scalar fields ip and ip to the renormalized expectation value of the stress 
tensor of quantum matter fields and their vacuum polarization effects computed by standard methods. 
The important qualitative conclusion of the study of the anomaly stress tensor in the Schwarzschild and 
de Sitter cases is that the diverging behavior proportional to (—K a K a )~ 2 = (—gtt)~ 2 is quite generic 
and appears in a wide range of states other than the Boulware state. A fine tuning of the state is 
necessary to remove all divergences of the stress tensor on the horizon in each case, which does not seem 
to be a priori warranted. 

Since the anomaly scalar fields have kinetic terms, they must also be treated as dynamical fields in 
their own right. In fact, their fluctuations induce the running of the classical cosmological term, and 
provide a mechanism for the vacuum energy to be dynamical, and dependent upon infrared effects, 
rather than a constant. The value GA = is an infrared fixed point of the renormalization group flow. 
Thus a dynamical relaxation of the vacuum energy to zero is possible by this mechanism of anomaly 
scalar field fluctuations [681 fTLH [T69] . 

The restoration of conformal invariance these fluctuations imply may be observable in the Cosmic 
Microwave Background. One form this conformal invariance could show up is in non-Gaussian cor- 



relations of a specific angular form (6.27), assuming that conformal invariance is realized on the flat 
spatial sections of a spatially hologeneous, isotropic FLRW cosmology. Another way of realizing con- 
formal invariance in de Sitter space is associated with the breakdown of the 0(4, 1) isometry group to 
0(3), relevant for cosmological dark energy in the present epoch. These different realizations lead to 
different cosmological models and different predictions for the non-Gaussianities which may allow us to 
differentiate models. 
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In de Sitter space the fluctuations of the anomaly scalars imply scalar degrees of freedom from 
fundamental theory without the ad hoc introduction of an inflaton field. Their fluctuations grow large 
at the de Sitter horizon as at the Schwarzschild horizon. This suggests that the horizon of a black hole 
should be viewed instead as the locus of a quantum phase transition induced by the anomaly scalar 
degrees of freedom, where the effective value of the gravitational vacuum energy density can change. In 
the EFT including the trace anomaly terms, A e ff becomes a dynamical condensate whose value depends 
upon the macroscopic boundary conditions at the horizon. By taking a positive value in the interior, the 
effective "repulsive" dark energy core can remove the black hole singularity of the classical Schwarzschild 
solution, replacing it with an horizon volume of de Sitter space. The resulting gravitational condensate 
star or gravastar model resolves all black hole paradoxes, and provides a testable non-singular quantum 
alternative to black holes as the final state of complete gravitational collapse. 

The cosmological dark energy of our universe may be a macroscopic finite volume condensate whose 
value depends not on microphysics but on the cosmological horizon scale. Finally, both the sensitivity 
of the trace anomaly terms to microphysics and the analogies with many body and condensed matter 
systems suggest that further development of the ideas presented here may lead to an entirely new basis 
for the microscopic degrees of freedom of quantum gravity and therefore of the constituents of spacetime 
itself [T70] . 



Acknowledgments 

I am very much indebted to all my collaborators, whose work is reviewed in this article: Paul R. 
Anderson, Ignatios Antoniadis, Maurizio Giannotti, Carmen Molina-Paris, Ruslan Vaulin, and especially 
Pawel O. Mazur, with whom the hypothesis of gravitational condensate stars was developed. I wish 
also to thank Michal Praszalowicz and the other organizers of the XLIX Krakow School of Theoretical 
Physics for the invitation to lecture at that school. Notes for these lectures eventually grew into this 
article. Finally I wish to thank the I. Antoniadis and the CERN Theoretical Physics Group for the 
Scientific Associateship and hospitality at CERN, Oct., 2009-Apr., 2010 where and when most of this 
article was written. 



120 



References 



[1] M. Planck, Ann. Physik 4, 553 (1901). 
[2] A. Einstein, Ann. Physik 22, 180 (1907). 
[3] R. C. Tolman, Phys. Rev. 55, 364 (1939); 

J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55. 374 (1939); 

B. K. Harrison, K. S. Thorne, M. Wakano, and J. A. Wheeler, Gravitation Theory and Gravitational 
Collapse (Univ. of Chicago Press, Chicago 1965); 

For more recent results sec e.g. I. Bombaci, Astron. Astrophys. 305, 871 (1996). 
[4] K. Schwarzschild, Sitzber. Preuss. Akad. Wiss. Berlin, Math.-Phys., 189 (1916). 
[5] A. Einstein, Ann. Math. 40, 922 (1939). 

[6] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, San Francisco 1973). 
[7] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity 
(Wiley, New York 1972). 

[8] C. Will, Theory and Experiment in Gravitational Physics (Cambridge Univ. Press, Cambridge 1993); Liv. 

Rev. Rel. 4, 4 (2001). 
[9] M. D. Kruskal, Phys. Rev. 119, 1743 (1960); 
G. Szekeres, Publ. Mat. Debrecen. 7, 285 (1960). 

[10] P. A. M. Dirac, Proc. Roy. Soc. Lon. A270, 354 (1962). 

[11] H. Rcissner, Ann. Phys. (Ger.) 50, 106 (1916); 

G. Nordstr0m, Proc. Kon. Ned. Akad., 20, 1238 (1918). 

[12] R. P. Kerr, Phys. Rev. Lett. 11, 237 (1963). 

[13] R. H. Boycr and R. W. Lindquist, J. Math. Phys. 8, 265 (1967); 
B. Carter, Phys. Rev. 174, 1559 (1968). 

[14] S. W. Hawking & G. F. R. Ellis, The Large Scale Structure of Space-Time (Cambridge Univ. Press, Cam- 
bridge 1973). 

[15] D. Christodoulou, Phys. Rev. Lett. 25, 1596 (1970); 

D. Christodoulou and R. Ruffini, Phys. Rev. D 4, 3552 (1971). 
[16] S. W. Hawking, Phys. Rev. Lett. 26, 1344 (1971). 

[17] J. M. Bardecn, B. Carter, and S. W. Hawking, Comm. Math. Phys. 31, 161 (1973). 
[18] L. Smarr, Phys. Rev. Lett. 30, 71 (1973). 

[19] J. D. Bckcnstcin, N. Cimento Lett. 4, 737 (1972); Phys. Rev. D7, 2333 (1973); ibid. D9, 3292 (1974). 
[20] S. W. Hawking, Nature 248, 30 (1974); Comm. Math. Phys. 43, 199 (1975). 
[21] J. Schwinger, Phys. Rev. 82, 664 (1951); ibid. 93, 615 (1954). 
[22] S. W. Hawking, Phys. Rev. D 13, 191 (1976). 

[23] D. J. Gross, M. J. Perry, and L. G. Yaffe, Phys. Rev. D 25, 330 (1982); 



121 



J. W. York, Phys. Rev. D 33, 2092 (1986). 

G. 't Hooft, Nucl. Phys. B256, 727 (1985); Nucl. Phys. Proc. Suppl. 68, 174 (1998). 

Ya. B. Zel'dovich and I. D. Novikov, Relativistic Astrophysics, Vol. 1 (Univ. of Chicago Press, Chicago 
1971) [Stars and Relativity (Dover, Mineola, NY 1996)]. 
S. W. Hawking, Comm. Math. Phys. 87, 395 (1982). 
S. W. Hawking, Phys. Rev. D 72, 084013 (2005). 



P. O. Mazur and E. Mottola, e-print arXiv: gr-qc/0109035 



P. O. Mazur and E. Mottola, Proc. Nat. Acad. Sci. 101, 9545 (2004). 

S. Carlip, Gen. Rel. Grav. 39, 1519 (2007); Int. J. Mod. Phys. D 17, 659 (2008), and references therein. 
L. Bombelli, R. K. Koul, J. Lee, and R. Sorkin, Phys. Rev. D 34, 373 (1986). 

V. P. Frolov and D. V. Fursaev, Phys. Rev. D 56, 2212 (1997); Class. Quant. Grav. 15. 2041 (1998). 
W. Israel, "Black Hole Thermodynamics," in Current Trends in Relativistic Astrophysics, L. Fernandez- 
Jambrina and L. M. Gonzazez-Romero, eds. (Springer, Berlin 2003); Led. Notes Phys. 617, 15 (2003). 
B. S. DeWitt, Phys. Rep. 19, 295 (1975). 

D. G. Boulware Phys. Rev. D 11, 1404 (1975); ibid. 13, 2169 (1976). 
S. M. Christcnsen and S. A. Fulling, Phys. Rev. D 15, 2088 (1977). 

R. C. Tolman, Phys. Rev. 35, 904 (1930); Relativity, Thermodynamics and Cosmology (Clarendon, Oxford 

1934) (Dover, Mineola, N. Y. 1987). 

W. G. Unruh, Phys. Rev. D 14 870 (1976). 

J. B. Hartle and S. W. Hawking, Phys. Rev. D13 2188 (1976). 

W. Israel, Phys. Lett. A57 107 (1976). 

S. W. Hawking, "The Path Integral Approach to Quantum Gravity," in General Relativity: An Einstein 

Centenary Survey, S. W. Hawking and W. Israel, eds. (Cambridge Univ. Press, Cambridge 1979). 

G. W. Gibbons and M. J. Perry, Phys. Rev. Lett. 36, 985 (1976); 

G. W. Gibbons and M. J. Perry, Proc. R. Soc. Lond. A358, 467 (1978). 

S. A. Fulling, J. Phys. A10, 917 (1977). 

G. W. Gibbons, "Quantum Fields in Curved Spacetime," in Ref. [JT] 

P. R. Anderson, W. A. Hiscock, and D. A. Samuel, Phys. Rev. D 51, 4337 (1995). 

R. Kubo, Jour. Phys. Soc. Japan 12, 570 (1957); 

P. C. Martin and J. Schwinger, Phys. Rev. 115, 1342 (1959). 

E. Mottola, Phys. Rev. D 33, 1616 (1986). 

F. London and H. London, Proc. R. Soc. Lon. A 149, 71 (1935); Physica 2, 341 (1935). 
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 108 1175 (1957). 

L. P. Gorkov, Sov. Phys. JETP 9, 1364 (1959); 

A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinskii, Methods of Quantum Field Theory in Statistical 
Physics (Prentice-Hall, Englewood Cliffs, NJ 1963) (Dover, New York 1975). 



122 



[51] M. F. Atiyah and R. Bott, Bull. Amer. Math. Soc. 72, 245 (1966); 

R. Bott, Mich. Math. J. 14, 231 (1967). 
[52] Y. Aharonov and D. Bohm, Phys. Rev. 115 485 (1959). 
[53] S. S. Chern, Ann. Math. 45, 747 (1944); ibid. 46, 674 (1945). 
[54] W. Pauli (unpublished); 

C. P. Enz and A. Thellung, Helv. Phys. Acta 33, 839 (1960); 

S. E. Rugh and H. Zinkernagel, Stud. Hist. Philos. Mod. Phys. 33, 663 (2002); 

N. Straumann, e-print arxiv: gr-qc/0208027 
[55] H. B. G. Casimir, Proc. Kon. Ned. Acad. Wetenschap. 51, 793 (1948); 

M. J. Sparnaay, Physica 24, 751 (1958). 
[56] D. Deutsch and P. Candelas, Phys. Rev. D 20, 3063 (1979). 
[57] M. Bordag, U. Mohideen, and V. M. Mostepanenko, Phys. Rept. 353, 1 (2001); 

K. A. Milton, The Casimir Effect: Physical Manifestations of Zero-Point Energy (World Scientific, Singa- 
pore 2001); 

M. Bordag, Advances in the Casimir Effect (Oxford Univ. Press, Oxford 2009). 
[58] A. G. Ricss et. al, Astron. J. 116, 1009 (1998); ibid. 607, 665 (2004); 

S. Perlmutter et. al, Astrophys. J. 517, 565 (1999); 

J. L. Tonry et. al, Astrophys. J. 594, 1 (2003). 
[59] W. de Sitter, Mon. Not. R. Astron. Soc 78, 3 (1917); 

A. S. Eddington, The Mathematical Theory of Relativity (Cambridge Univ. Press, Cambridge 1924). 
[60] E. Mottola, Phys. Rev. D 31, 754 (1985). 
[61] E. Mottola, Phys. Rev. D 33, 1616 (1986). 
[62] P. O. Mazur and E. Mottola, Nucl. Phys. B278 (1986). 

[63] E. Mottola, in Physical Origins of Time Asymmetry, J. J. Halliwell et al. eds. (Cambridge Univ. Press, 
Cambridge 1993). 

[64] N. C. Tsamis and R. P. Woodard, Nucl. Phys. B474, 235 (1996); Ann. Phys. 253, 1 (1997); ibid. 267, 145 
(1998). 

[65] V. Mukhanov, L. R. W. Abramo, and R. Brandenberger, Phys. Rev. Lett. 78, 1624 (1997); 

L. R. Abramo, R. H. Brandenberger, and V. F. Mukhanov, Phys. Rev. D 56, 3248 (1997). 
[66] L. R. Abramo and R. P. Woodard, Phys. Rev. D 60, 044010 (1999); ibid. 65, 063516 (2002), and references 

therein. 

[67] A. Ishibashi and R. W. Wald, Class. Quant. Grav. 23, 235 (2006). 
[68] I. Antoniadis, P. O. Mazur and E. Mottola, New J. Phys. 9, 11 (2007). 

[69] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, aand E. Mottola, Phys. Rev. Lett. 67, 2427 (1991); 
Phys. Rev. D 45, 4659 (1992); 

Y. Kluger, E. Mottola, and J. M. Eisenberg, Phys. Rev. D 58, 125015 (1998). 



123 



[70] S. Habib, C. Molina-Pan's, and E. Mottola, Phys. Rev. D 61, 024010 (2000). 
[71] O. Nachtman, Comm.. Math. Phys. 6, 1 (1967); 

N. A. Chernikov and E. A. Tagirov, Ann. Poincare Phys. Theor. A9, 109 (1968) 

T. S. Bunch and P.C.W. Davies, Proc. Roy. Soc. Lond. A360, 117 (1978). 
[72] L. H. Ford and L. Parker, Phys. Rev. D 16, 245 (1977); 

B. Allen, Phys. Rev. D 32, 3136 (1985); 

B. Allen and A. Folacci, Phys. Rev. D 35, 3771 (1987). 
[73] B. L. Beers and H. H. Nickle, Lett. N. Cimento 4, 320 (1972). 
[74] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2752 (1977). 

[75] I. Antoniadis and E. Mottola, preprint CERN-TH-4605/86; Jour. Math. Phys. 32, 1037 (1991). 
[76] B. Allen and M. Turyn, Nucl. Phys. B292, 813 (1987). 

[77] E. G. Floratos, J. Iliopoulos, and T. N. Tomaras, Phys. Lett. B197, 373 (1987). 
[78] N. P. Mhyrvold, Phys. Rev. D 28, 2439 (1983). 
[79] L. Ford, Phys. Rev. D 31, 710 (1985). 

[80] A. M. Polyakov, Nucl. Phys. B797, 199 (2008); ibid. B834, 316 (2010). 

[81] S. L. Adler, Rev. 177, 2426 (1969); Lectures on Elementary Particle and Quantum Field Theory, Vol. 1, S. 

Deser, M. Grisaru, and H. Pendleton cds. (MIT Press, Cambridge 1970). 
[82] J. S. Bell and R. Jackiw, N. Cimento A 60, 47 (1969). 

[83] R. Jackiw, in Lectures on Current Algebra and Its Applications, S. B. Treiman, R. Jackiw, and D. J. Gross, 

eds. (Princeton Univ. Press, Princeton 1972). 
[84] A. D. Dolgov and V. I. Zakharov, Nucl. Phys. B27, 525 (1971). 

[85] J. Hofejsf, Czech. J. Phys. B 35, 820 (1985); Phys. Rev. D 32, 1029 (1985); J. Phys. G: Nucl. Phys. 12, 

L7 (1986); Czech. J. Phys. 42, 345 (1992). 
[86] L. Rosenberg, Phys. Rev. 129, 2786 (1963). 
[87] J. Hofejsf, Czech. J. Phys. 42, 241 (1992). 

[88] M. Giannotti and E. Mottola, Phys. Rev. D 79, 045014 (2009). 

[89] See e.g. R. A. Bertlmann, Anomalies in Quantum Field Theory (Clarendon Press, Oxford 1996). 
[90] A. Smailagic and E. Spallucci, Phys. Rev. D 62, 045010 (2000). 
[91] L. S. Brown, N. Cimento, 29, 619 (1963). 

[92] S. R. Coleman and B. Grossman, Nucl. Phys. B203, 205 (1982). 

[93] C. Adam, R. A. Bertlmann, and P. Hofer, Riv. N. Cimento, 16N8, 1 (1993); 

R. A. Bertlmann and E. Kohlprath, Ann. Phys. 288, 137 (2001). 
[94] W. A. Bardeen, H. Fritzsch and M. Gell-Mann, in Scale and Conformal Symmetry in Hadron Physics, R. 

Gatto, ed. (Wiley, New York 1973); 

T. P. Cheng and L. F. Li, Gauge Theory of Elementary Particle Physics, (Oxford Univ. Press, Oxford 
1984). 



124 



[95] S. Weinberg, The Quantum Theory of Fields, Vol. 2, (Cambridge Univ. Press, Cambridge 1995). 

[96] G. 't Hooft, Recent Developments in Gauge Theories, 1979 NATO Cargese Advanced Study Institutes, 

Series B, Physics, vol. 59 (Plenum, New York 1980); See also Ref. 95], Sec. 22.5, pp. 389-396. 
[97] G. Veneziano, Nucl. Phys. B159, 213 (1979); 
P. Di Vecchia, Phys. Lett. 85B, 357 (1979); 
P. Di Vecchia and G. Veneziano, Nucl. Phys. B171, 253 (1980) 
[98] A. Aurilia, Y. Takahashi, and P. K. Townsend, Phys. Lett. 95B, 265 (1980). 
[99] M. S. Chanowitz and J. R. Ellis, Phys. Rev. D 7, 2490 (1973). 
[100] S. L. Adler, J. C. Collins and A. Duncan, Phys. Rev. D 15, 1712 (1977). 
[101] I. T. Drummond and S. J. Hathrell, Phys. Rev. D 21, 958 (1980). 
[102] See e.g. P. Achard et. at, L3 Collaboration, Phys. Lett. B623, 26 (2005). 

[103] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge Univ. Press, Cambridge, 

1982), and references therein. 
[104] J. Wess and B. Zumino, Phys. Lett. B 37, 95 (1971); 

E. Witten, Nucl. Phys. B156 269 (1979); Nucl. Phys. B B223 422 (1983); 
W. A. Bardeen and B. Zumino, Nucl. Phys. B244, 421 (1984). 

[105] A. M. Polyakov, Phys. Lett. B103, 207 (1981). 

[106] I. Antoniadis, P. O. Mazur, and E. Mottola, Nucl. Phys. B388 627 (1993). 
[107] P. O. Mazur and E. Mottola, Phys. Rev. D 64, 104022 (2001). 

[108] V. G. Knizhnik, A. M. Polyakov, and A. B. Zamolodchikov, Mod. Phys. Lett. A 3 819 (1988); 

F. David, Mod. Phys. Lett. A 3, 1651 (1988); 

J. Distler and H. Kawai, Nucl. Phys. B 321, 509 (1989). 
[109] See e.g. J. Ambj0rn, B. Durhuus, and T. Jonsson, Quantum Geometry: A Statistical Field Theory Approach 
(Cambridge Univ. Press, Cambridge 1997); 

J. Ambj0rn, M. Carfora, and A. Marzuoli, The Geometry of Dynamical Triangulations (Springer, Berlin 
1997). 

[110] S. Cattcrall and E. Mottola, Phys. Lett. B467, 29 (1999); Nucl. Phys. Proc. Suppl. 83, 748 (2000). 
[Ill] D. Capper and M. J. Duff, N. Cimento A23 173 (1974); Phys. Lett. A53 361 (1975); 
S. Deser, M. Duff and C. J. Isham, Nucl. Phys. Bill 45 (1976); 

M. Duff, Nucl. Phys. B125 334 (1977); Nucl. Phys. B234 269 91984); Class. Quant. Grav. 11, 1387 (1994). 
[112] A. A. Starobinsky, Phys. Lett. B91, 99 (1980); 
A. Vilenkin, Phys. Rev. D 32, 2511 (1985); 

S. W. Hawking, T. Hertog, and H. Reall, Phys. Rev. D 63, 083504 (2001). 
[113] R. J. Riegert, Phys. Lett. B134, 56 (1984); 

E. S. Fradkin and A. A. Tseytlin, Phys. Lett. B134, 187 (1984). 
[114] I. Antoniadis and E. Mottola, Phys. Rev. D 45, 2013 (1992). 



125 



115] E. Mottola and R. Vaulin, Phys. Rev. D 74, 064004 (2006). 

116] R. P. Feynman, F. B. Morinigo, and W. G. Wagner, Feynman Lectures on Gravitation (Univ. of Bangalore 

Press, New Delhi 1995). 
117] P. Jordan, Z. Phys. 157, 112 (1959); 

C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961). 
118] J. F. Donoghue, E. Golowich, and B. R. Holstein, Dynamics of the Standard Model (Cambridge Univ. Press, 

Cambridge 1992 ). 
119] J. F. Donoghue, Phys. Rev. D 50, 3874 (1994). 
120] K. G. Wilson and J. Kogut, Phys. Rep. C12, 75 (1974); 

K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 
121] B. P. Jensen, J. G. Mc Laughlin, and A. C. Ottewill, Phys. Rev. D 45, 3002 (1992). 
122] M. R. Brown and A. C. Ottewill, Phys. Rev. D 31, 2514 (1985); ibid. 34, 1776 (1986). 
123] D. N. Page, Phys. Rev. D 25, 1499 (1982). 

124] M. R. Brown, A. C. Ottewill, and D. N. Page, Phys. Rev. D 33, 2840 (1986). 

125] R. Balbinot, A. Fabbri, and I. L. Shapiro, Phys. Rev. Lett. 83, 1494 (1999); Nucl. Phys. B559, 301 (1999). 
126] P. R. Anderson, E. Mottola, and R. Vaulin, Phys. Rev. D 76, 124028 (2007). 
127] L. S. Brown and J. P. Cassidy, Phys. Rev. D 16, 1712 (1977). 
128] S. D. Odintsov, Z. Phys. C 54, 531 (1992); 

I. Antoniadis and S. D. Odintsov, Phys. Lett. B 343, 76 (1995). 
129] I. Antoniadis, P. O. Mazur, and E. Mottola, Phys. Lett. B444 284 (1998). 
130] I. Antoniadis, P. O. Mazur, and E. Mottola, Phys. Rev. D 55, 4756; ibid. 55, 4770 (1997). 



131] I. Antoniadis, P. O. Mazur, and E. Mottola, e-print arXiv: astro-ph/9705200 



132] V. F. Mukhanov, H. A. Feldman, and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 

133] See e.g. A. R. Liddle and D. H. Lyth, Cosmological Inflation and Large-Scale Structure (Cambridge Univ. 

Press, Cambridge 2000). 

134] N. Bartolo, E. Komatsu, S. Matarrese, and A. Riotto, Phys. Rept. 402, 103 (2004). 

135] E. R. Harrison, Phys. Rev. D 1, 2726 (1970); 

Ya. B. Zel'dovich, Mon. Not. R. Astron. Soc. 160, IP (1972). 

136] I. Antoniadis, P. O. Mazur, and E. Mottola, Phys. Rev. Lett. 79, 14 (1997). 

137] A. M. Polyakov, Sov. Phys. JETP Lett. 12 381 (1970). 

138] P. J. E. Peebles, Principles of Physical Cosmology (Princeton Univ. Press, Princeton 1993). 

139] P. R. Anderson, C. Molina-Paris, and E. Mottola, Phys. Rev. D 67, 024026 (2003); ibid. 72, 043515 (2005). 

140] P. R. Anderson, C. Molina-Paris, and E. Mottola, Phys. Rev. D 80, 084005 (2009). 

141] G. T. Horowitz and R. M. Wald, Phys. Rev. D 21, 1462 (1980); ibid. 25, 3408 (1982). 

142] A. A. Starobinsky, Piz. Eksp. Teor. Fiz. 34, 460 (1981) [JETP Lett. 34, 438 (1981)]. 

143] J. M. Bardecn, Phys. Rev. D 22, 1882 (1980). 



126 



144] J. M. Stewart, Class. Quant. Grav. 7, 1169 (1990). 

145] P. O. Mazur and E. Mottola, in Quantum field theory under the influence of external conditions QFEXT03, 

K. A. Milton, ed. (Rinton, Princeton 2004). 
146] G. Chaplinc, E. Hohlfield, R. B. Laughlin, and D. Santiago, Phil. Mag. B81, 235 (2001); Int. Jour. Mod. 

Phys. A18, 3587 (2003). 
147] R. B. Laughlin, Int. Jour. Mod. Phys. A18, 831 (2003). 
148] E. P. Gross, J. Math. Phys. 4, 195 (1963). 

149] J. D. Anderson, Fundamentals of Aerodynamics, 3 rd ed., (McGraw-Hill, New York 2001); Physics Today 
58, 42 (2005). 

150] W. Israel, N. Cimento B44, 1; ibid. B48, 463 (1966). 
151] H. A. Buchdahl, Phys. Rev. 116, 1027 (1959). 

152] C. Cattoen, T. Faber, and M. Visser, Class. Quant. Grav. 22, 4189 (2005). 

153] M. A. Abramowicz, W. Kluzniak, and J.-P. Lasota, Astron. Astroph., 396, L31 (2002). 

154] J. D. Bjorken, Phys. Rev. D 67, 043508 (2003). 

155] P. A. M. Dirac, Proc. Roy. Soc. Lond. A268, 57 (1962). 

156] T. H. Boyer, Phys. Rev. 174, 1764 (1968). 

157] A. D. Sakharov, Sov. Phys. JETP 22, 241 (1966); 

E. B. Glincr, Sov. Phys. JETP 22, 378 (1966). 
158] M. Visser and D. L. Wiltshire, Class. Quant. Grav. 21, 1135 (2004). 
159] I. Dymnikova and E. Galaktionov, Class. Quant. Grav. 22, 2331 (2005). 
160] B. Carter, Class. Quant. Grav. 22, 4551 (2005) 
161] F. S. N. Lobo, Class. Quant. Grav. 23, 1525 (2006). 

162] A. DeBenedictis, D. Horvat, S. Ilijic, S. Kloster, and K. S. Viswanathan, Class. Quant. Grav. 23, 2303 
(2006). 

163] C. B. M. H. Chirenti and L. Rezzolla, Class. Quant. Grav. 24, 4191 (2007). 

164] V. Cardoso, P. Pani, M. Cadoni, and M. Cavaglia, Phys. Rev. D 77, 124044 (2008); Class. Quant. Grav. 
25, 195010 (2008). 

165] C. B. M. H. Chirenti and L. Rezzolla, Phys. Rev. D 78, 084011 (2008). 

166] P. Rocha, R. Chan, M. F. A. da Silva, and A. Wang, J. Cosmo. Astrop. Phys. 11, 010 (2008). 



167] J. Barbieri and G. Chapline, Phys. Lett. B590, 8 (2004); e-print arXiv:0806.1550 and arViv:0812.0833. 
168] N. Yunes and F. Pretorms, Phys. Rev. D 80, 122003 (2009). 

169] E. Mottola, e-Print arXiv:1006.356, in Quantum field theory under the influence of external conditions 

QFEXT09 , K. A. Milton and M. Bordag, eds. (World Scientific, Singapore 2010). 
170] P. O. Mazur, Acta Phys. Polon. B 27 1849 (1996); e-print arXiv: |hep-th/9712208"l 



127 



